(Leave  blank) 


AD 


Award  Number:  DAMD17-03-1-0363 


TITLE:  Therapeutic  vascular  targeting  and  irradiation: 

correlation  of  MRI  tissue  changes  at  cellular  and 
mol  ecul  ar  levels  to  optimizing  outcome 


PRINCIPAL  INVESTIGATOR:  Dawen  Zhao,  M.D.,  Ph.D. 


CONTRACTING  ORGANIZATION:  University  of  Texas 

Southwestern  Medical  Center  at  Dal  I  as 
Dal  I  as,  TX  753  90-  9058 


REPORT  DATE:  June  2008 


TYPE  OF  REPORT:  Final 


PREPARED  FOR:  U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


DISTRIBUTION  STATEMENT:  (Check  one) 

X  Approved  for  public  release;  distribution  unlimited 

Distribution  limited  to  U.S.  Government  agencies  only; 
report  contains  proprietary  information 

The  views^  opinions  and/or  findings  contained  in  this  report  are 
those  of  the  author  (s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position^  policy  or  decision  unless  so 
designated  by  other  documentation. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  074-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of 
Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503 

1.  Agency  Use  Only  (Leave  blank) 

2.  Report  Date 

30-06- 2008 

3.  Report  Type  and  Period  Covered  (i.e.,  annual  1  Jun  00  -  31  May  01) 

Final,  1  j  UN  2003  -  31  MAY  2008 

4.  Title  and  Subtitle 

5.  Award  Number 

Therapeutic  vascular  targeting  and  irradiation:  correlation  of  MRI  tissue 
changes  at  cellular  and  molecular  levels  to  optimizing  outcome 


6.  Author(s) 

Dawen  Zhao,  M.D.,  Ph.D. 


DAMD17-03-1-0363 


7.  Performing  Organization  Name  (Include  Name,  City,  State,  Zip  Code  and  Email  for  Principal 
Investigator) 

U.T.  Southwestern  Medical  Center  at  Dallas 
5323  Harry  Hines  Blvd. 

Dallas,  TX  75390-9058 

E-mail:  dawen.zhao@UTSouthwestem.edu 


8.  Performing  Organization  Report  Number 
(Leave  Biank) 


9.  Sponsoring/Monitoring  Agency  Name  and  Address 

U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


10.  Sponsoring/Monitoring  Agency  Report 
Number  (Leave  Biank) 


11.  Supplementary  Notes  (i.e.,  report  contains  color  photos,  report  contains  appendix  in  non-print  form,  etc.) 

12a.  Distribution/Availability  Statement  (check  one) 

X  Approved  for  public  release;  distribution  unlimited 

Distribution  limited  to  U.S.  Government  agencies  only  -  report 
contains  proprietary  information 

12b.  Distribution  Code 

(Leave  Blank) 

13.  Abstract  (Maximum  200  Words)  (abstract  should  contain  no  proprietary  or  confidential  information) 

Vascular  targeting  agents  (VTA)  are  new  types  of  anticancer  drugs  that  act  on  existing  tumor  vasculature, 
causing  vascular  disruption,  which  ultimately  leads  to  extensive  ischemic  tumor  cell  death.  Interesting  findings 
showed  that  VTA  killed  cells  predominantly  in  the  more  hypoxic  tumor  center,  as  a  consequence  of  hemorrhagic 
necrosis  after  vascular  collapse,  whereas  the  better  perfused  peripheral  rim  was  less  affected.  This  apparently 
limits  the  effectiveness  of  such  agents  and  rapid  re-growth  of  tumor  residues  occurs.  However,  these  findings 
suggest  the  possibility  and  promise  of  a  combination  of  VTA  with  treatments  specifically  targeting  the  viable 
tumor  rim.  Radiation  can  certainly  be  expected  to  be  most  effective  against  the  well-perfused  and  oxygenated 
cell  populations  at  the  peripheries  of  the  tumors. 

One  major  goal  of  this  project  is  to  fiilly  understand  and  precisely  assess  the  dynamic  changes  in  blood 
perfusion  and  oxygenation  after  VTA,  so  that  we  may  predict  response  and  optimize  the  therapy.  I  propose  to 
use  in  vivo  MRI  to  measure  and  assess  physiological  changes,  e.g.  tumor  blood  perfusion  and  dynamic  tissue 
oxygenation,  in  the  tumors  before  and  after  treatment.  I  believe  non-invasive  MRI  approaches  may  provide  a 
valuable  prognostic  tool  to  predict  the  response  of  specific  breast  tumors  to  VTA. 


14.  Subject  Terms  (keywords  previously  assigned  to  proposal  abstract  or  terms  which  apply  to  this  award) 

Vascular  targeting.  Magnetic  resonance  imaging  (MRi) ,  p02. 

Perfusion,  irradiation 


15.  Number  of  Pages  (count  all 
pages  including  appendices)  94 


16.  Price  Code 


17.  Security  Classification  of 
Report 

Unclassified 


18.  Security  Classification  of 
this  Page 

Unclassified 


19.  Security  Classification  of 
Abstract 

Unclassified 


20.  Limitation  of  Abstract 


Unlimited 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-18 
298-102 


Table  of  Contents 


Page 


Introduction . 4 

Body . 4-12 

Key  Research  Accomplishments . 12-13 

Reportable  Outcomes . 13-14 

Conclusions . 14 

References . 15-16 

Appendices . 17-94 


Introduction: 


Targeting  tumor  vasculature  promises  effective  cancer  therapy  (1,  2).  It  avoids  issues  of  drug  delivery 
and  is  potentiated  by  massive  downstream  effects  where  one  blood  vessel  may  supply  the  nutrients  for 
thousands  of  tumor  cells.  Thus,  disrupting  the  vascular  supply  should  generate  magnified  tumor  cell  kill. 
Thus,  inhibiting  the  growth  of  new  blood  vessels,  i.e.,  antiangiogenesis,  should  prevent  growth  and 
metastasis  of  the  primary  tumor  (1,  2).  In  addition  to  the  focus  on  the  antiangiogenic  approaches,  vascular 
targeting,  directly  attacking  the  existing  neovasculature,  is  an  alternative  strategy  against  the  tumor  blood 
vessel  network.  Tubulin  binding  agents,  e.g.,  combretastatin  A-4-phosphate  (CA4P)  represent  one  kind  of 
vascular  targeting  agent  (VTA)  (3-5).  Promising  preclinical  studies  have  shown  that  such  agents  selectively 
cause  tumor  vascular  shutdown  and  subsequently  trigger  a  cascade  of  tumor  cell  death  in  experimental 
tumors  (4-6).  However,  survived  tissues  in  a  thin  viable  rim  of  tumor  usually  re-grow  in  spite  of  induction 
of  massive  necrosis.  Thus,  a  combination  of  VTAs  with  additional  conventional  therapeutic  approaches  will 
be  required  (7,  8).  To  better  understand  the  mode  of  action,  and  hence,  optimize  such  combinations,  we  plan 
to  apply  in  vivo  MR  imaging  approaches  to  monitoring  physiological  changes  in  response  to  VTA 
administration.  Dynamic  contrast  enhanced  (DCE)  MRI  based  on  the  transport  properties  of  gadolinium- 
DTPA  (Gd-DTPA)  is  the  most  commonly  used  imaging  approach  to  study  tumor  vascular  perfusion  and 
permeability  (9,  10).  For  combination  with  radiotherapy,  measurement  of  tumor  oxygen  dynamics  will  be 
especially  important  since  hypoxia  affects  radiation  response.  By  applying  FREDOM  (Fluorocarbon 
Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping)  MRI  (11),  dynamic  tumor 
oxygenation  can  be  monitored  following  the  treatment  with  CA4P.  Based  on  pathophysiological  changes 
monitored  by  MRI,  optimum  scheme  of  the  combined  radiation  and  CA4P  will  be  designed  and 
experimental  treatment  will  be  performed  on  the  syngeneic  rat  breast  tumors. 

Successful  completion  of  this  project  will  confirm  the  potential  of  this  new  therapeutic  approach  to 
breast  cancer.  It  will  lay  the  foundation  for  future  clinical  application  for  treating  breast  cancer  patients. 

Body: 

The  Statement  of  Work  in  this  project  had  two  major  tasks: 

Task  1.  To  assess  vascular  and  oxygen  dynamics  in  response  to  VTA,  Months  1-18. 

a.  In  vivo  MRI  assessment  of  vascular  and  oxygen  dynamics  in  response  to  VTA,  Months  1-18 

b.  To  study  morphological  and  biological  changes  of  tumor  vasculature  and  hypoxia  at  cellular  and 
molecular  levels  in  response  to  VTA,  Months  1-18. 

Task  1  was  completed  during  the  Years  1  and  2. 

In  accordance  with  Statement  of  Work  for  Year  1  in  this  project,  I  have  applied  functional  MRI 
approaches  including  DCE  (dynamic  contrast  enhanced)  MRI,  BOLD  (Blood  oxygen  level  dependent)  MR 
and  '^F  FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping) 
NMR  oximetry  to  extensively  studying  tumor  vasculature  and  oxygenation  in  subcutaneous  model  of  13762 
rat  breast  tumors.  For  the  untreated  control  tumors,  I  designed  a  novel  MRI  approach  to  study  potential 
correlation  between  tumor  vasculature  evaluated  by  'h  DCE,  BOLD  MRI  and  tumor  oxygenation  acquired 
by  *^F  FREDOM  oximetry.  A  significant  correlation  was  found  between  changes  in  tissue  p02  and  BOLD 
response  accompanying  oxygen  intervention  in  the  total  of  nine  13762NF  breast  tumors  (R  >  0.9,Fig.  1)(1 1- 
14). 

For  the  treated  tumors,  integrated  functional  MRI  techniques  including  DCE,  BOLD  MRI  and  '^F  MRI 
oximetry  were  applied  to  detect  early  changes  in  tumor  tissue,  vasculature  and  oxygenation  after  a  vascular 
targeting  agent,  Combretastatin  A4  phosphate  (CA4P;  OXiGENE,  Inc,  Waltham,  MA)  injection.  In  terms  of 
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CA4P  dosage,  I  tried  several  doses  ranging  from  30mg/kg  to  150mg/kg.  Our  results  showed  that  even  a 
lower  dose  at  30mg/kg  produeed  significant  effect  on  tumor  vasculature.  So  the  dose  of  30mg/kg  i.p.  was 
used  for  the  whole  study.  DCE  MRI  showed  a  ~70%  decrease,  averaged  over  a  group  of  13762NF  rat  breast 
tumors  (n  =  5),  in  perfusion/permeability  (initial  area  under  signal-intensity  curve  (lAUC))  at  2  h.  The 
lAUC  recovered  fully  after  24  h  in  a  thin  peripheral  region,  but  not  the  tumor  center  (Fig.  2)  (15,  16). 


Mean  BOLD  (%) 

Figure  1  A  significant  linear  correlation  was  found  between  mean  increase  in  '^F  tissue  pOi  (ApOi)  and 
'H  bold  signal  increase  in  the  nine  tumors  with  respect  to  oxygen  intervention  (R  >  0.9;  p  < 
0.001)(13). 
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Figure  2.  Dynamic  contrast  enhanced  (DCE)  MRI  with  respeet  to  CA4P  (30  mg/kg)  treatment.  A) 
Normalized  T1 -weighted  eontrast  enhaneed  images  were  acquired  23  s  after  a  bolus  injection  of  contrast 
agent  pre,  and  2  h  and  24  h  after  treatment.  Significantly  less  signal  enhaneement  was  observed  for  the 
whole  tumor  region  2  h  after  treatment.  A  full  reeovery  was  apparent  in  the  tumor  rim  24  h  post-treatment. 
B)  lAUC  frequency  of  DCE  MRI  on  voxel  by  voxel  basis  obtained  on  the  same  tumor  at  pre,  2  h  and  24  h 
after  CA4P.  Compared  with  pre,  eurves  of  2  h  and  24  h  after  showed  a  signifieant  inerease  in  the  number 
of  voxels  with  no  signal  enhaneement  (typically  lAUC  <  0.05).  However,  the  frequency  of  highly 
enhancing  voxels  (lAUC  >  1.5)  at  24  h  reeovered  to  the  pre-treatment  level  (15). 

FREDOM  MRI  further  demonstrated  that  a  signifieant  decrease  in  the  mean  p02  for  the  group  of  7 
tumors  was  found  at  90  min  after  CA4P  (23  ±  5  vs.  9  ±  3  torr;  p  <  0.05)  and  a  further  deerease  at  2  h  (2  ± 
2  torr).  The  mean  p02  increased  to  15  ±  4  torr  at  24  h,  whieh  was  still  significantly  lower  than  pretreated 
baseline  (p  <  0.05).  More  interestingly,  the  initial  ehanges  in  p02  in  central  and  peripheral  regions  were 
parallel,  but  by  24  h  post  treatment  significant  difference  was  apparent:  p02  in  the  periphery  improved 
signifieantly,  while  the  center  remained  hypoxie  (Figs.  3  and  4)  (15). 
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Figure  3  A)  Distribution  of  hexafluorobenzene  (HFB)  in  a  representative  tumor.  An  overlay  of  ’^F  signal 
density  on  the  anatomie  image  indieates  HFB  in  both  peripheral  and  eentral  regions.  B)  p02  maps  obtained 
from  the  same  tumor  eomprise  two  separated  groups  of  voxels,  whieh  eorrespond  to  the  locations  of  HFB 
on  the  anatomie  image.  Twenty  five  individual  voxels  were  traeeable  from  the  pretreated  baseline  to  2  h 
after  CA4P  with  oxygen  breathing.  Thirty  one  voxels  could  be  followed  24  h  post-treatment.  Signifieant 
deerease  in  p02  was  evident  for  all  the  individual  voxels  after  CA4P,  and  p02  did  not  respond  to  oxygen 
inhalation  after  2  h.  The  24  h  maps  showed  improved  p02  and  significant  response  to  oxygen  breathing  in 
the  peripheral  region  (right),  but  not  in  the  central  region  (left)  (15). 


6 


Figure  4  Mean  pOi  curves  are  shown  for  the  peripheral  (■)  and  central  (o)  voxels  of  the  tumor  shown  in 
Fig.  3.  Significant  decrease  in  p02  was  found  as  early  as  30  min  after  CA4P  (30  mg/kg)  for  both 
peripheral  and  central  tumor.  *  p  <  0.05  from  baseline  air,  ^  p  <0.05  from  24  h  air,  ^  p  <  0.05  from 
periphery  (15). 

DCE  MRI  data  of  tumor  vascular  perfusion/permeability  change  in  response  to  CA4P  is  summarized 
in  Table  1,  and  MRI  p02  data  is  shown  in  Table  2  (15). 


Table  1  Normalized  lAUC  by  dynamic  contrast  enhanced  MRI 


Groups 

Mean  lAUC 

baseline 

2h 

24  h 

CA4P  (n  =  6) 

periphery 

0.88  ±0.06 

0.31  ±0.08*^ 

0.86  ±0.06 

(30  mg/kg) 

center 

0.44  ±0.11^ 

0.13  ±0.08*^ 

0.12  ±0.05*^^ 

Control  (n  =  3) 

periphery 

1.09  ±0.04 

0.93  ±0.13 

0.86  ±0.05 

(saline) 

center 

0.52  ±0.19^ 

0.46  ±  0.05^ 

0.44  ±  0.03^ 

lAUC;  initial  area  under  signal-intensity  curve; 

*  “h  1* 

p  <  0.05  from  baseline;  p  <  0.05  from  periphery;  ’^p  <  0.05  from  control. 


In  addition  to  the  rat  breast  tumor,  I  have  also  performed  the  MRI  and  bioluminescent  imaging 
studies  on  human  breast  MDA-MB-23 1  tumors.  Similar  to  the  observations  in  the  rat  tumor,  DCE  MRI 
revealed  significant  reduction  in  tumor  vascular  perfiision/permeability  in  response  to  CA4P  in  the  human 
tumors  (Fig.  5)(17-19). 
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Table  2  Tumor  oxygen  dynamics  assessed  by  ’^F  MRI  with  respect  to  CA4P  treatment 


Groups 

Case 

no. 

Size 

3 

em 

PO2  (mean  ±  se  torr) 

baseline 

Air 

Time  eourse 

30  min 

60  min 

90  min 

2h 

oxygen 

24  h  base 

oxygen 

1 

2.0 

36  ±  1 

NA 

NA 

NA 

2  ±3* 

* 

5±2 

NA 

NA 

2 

0.9 

12  ±2 

NA 

NA 

NA 

-3  ±2* 

-3  ±2* 

NA 

NA 

CA4P 

3 

1.2 

13  ±0 

35±  11 

11  ±3 

9±4 

1  ±2* 

* 

2±2 

* 

5±0 

30  ±6^ 

(30  mg/kg) 

4 

0.7 

43  ±5 

29  ±3* 

21  ±3* 

20  ±2* 

* 

13  ±  1 

16  ±5* 

20  ±2* 

48  ±7^ 

(n=7) 

5 

0.4 

12  ±2 

5  ±3* 

9±3 

4  ±2* 

2±  1* 

* 

7±3 

NA 

NA 

6 

0.7 

15±  1 

19±3 

* 

9±2 

4  ±2* 

1  ±2* 

* 

0±3 

12  ±  1 

73  ±  12^ 

7 

0.3 

28  ±6 

29  ±2 

16±4 

6  ±2* 

1  ±  1* 

* 

7±2 

23  ±2* 

79  ±6^ 

Mean 

0.9  ±0.2 

23  ±5 

23±  13 

13  ±5 

9  ±3* 

2  ±2* 

* 

5±3 

15±4 

58±  ir 

8 

0.6 

34  ±2 

47  ±5 

40  ±5 

36±  10 

31±8 

88±  10* 

37  ±0 

86  ±5^ 

Saline 

9 

0.2 

47  ±2 

62  ±8 

65  ±6 

64  ±5 

61±4 

114±8* 

67  ±  1 

159  ±7^ 

(n=3) 

10 

0.6 

13  ±  1 

16±6 

10±4 

9±4 

12  ±6 

35±  10* 

14  ±  1 

78±  15^ 

Mean 

0.5  ±0.1 

31±  10 

42±  14 

38±  16 

36±  16 

35±  14 

79  ±  23* 

39±  15 

108  ±26^ 

p  <  0.05  from  baseline;  ^  p  <  0.05  from  24  h  air  (24  h  after  injection  of  CA4P);  NA:  no  measurement. 


T1 -weighted  T2-weighted  T1  eontrast 
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Figure  5.  DCE  MRI  monitoring  of  tumor  response  to  CA4P.  Top  row:  Conventional  Tl-  and  T2- 
weighted,  and  Tl-  weighted  eontrast  enhaneed  MR  images  of  a  nude  mouse  with  MDA-MB-231-luc 
mammary  tumor.  Bottom  row:  Dynamic  contrast  enhance  MRI  was  performed  in  the  mouse  before  (left) 
and  2  h  after  (right)  CA4P  (120  mg/kg)  i.p.  injection.  A  normalized  contrast  enhanced  image  at  30  s  after 
a  bolus  injection  of  Gd-DTPA-BMA  acquired  before  and  2  h  after  treatment,  respectively  is  superimposed 
on  the  Tl -weighted  image.  Significantly  decreased  signal  enhancement,  compared  to  pretreatment,  was 
observed  2  h  after  i.p.  injection  of  CA4P  (17). 
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I  have  also  performed  extensive  histological  studies  on  tumor  blood  vessels,  perfusion,  and  hypoxia  in 
breast  tumors  (Figs.  6  and  7).  Immunohistochemical  studies  of  tumor  hypoxia  and  vasculature  using 
hypoxic  marker  pimonidazole  and  endothelium  marker  CDS  1  showed  that  a  higher  microvascular  density 
(MVD)  and  lower  labeling  index  of  pimonidazole  in  tumor  periphery  than  the  central  regions. 


■  , 

.1 

t 

pimonidazole 

D31 


Figure  6  Comparison  of  pimonidazole  and  CD31  in  a  representative  tumor.  Light  immunostaing  for 
CDS  1  (left)  indicated  extensive  vasculature.  Fluorescent  staining  for  hypoxia  (pimonidazole,  green)  was 
observed  at  some  distance  from  blood  vessels  stained  for  CDS  1  (red). 


Figure  7  Perfusion  marker  Hoechst  staining  pre,  2  h,  and  24  h  post  CA4P  (left).  Vascular  endothelium 
of  the  same  field  was  immunostained  by  anti-CD31  (red,  middle).  A  good  match  (right)  between 
Hoechst  and  anti-CDSl  stained  vascular  endothelium  was  found  in  the  pretreated  tumor.  Two  hours 
after  treatment,  significant  reduction  in  perfused  vessels  was  detected,  followed  by  a  recovery  at  24  h 
point.  Bar  =100  pm  (15). 
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Task  2.  Experimental  tumor  therapy,  Month  19-36 

a.  Learn  and  become  proficient  in  operating  state  of  the  art  irradiation  system  (AccuRay). 

Completed. 


b.  Design  therapeutic  protocol  based  on  the  MR  and  histological  findings:  compare  the  order 

and  timing  of  the  combined  therapy  (CA4-P  lOOmg/kg,  i.p.,  irradiation  30  Gy  single  dose).  Completed. 


Doses  for  CA4P  and  Radiation: 

CA4P  dose:  Based  on  the  MRI  studies,  a  dose  of  30  mg/kg  induced  significant  reduction  in  tumor 
vascular  perfusion/permeability  and  tissue  p02  (20).  Therefore,  we  decided  to  use  this  dose  instead  of  the 
proposed  100  mg/kg  to  investigate  experimental  treatment. 

Radiation  dose:  While  a  30  Gy  single  dose  was  proposed  originally,  we  found  this  dose  was  well  over  the 
TCD50  for  the  proposed  13762NF  breast  tumor.  A  10  Gy  single  dose  was  then  investigated.  Results 
showed  that  this  dose  also  significantly  inhibited  tumor  growth  (Fig.  8).  To  study  potential  effects  by 
adding  CA4P,  we  decided  to  further  lower  the  radiation  dose  to  5  Gy  (21). 
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Figure  8.  Tumor  growth  was  significantly  delayed  by  a  single  dose  of  10  Gy  radiation  (•)  compared  to 
control  tumors  (□). 

Order  and  timing  of  the  combination:  Our  MRI  results  from  Years  1  and  2  have  shown  that  tumor 
blood  perfusion/permeability  decreased  significantly  to  ~30%  of  baseline  pretreatment  level  at  2  h  after 
CA4P  (30  mg/kg,  i.p.)  infusion,  which  recovered  frilly  after  24  h  in  a  thin  peripheral  region,  but  not  the 
tumor  center.  More  importantly,  dynamic  tumor  regional  pOi,  which  is  well  recognized  to  correlate 
closely  with  radiation  outcome,  was  evaluated  by  *^F  MRI.  Tumor  pOi  was  found  to  decline  within  60 
min,  become  significantly  lower  at  90  min,  and  decrease  further  at  2  h  after  CA4P  infusion.  Some  regional 
recovery  was  seen  24  h  later  but  the  p02  was  still  significantly  lower  than  the  pretreatment  level. 
However,  oxygen  breathing  at  this  point  modified  tumor  pOi  significantly,  which  resulted  in  essential 
elimination  of  tumor  hypoxia.  All  the  MRI  findings  have  been  confirmed  by  histological  and 
immunohistological  studies.  These  results  have  recently  been  published  (20). 

c  and  d.  Evaluate  tumor  growth  delay  after  the  vascular  targeting  treatment  and/or  irradiation. 

Completed 
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Based  on  these  observations,  we  proposed  to  administer  CA4P  (30  mg/kg)  on  Day  1,  and  5  Gy  radiation 
on  Day  2,  while  having  animals  breathe  100%  O2  from  20  min  before  to  the  end  of  radiation.  Previous 
studies  by  others  have  demonstrated  that  administration  of  VTA  1  h  post  radiation  produced  better 
improvements  in  tumor  response  than  other  combination  schemes.  Here,  we  plan  to  test  and  compare  our 
combination  approach  with  other  possible  combinations  on  the  13762NF  tumors.  The  preliminary 
treatment  data  has  been  presented  at  national  and  international  conferences  (Fig.  9)  (19,  21). 

Animals  bearing  pedicle  13762NF  tumors  were  grouped  as: 

1 .  Control  without  treatment  (n  =  6); 

2.  CA4P  (30  mg/kg,  i.p.  single  dose)  alone  (n  =  6) 

3.  Radiation  alone  (5  Gy  single  dose,  n  =  6) 

4.  Radiation  (5  Gy)  +  O2  (n  =  5),  The  animals  started  to  breathe  oxygen  (100%  O2  +  1% 

isofluorane)  20  min  before  receiving  a  5  Gy  radiation  delivered  by  Accuray  system. 

5.  Radiation  (5  Gy)  +  CA4P  (1  h  post  Rx,  30  mg/kg,  i.p.  n  =  6) 

6.  Radiation  with  O2  (5  Gy)  +  CA4P  (1  h  post  Rx,  30  mg/kg;  n  =  5);  The  animals  started  to 

breathe  oxygen  (100%  O2  +  1%  isofluorane)  20  min  before  radiation 

7.  CA4P  (30  mg/kg)  +  radiation  (24  h  later,  5  Gy;  n  =  6) 

8.  CA4P  (30  mg/kg)  +  radiation  with  O2  (24  h  later,  5  Gy;  n  =  6).  The  animals  started  to  breathe 
oxygen  (100%  O2  +  1%  isofluorane)  20  min  before  receiving  a  5  Gy  radiation  delivered  by  Accuray 
system. 


IR  +  O2 

CA4P  +  24h  IR 
IR 

(IR  +  O2)  +  Ih  CA4P 
IR  +  Ih  CA4P 


CA4P  +  24h  (IR  +  O2) 


Time  (day) 

Figure  9.  Growth  delay  versus  time  curve  for  the  13762NF  tumors.  Significant  growth  inhibition  was 
observed  in  the  CA4P  +  24h  (IR  +  O2)  treated  group. 
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Tumor  volume  change  (normalized  mean  ±  s.e.)  versus  time  curve  was  plotted,  as  shown  in  Fig.  9. 
Most  tumors  in  radiation  alone  or  the  combination  groups  became  severely  ulcerated,  which  led  to 
termination  of  this  study  on  Day  17.  The  results  showed  that  a  single  dose  (5  Gy;  Group  3)  radiation 
alone  inhibited  tumor  growth  significantly  (p<0.05).  The  relatively  higher  radiosensitivity  of  this  tumor 
line  is  in  a  good  agreement  with  the  lower  hypoxia  fractions  observed  by  MRI.  While  significant  growth 
delay  was  achieved  during  the  first  3  days  after  CA4P  (30  mg/kg,  i.p.;  Group  2),  tumors  started  to  re¬ 
grow  and  caught  up  with  the  Control  group  rapidly  on  Day  7.  This  observation  is  in  common  with  other 
studies  (7,  8).  Oxygen  breathing  didn’t  improve  radiosensitivity  in  tumors  of  Group  4.  This  result  again 
supported  our  MRI  observation  that  this  tumor  line  is  relatively  better  oxygenated,  which  typically  has  a 
hypoxia  fractions  (<  5  torr)  less  than  20%  (20).  Similar  observations  has  been  reported  in  our  previous 
studies  in  the  Dunning  prostate  HI  tumors  that  have  a  hypoxia  range  similar  to  the  breast  I3762NF  tumors 
in  this  study  (22).  The  approach  with  CA4P  1  h  post  radiation  (IR  -I-  Ih  CA4P;  Group  5)  showed  no 
beneficial  effects  over  the  IR  alone  by  Day  10.  However,  tumors  in  this  group  seemed  to  stop  growing 
after  Day  10  while  the  IR  alone  tumors  started  to  grow  rapidly  on  Day  10.  Unfortunately,  longer  term  of 
growth  delay  in  these  tumors  could  not  be  achieved  because  of  tumor  ulceration.  Again,  tumors  in  Group 

6  didn’t  benefit  from  addition  of  oxygen  breathing.  Our  MRI  data  showed  that  significantly  increased 
hypoxic  fractions  induced  by  CA4P  was  still  observed  24  h  later  (20).  Thus,  delivery  of  radiation  at  this 
time  point  will  not  induced  significant  tumor  growth  delay.  This  has  been  proved  in  the  tumors  of  Group 

7  when  compared  to  tumors  in  Group  5  or  6.  However,  oxygen  breathing  24  h  after  CA4P  administration 
was  found  to  essentially  eliminate  tumor  hypoxia  in  the  survived  peripheral  rim  (20),  which  will  improve 
radiosensitivity.  Indeed,  the  tumors  in  Group  8  showed  significantly  slower  growth  rate  from  Day  3  than 
any  other  group  (p  <  0.05). 


Key  Research  Accomplishments; 

•  Assessment  of  tumor  perfusion  and  oxygenation  in  the  untreated  breast  tumors  by  in  vivo  MR 
approaches 

a.  Heterogeneous  vasculature  and  oxygenation  within  a  tumor  or  between  individual  tumors  has 
been  evaluated  in  the  13762NF  rat  breast  tumors  using  integrated  proton  and  fluorine  MR 
techniques. 

b.  Results  of  'H  BOLD  and  DCE  MRI,  providing  information  about  qualitative  vascular 
oxygenation  and  perfusion. 

c.  Significant  correlation  has  been  found  between  changes  in  tissue  pOi  monitored  by  '^F  NMR 
oximetry  and  *H  BOLD  response  accompanying  oxygen  intervention. 

•  Assessment  of  dynamic  perfusion  and  oxygenation  in  the  breast  tumors  in  response  to  a  vascular 
targeting  agent,  Combretastatin  A4  phosphate,  by  in  vivo  MR  approaches 

a.  Significant  reduction  in  tumor  perfusion  was  found  at  early  stage  of  treatment  with  CA4P  (2hr) 
by  *H  DCE  MRI  and  diffusion-weighted  MRI. 

b.  Tumor  p02  started  to  drop  significantly  at  90  min  and  continued  to  decrease  to  a  lowest  value  at 
2hr  after  CA4P. 

c.  Tumor  p02  improved  at  24hr  later,  but  it  was  still  lower  than  the  pretreated  baseline  p02;  pure 
oxygen  inhalation  at  this  time  point  significantly  improved  the  p02,  which  can  be  useful  in  terms 
of  timing  for  a  combination  with  irradiation  in  my  Task  2  study  of  this  project. 
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•  Correlation  of  MR  findings  with  biological  studies 

a.  Immunohistochemical  studies  of  tumor  hypoxia  and  vasculature  using  hypoxic  marker 
pimonidazole  and  endothelium  marker  CDS  1 . 

b.  Consistent  with  MRI  findings,  histological  study  of  tumor  perfusion  using  Hoechst  dye  33342 
showed  a  significant  reduction  in  perfused  vessels  at  2hr  after  CA4P,  which  recovered  24  h  later. 

•  Experimental  therapy 

a.  Based  on  in  vivo  study  of  tumor  physiological  dynamics  evaluated  by  MRI,  we  designed  a 
treatment  scheme  to  administer  CA4P  24  h  before  a  single  dose  radiation  plus  oxygen  inhalation. 

b.  The  results  showed  significantly  slower  tumor  growth  in  this  treatment  group  compared  to  those 
in  other  groups  with  different  combination  order  or  timing. 


Reportable  Outcomes 


Abstracts  (Published  Conference  Proceedings): 
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microcirculation  and  oxygenation  using  a  combination  of  BOLD,  DCE  and  MRI.  iP^  ISMRM,  #222, 
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3.  Zhao,  D.,  Jiang,  L.,  and  Mason,  R.P.  MRI  evaluation  of  tumor  physiological  response  to 
combretastatin  A4  phosphate:  correlation  with  a  combined  radiation  response.  Joint  ISMRM-ESMRMB, 
Berlin,  Germany,  May  2007. 
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4.  Zhao,  D.,  Jiang,  L.,  Adam,  A.,  Hahn,  E.W.,  and  Mason,  R.P.  Tumor  physiological  response  to 
antivascular  agent  combretastatin  A4  phosphate  assessed  by  magnetic  resonance  imaging.  96‘^  AACR, 
#3799,  Anaheim,  CA,  Apr  2005. 
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6.  Zhao,  D.,  Richer,  E.,  Liu,  Li.,  Ya  Ren,  Slavine,  N.,  Shay,  J.  W.,  Antich,  P.  P.,  Mason,  R.  P.  In  vivo 
monitoring  of  antivascular  effects  of  combretastatin  A4  phosphate  in  a  breast  cancer  xenograft  model. 

AACR,  Washington,  DC,  Apr  2006. 
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enhances  therapeutic  response  to  Combretastatin  A4  phosphate.  DoD  Era  of  Hope,  Baltimore,  MD, 
Jun  2008. 
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1.  Zhao,  D.,  Jiang,  L.,  and  Mason,  R.P.  Measuring  changes  in  tumor  oxygenation.  Methods  EnzymoL 
386,378-418,  2004. 
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2.  Zhao,  D.,  Jiang,  L.,  Hahn,  E.W.,  and  Mason,  R.P.  Tumor  physiological  response  to  eombretastatin 
A4  phosphate  assessed  by  MRI.  Int.  J.  Radiat.  Oncol.  Biol.  Phys.  62,  872-80,  2005. 

3.  Zhao,  D.,  Richer,  E.,  Antieh,  P.  P.,  and  Mason,  R.  P.  Antivaseular  effeets  of  eombretastatin  A4 
phosphate  in  breast  eancer  xenograft  assessed  using  dynamic  bioluminescence  imaging  and 
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Manuscript  submitted: 

1.  Zhao,  D.,  Jiang,  L.,  Hahn,  E.  W.,  and  Mason,  R.  P.  Evaluation  of  tumor  oxygenation  by  integrated 
'H  bold  and  MRI.  Magn.  Reson.  Med.  submitted  2008. 

Manuscript  in  preparation: 

1.  Zhao,  D.,  Chang,  K.,  Jiang,  L.,  Mason,  R.P.  In  vivo  monitoring  of  tumor  physiological  response  to 
eombretastatin  A4  phosphate:  eorrelation  with  a  eombined  radiation  response.  Radiat.  Res. 


Conclusion: 

Based  on  the  data  of  in  vivo  tumor  perfusion  and  oxygenation  dynamics  in  response  to  the  vascular 
targeting  agent,  eombretastatin  A-4-phosphate  (CA4P)  evaluated  by  MRI,  we  successfully  designed  a 
seheme  to  combine  the  radiation  treatment  and  CA4P  to  treat  breast  tumors.  This  is  the  major  goal  of  the 
proposed  project.  Moreover,  the  pathophysiologieal  information  will  be  especially  useful  for  designing  a 
eomplicated  scheme,  which  usually  involves  combination  of  fractionated  radiation  and  multiple  dose  of 
systemic  chemotherapy  at  elinical  settings.  1  am  eonfident  that  the  proposed  projeet  will  be  fulfilled  by  the 
next  term. 
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ABSTRACT: 

Tumor  oxygenation  has  been  widely  recognized  as  a  potent  factor  influencing  tumor  response  to  various 
therapies,  especially  radiotherapy.  Recognizing  the  intimate  interplay  of  tumor  oxygenation  and 
vascularization,  we  have  initiated  investigations  to  compare  regional  changes  in  tissue  p02  with 
vascularity  in  syngeneic  breast  NF13762  carcinomas  by  integrated  'H  and  MRI.  Tissue  pOj  dynamics 

were  assessed  using  FREDOM  (Fluorocarbon  Relaxometiy  using  Echo  planar  irnaging  for  Dynamic 
Oxygen  Mapping)  with  hexafluorobenzene,  as  the  reporter  molecule.  Dynamic  Contrast  Enhanced  (DCE) 
H  MRI  based  on  exogenous  Gd-DTPA  and  Blood  Oxygen  Level  Dependent  (BOLD)  'H  MRI  based  on 
endogenous  contrast  deoxyhemoglobin  were  used  to  interrogate  vascular  characteristics.  FREDOM 
revealed  considerable  intra  tumoral  heterogeneity  in  the  distribution  of  pOz  values.  A  mean  pOz  increased 
significantly  fi-om  a  baseline  1 8  ±4  (se)  torr  to  a  maximum  value  78  ±1 6  torr  and  a  mean  HF<1 0  mmHg 
decreased  fi-om  23  %  to  4  %,  in  response  to  oxygen  breathing.  The  increase  in  tissue  pOz  (ApOz) 
correlated  closely  with  BOLD  response  to  oxygen  (r  >0.9,  p<0.001).  However,  DCE  MRI  data  (lAUC) 
showed  no  correlation  with  T  pOz  or  BOLD  data.  The  non-invasive  BOLD  and  DCE  MRI  are  capable  of 
qualitatively  measuring  blood  oxygenation  and  perfosion.  Combination  with  infonnation  on  quantitative 
tissue  oxygenation  by  F  MR  will  provide  better  understanding  of  tumor  physiology  and  response  to 
intervention.  Supported  by  DOD  Breast  Cancer  DAMD  170310363  andNCIROl  EB2762  and P20 
CA86354. 
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Evaluation  of  breast  tumor  microcirculation  and  oxygenation  using  a  combination  of  BOLD,  DCE 

and  19F  MM 


D.  Zhao',  L.  Jiang',  A.  Constantinescu',  E.  W.  Hahn',  R.  P.  Mason' 

‘Radiology,  UT  Southwestern  Medical  Center,  Dallas,  TX,  United  States 

P™B»«“ion  and  metastasis,  as  well  as  response 
(*)•  Recognizing  the  mtimate  interplay  of  tumor  oxygenation  and  blood  flow,  we  have  initiated  investigations  to 
coi^are  regional  changes  in  tissue  pOj  with  vascularity.  We  have  recently  established  a  novel  magnetic  resonance  approach  to  measurins 
helTrl  oxygen  tmsion  FREDOM  (Fluorocarbon  Relaxometiy  using  Echo  planar  imaging  for  Dynamic  Oi^gen  Mapping)  with 
he^fluorobenzene,  as  the  reporter  molecule.  This  technique  allows  us  to  not  only  simultaneously  exaiLe  multiple  speS  Stbm 

Sjt?  ITT  individual  locations  with  respect  to  intervention.  Dynamic  ConLt  Enhanced  (DCE) 

H  MRI  bMed  on  exogenous  Gd-DTPA  and  Blood  Oxygen  Level  Dependent  (BOLD)  based  on  endogen^s  contrast  deoxyhemoS  2 
each  sensitive  to  vascular  chaiactenstics.  Here,  we  apply  these  MRI  approaches  to  evaluate  tunfor  oxygenation  and  vas^rirTId 
investigate  the  potential  correlations  among  data  set  acquired  by  each  technique. 

ma'cJe*?*'  carcinojMs  were  implanted  in  skin  pedicles  on  the  fbieback  of  female  Fisher  rats.  When  the  tumors 

T  r  “  ~  ^  ^  4.7  T  Varian  system.  Each  rat  was  maintained  ZIr 

general  anesthesia  (air  and  l/o  isoflurane).  Hexafluorobenzene  (50  |il)  was  injected  into  both  central  and  peripheral  regions  in  a  single 
central  plane  of  the  tumor  comnal  to  the  rat’s  body.  A  tunable  ('H/'’F)  single-tum  solenoid  coil  (2  cm  in  diamLr  matcKo”he  Sor 
sae)  was  pkced  ^ound  the  turaor-be^g  pedicle.  A  smgle  2mm  slice  parallel  to  the  rat  body  containing  the  strongest  fluorine  signal  was 
chosen  for  the  foDowng  H  DCE  and  BOLD  and  *F  pOj  studies.  The  transverse  relaxation  rate  R2*  was  measured  using  multSadient 
?n  (TR-195nis,  TE=7ms  and  spacing  =6ms)  during  air  or  oxygen  breathing.  After  air  equilibration,  a  series  of 

spm  echo  planar  images  (constant  recovery  time  x  =  500  ms  (=  TR)  and  TE=53  ms)  obtained  for  BOLD  response  measurements  during 

S  EW  of™  wftffi  <»^y8enation  was  assessed  on  thT  same  2nmi  slice  us4  ”F 

Tfrf  ^  ^  tune  resolution.  A  series  of  pOj  maps  was  acquired  over  a  period  of  60  min  with  respiratory  challenge 

^  pOi  was  estimated  using  the  relationship:  pOz  (toir)  =  (Rl-0.0836)/0.001876  (2).  Finally,  Dynamic  Centrist 

Enhaiiced  (DC^  was  performed  on  the  2nim  slice  usmg  Tl-weighted  spin  echo  sequence  (TR=  ISOms,  TE=  18m)  after  a  bolus 
mjection  of  Gd-DTPA-BMA  (0. 1  mmoMcg,  Omniscan)  through  a  tail  vein  catheter.  All  data  analysis  was  based  on  pixel  by  ^Sbasls 

Results:  Each  technique  demonstrates  intra-tumoral  heterogeneity.  As  shown  in  Table  1,  tumor  pOiand  BOLD  SI  increased  and  R2* 
decreased  in  response  to  respiratory  chaUenge,  and  tumors  with  higher  initial  pOi  and  less  hypoxic  fraction  (HF,o)  had  higher  challenged 
pOa  values  (r  >0.8,  p<0  01)  A  significant  correlation  was  ftmnd  between  ApOj  and  BOLD  response  (r  >0.9,  p<0.001.  Fig  1)  while  a  lick 
of  correlation  between  baselme  pOj  and  the  BOLD  SI.  There  was  no  general  correlation  among  pOj ,  and  lAUC. 

Table  1.  Results  of  diverse  MR  approaches _ 


Case 

no. 

Air 

‘^F  pOaCtorr) 
HF,o  (%)  A 

BOLD 

(%) 

Air 

RZ-^  (s-') 
Oxygen 

A 

DCE 

lAUC 

I 

36.3 

0 

136 

14.1 

46.1 

42.6 

-3.5 

0.66 

2 

16.6 

23 

64.5 

9.3 

54.7 

52.6 

-2.1 

1.00 

3 

12.5 

31 

49.6 

7.5 

69.4 

63.8 

-5.6 

0.58 

4 

12.8 

29 

67.7 

12.3 

115.6 

101.6 

-14.0 

0.49 

5 

11.5 

33 

21.5 

2.9 

64.1 

66.7 

2.6 

0.40 

6 

12.3 

21 

19.8 

2.3 

49.9 

47.6 

-2:0 

0.31 

7 

24.3 

25 

47.2 

4.5 

36.2 

33.8 

-2.4 

0.61 

8 

NA 

NA 

NA 

6.1 

90.6 

86.6 

-4.0 

0.87 

Mean 

18.0 

23 

58 

7.4 

65.8 

61.9 

-3.9 

0.62 

7  24.3  25  47.2  4.5  36.2  33.8  -2.4  0.61  0  50  100  150 

8  NA  NA  NA  6.1  90.6  86.6  -4.0  0.87  Ap02(tonr) 

Mean!  18.0  23  58  |  7.4  65.8  61.9  -3,9  062  4^'®’  *  ■  Correlation  between  Ap02  and  BOLD 

*at  t^or  pOj  has  prognostic  value  in  the  clinic  (1).  Using  the  FREDOM  approach  we  are  able  to 
detect  mtra-tumor  differences  m  o^genation.  One  nught  anticipate  that  non-invasive  vascular  dynamics  could  provide  surrogate  markers 
for  tunior  oxypnation.  While  we  have  found  that  BOLD  and  DCE  provide  a  qualitative  indication  of  tumor  vascular  dynamicf  oul^uhs 
Mggest  a  lack  of  correlation  with  tumor  pOz  itself  However,  there  was  a  strong  correlation  between  BOLD  response  and  ApO, 
accompanying  oxygen  challenge.  ^ 
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MRI  evaluation  of  tumor  physiological  response  to  combretastatin  A4  phosphate:  correlation  with  a  combined  radiation 

response 


D.  Zhao\  L.  Jiang\  and  R.  P.  Mason^ 

^Radiology,  University  of  Texas  Southwestern  Medical  Center,  Dallas,  Texas,  United  States 


Introduction:  The  vascular  targeting  agent,  combretastatin  A-4-phosphate  (CA4P)  causes  tumor  vascular  shutdown  inducing  massive  cell  death. 
Although  massive  necrosis  can  be  induced,  tumors  usually  regrow  from  a  thin  viable  rim.  Thus,  a  combination  of  VTAs  with  additional  conventional 
therapeutic  approaches,  e.g.,  radiation,  will  be  required  (1).  For  combination  with  radiotherapy,  measurement  of  tumor  oxygen  dynamics  will  be 
especially  important,  since  reduced  perfusion  can  induce  hypoxia,  potentially  modulating  radiation  response.  Thus,  we  have  assessed  dynamic 
changes  in  tumor  oxygenation  as  compared  with  vascular  perfusion/permeability  after  CA4P  treatment  by  combining  and  MRI  (2).  Based  on 
pathophysiological  changes  monitored  by  MRI,  optimum  scheme  of  the  combined  radiation  and  CA4P  treatment  was  designed  and  experimental 
treatment  was  initiated  on  rat  breast  tumors. 


Materials  and  Methods:  Rat  mammary  carcinoma  13762NF  was  implanted  syngeneically  in  a  skin  pedicle  surgically  created  on  the  foreback  of 
Fisher  344  adult  female  rats  and  allowed  to  grow  to  ~  I  cm  diameter. 

MRI  study:  MRI  studies  were  performed  using  a  4.7  T  Varian  Inova  imaging  system.  Each  rat  was  maintained  under  general  anesthesia  (air  and  1% 
isoflurane).  A  tunable  (^H/^^F)  volume  RF  coil  was  placed  around  the  tumor-bearing  pedicle.  MRI  R2*  maps  were  obtained  before  and  2  h  after 
CA4P  (30  mg/kg,  i.p.,  OXiGENE,  Inc.)  by  gradient  echo  sequence  (GEMS)  with  8  echoes  (TR  =195  ms,  TE  =  7  ms  and  spacing  =  6  ms).  Dynamic 
contrast  enhanced  (DCE)  MRI  using  a  Ti-weighted  spin  echo  sequence  (TR  =  70  ms,  TE  =  12  ms)  based  on  i.v.  bolus  injection  of  Gd-DTPA-BMA 
through  a  tail  vein  catheter  was  also  acquired  before,  2  h  and  24  h  after  CA4P.  Eor  ^^E  NMR  oximetry,  hexafluorobenzene  (50  ql)  was  injected 
directly  into  the  tumor  along  two  or  three  tracks  in  a  single  central  plane  of  the  tumor  using  a  fine  sharp  needle  (32G),  as  described  in  detail 
previously  (25).  FREDOM  (Eluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping)  MRI  was  performed  to  acquire  a 
series  of  p02  maps  under  air  or  oxygen  breathing  before  and  at  different  time  points  after  CA4P.  Data  analysis  was  carried  out  on  a  voxel-by- voxel 
basis  with  IDL  based  house  made  software. 

Experimental  treatment:  Animals  bearing  13762NE  tumors  (n  =  30)  were  grouped  as:  1)  control  without  treatment;  2)  CA4P  (30  mg/kg,  i.p.)  alone; 
3)  Radiation  (IR)  alone  (5  Gy  single  dose);  4)  IR  (5  Gy)  +  CA4P  (1  h  post  IR,  30  mg/kg,  i.p.);  5)  CA4P  (30  mg/kg)  +  IR  plus  O2  (24  h  post  CA4P,  5 
Gy).  The  animals  started  to  breathe  oxygen  (100%  02+1%  isoflurane)  20  min  before  receiving  a  5  Gy  IR  delivered  by  Accuray  system. 
Immunohistochemistry:  Immuno staining  for  Hoechst  33342  (perfusion  marker)  and  CD31  (vascular  endothelium)  was  performed  to  correlate  with 
imaging  findings. 


Results:  MRI  showed  that  tumor  blood  perfusion/permeability  by  DCE  MRI  decreased  significantly  to  -30%  of  baseline  pretreatment  level  at  2  h 

after  CA4P  (30  mg/kg,  i.p.)  infusion,  which  recovered  fully  after  24  h  in  a  thin  peripheral  region,  but  not  the  tumor  center.  Analysis  of  R2*  maps 
revealed  significantly  increased  values  after  2  h,  compared  to  pretreatment  (88.6  v.s.  85.1  s  \  p  <  0.05).  Tumor  p02by  ^^E  MRI  was  found  to  decline 
within  60  min,  become  significantly  lower  at  90  min,  and  decrease  further  at  2  h  after  CA4P  infusion.  At  this  time  there  was  no  response  to  breathing 
O2.  Some  regional  recovery  was  seen  24  h  later,  but  the  p02  was  still  significantly  lower  than  the  pretreatment  level.  However,  oxygen  breathing  at 
24  h  point  modulated  tumor  p02  significantly,  which  resulted  in  essential  elimination  of  tumor  hypoxia  (Eig.  1).  Correlating  well  with  MRI 
observations,  the  tumors  of  Group  5  with  radiation  plus  oxygen  24  h  post  CA4P  showed  a  significantly  prolonged  growth  delay  (  p  <  0.05,  Eig.  2). 
Histological  data  using  the  perfusion  marker,  Hoechst  33342,  confirmed  a  significant  decrease  in  perfused  vessels  2  h  after  CA4P,  while  there  was 
recovery  evident  at  24  h. 


Baseline  ^  CA4P  (30  mg/kg)  CA4P24hr 


Fl§’  hdh  lAIK  ujomfouu^  o^fnmoi.  oxXSgu  qXugiiJicg  lu  i.G8bou8G  fo  cVIb 


IR 

IR  (  LAW 
CA4P+IR+0, 


^  0  5  7  10  14 

XuDi  (daY) 

Fig.  2.  Growth  delay  versus  tinie  curve  for  the  13762NF  tumors. 


Discussion:  There  is  a  distinct  similarity  between  the  results  of  the  p02  measurements  and  the  more  traditional  DCE,  but  the  quantitative  p02  values 
provide  the  potential  for  exploiting  synergy  with  other  oxygen  dependent  therapies.  The  observations  further  demonstrate  the  value  of  FREDOM  in 
assessing  dynamic  changes  in  regional  tumor  p02  in  vivo  in  response  to  intervention.  We  believe  that  dynamic  measurements  are  particularly 
valuable  for  understanding  the  mode  of  action  of  therapeutic  response  to  VTAs.  Most  significantly,  these  measurements  lay  a  foundation  to  optimize 
the  timing  of  combination  therapy  involving  fractionated  radiotherapy  and  multiple  doses  of  VTAs. 

References:  1.  Stevenson,  J.  P.,  Rosen,  M.,  Sun,  W.,  Gallagher,  M.,  et  al.  J  Clin  Oncol,  21:  4428-4438,  2003. 
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IN  VIVO  MULTIMODAL  IMAGING  OF  ANTIVASCULAR  EFFECT  OF  COMBRETASTATIN  A4  PHOSPHATE 
ENHANCES  TUMOR  RESPONSE  TO  RADIATION 

Dawen  Zhao,  M.D.,  PhD.,  Karen  Chang,  PhD.,  Edmond  Richer,  PhD.,  Peter  Antich,  PhD.,  Ralph  P.  Mason, 
PhD. 

Department  of  Radiology,  University  of  Texas  Southwestern  Medical  Center,  Dallas,  TX. 

The  vascular  disrupting  agent  (VDA),  Combretastatin  A-4-phosphate  (CA4P)  causes  tumor  vascular 
shutdown  inducing  massive  cell  death.  We  have  applied  magnetic  resonance  imaging  (MRI)  and 
bioluminescent  imaging  (BLI)  to  evaluate  pathophysiological  response  to  CA4P  in  breast  tumors.  Although 
massive  necrosis  can  be  induced  by  CA4P,  tumors  usually  regrow  from  a  thin  viable  rim.  Thus,  a  combination 
of  VDAs  with  additional  conventional  therapeutic  approaches,  e.g.,  radiation,  will  be  required.  Based  on 
imaging  data  of  temporal  changes  in  tumor  perfusion  and  especially  hypoxiation,  induced  by  CA4P,  we  have 
enhanced  the  combined  treatment  with  radiation  in  a  breast  tumor  model. 

For  BLI  study  of  acute  vascular  perfusion  change,  MDA-MB  231  human  breast  tumors  were  infected  to 
stably  express  firefly  luciferase  and  highly  expressing  clones  isolated  and  implanted  subcutaneously  in  the 
thigh  of  a  nude  mouse.  CA4P  (120  mg/kg;  OXiGENE,  Inc.  Waltham,  MA)  was  injected  i.p.  immediately  after 
baseline  BLI  and  then  2  h  and  24  h  later  the  BLI  time  course  was  repeated.  In  the  CA4P  treated  group,  the 
detected  light  emission  decreased  between  50%  and  90%  and  time  to  maximum  was  significantly  delayed. 
Twenty-four  hours  later,  there  was  some  recovery  of  light  emission.  Comparable  vascular  changes  were 
evidenced  by  dynamic  contrast  enhance  (DCE)  ^H  MRI,  and  further  validated  by  histological  study.  Tumor 
hypoxiation  induced  by  CA4P  in  syngeneic  13762NF  rat  breast  tumors  was  evaluated  by  using  the  ^®F  NMR 
oximetry  approach,  FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen 
Mapping).  Tumor  p02was  found  to  decline  within  60  min,  become  significantly  lower  at  90  min,  and  decrease 
further  at  2  h  after  CA4P  infusion.  At  this  time  there  was  no  response  to  breathing  O2.  Some  regional  recovery 
was  seen  24  h  later,  but  the  PO2  was  still  significantly  lower  than  the  pretreatment  level.  Oxygen  breathing  at 
the  24  h  point  modulated  tumor  PO2  significantly,  which  essentially  resulted  in  elimination  of  tumor  hypoxia. 
Thus,  a  combination  treatment  with  radiation  plus  oxygen  24  h  post  CA4P  was  tested  and  compared  with 
alternatives.  The  results  showed  that  the  optimized  combination  scheme  led  to  a  significantly  slower  growth 
rate  than  any  other  group  (p  <  0.05). 

Both  BLI  and  MRI  revealed  tumor  vascular  shutdown  after  CA4P  treatment.  We  believe  that  dynamic 
measurements  are  particularly  valuable  for  understanding  the  mode  of  action  of  therapeutic  response  to  VTAs. 
Most  significantly,  these  measurements  lay  a  foundation  to  optimize  the  timing  of  combination  therapy 
involving  fractionated  radiotherapy  and  multiple  doses  of  VTAs. 
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3799  Tumor  physiological  response  to  antivascular  agent  combretastatin  A4  phosphate 
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Medical  Ctr.,  Dallas,  TX. 

Tumor  oxygenation  has  been  widely  recognized  as  a  potent  factor,  which  influences  tumor  response 
to  chemo-  and  radio-therapy.  The  vascular  targeting  agent,  combretastatin  A-4-phosphate  (CA4P) 
causes  tumor  vascular  shutdown  and  a  combination  of  CA4P  with  radiotherapy  has  shown  enhanced 
tumor  response.  We  have  recently  established  a  novel  magnetic  resonance  approach  to  measuring 
regional  tumor  oxygen  tension  FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for 
Dynamic  Oxygen  Mapping)  with  hexafluorobenzene,  as  a  reporter  molecule.  Here,  the  effect  of 
CA4P  on  real-time  tumor  oxygenation  and  perfusion/permeability  was  monitored  in  vivo  in  NF 13762 

rat  breast  carcinomas  by  ^^F  FREDOM  and  ^H  dynamic  contrast  enhanced  (DCE)  magnetic 
resonance  imaging  (MRI).  Syngeneic  NF  13762  rat  breast  carcinomas  were  implanted  in  a  skin 
pedicle  on  the  foreback  of  Fisher  344  adult  female  rats.  When  tumors  reached  ~  1  cm  diameter  (~  0.6 

cm^)  MRI  was  performed  using  a  4.7  T  horizontal  bore  magnet  with  a  Varian  Unity  Inova  system.  A 
series  of  p02  maps  acquired  at  10,  30,  60,  90,  and  120  minutes  after  i.p.  CA4P  (30  mg/kg; 

OXiGENE,  Inc.  Waltham,  MA)  showed  a  significant  decrease  in  p02  of  some  tumors  as  early  as  30 
min  post  injection.  For  all  tumors  there  was  a  significant  drop  in  tumor  p02  within  90  min  after 
treatment  (mean  baseline  p02  =  23  torr  to  9  torr,  p  <  0.05)  and  a  further  decrease  was  observed  at  2  h 
(mean  =  2  torr;  p  <  0.01).  Intriguingly,  the  initial  changes  in  p02  in  central  and  peripheral  regions 
were  parallel,  but  by  24  h  post  treatment  significant  difference  was  apparent:  the  p02  in  periphery 

improved  significantly,  while  center  remained  hypoxic.  These  data  are  consistent  with  DCE  MRI, 
which  revealed  a  -70%  decrease  in  perfusion/permeability  (initial  area  under  signal-intensity  curve 
(lAUC))  at  2  h  after  30  mg/kg  CA4P  (p  <  0.001).  lAUC  fully  recovered  in  a  thin  peripheral  region, 
but  not  the  tumor  center  24  h  post-treatment.  Vascular  perfusion  marker  Hoechst  33342  staining 
confirmed  that  the  number  of  perfused  vessels  decreased  significantly  at  2  h,  and  then  recovered  at 
24  h.  We  conclude  that  magnetic  resonance  imaging  can  monitor  in  vivo  dynamic  changes  in  tumor 
vasculature  and  oxygenation  following  VTA  administration.  These  results  provide  the  first  insight 
into  regional  tumor  oxygen  dynamics  in  response  to  CA4P  in  a  syngeneic  rat  tumor.  While  dynamic 
contrast  MRI  indicates  vascular  shut  down,  the  critical  p02  measurements  are  potentially  more 

important  for  optimizing  therapeutic  combination  of  VTAs  with  conventional  therapy,  in  particular, 
irradiation.  Supported  by  POD  Breast  Cancer  IDEA  Award  (DAMP  170310363).  in  conjunction 
with  NCI  ROl  CA79515/EB002762  and  the  Cancer  Imaging  Program,  a  P20  Pre-ICMIC  CA86354. 
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IN  VIVO  MRI  MONITORING  OF  BREAST  TUMOR  RESPONSE  TO  THE 
VASCULAR  TARGETING  AGENT  COMBRETASTATIN  A4  PHOSPHATE 
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Hypoxic  cells  in  tumors  has  long  been  recognized  as  a  significant  factor  influencing 
tumor  response  to  chemo-  and  radio-therapy.  The  vascular  targeting  agent,  combretastatin 
A-4-phosphate  (CA4P)  causes  tumor  vascular  shutdown  and  a  combination  of  CA4P  with 
radiotherapy  has  shown  enhanced  tumor  response.  For  the  potential  combination  with 
radiotherapy,  measurement  of  tumor  oxygen  dynamics  will  be  especially  important.  Here, 
we  applied  a  novel  magnetic  resonance  imaging  (MRI)  approach  to  measuring  regional 
tumor  oxygen  tension:  FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging 
for  Dynamic  Oxygen  Mapping)  with  hexafluorobenzene,  as  a  reporter  molecule. 
FREDOM  allowed  us  to  quantitatively  determine  the  in  situ,  in  vivo  dynamic  action  of 
CA4P  on  real-time  tumor  oxygenation. 

Syngeneic  NF 13762  rat  breast  carcinomas  were  implanted  in  a  skin  pedicle  on  the 
foreback  of  Fisher  344  adult  female  rats.  When  tumors  reached  ~  1  cm  diameter  (~  0.6 
cm3)  MRI  was  performed  using  a  4.7  T  horizontal  bore  magnet  with  a  Varian  Unity 
Inova  system.  A  series  of  p02  maps  acquired  at  10,  30,  60,  90,  and  120  minutes  after  i.p. 
CA4P  (30  mg/kg;  OXiGENE,  Inc.  Waltham,  MA)  showed  a  significant  decrease  in  p02 
of  some  tumors  as  early  as  30  min  post  injection.  For  all  tumors  there  was  a  significant 
drop  in  tumor  p02  within  90  min  after  treatment  (mean  baseline  p02  =  23  torr  to  9  torr,  p 
<  0.05)  and  a  further  decrease  was  observed  at  2  h  (mean  =  2  torr;  p  <  0.01).  Intriguingly, 
the  initial  changes  in  p02  in  central  and  peripheral  regions  were  parallel,  but  by  24  h  post 
treatment  significant  difference  was  apparent:  the  p02  in  periphery  improved 
significantly,  while  center  remained  hypoxic.  In  contrast  to  tumor  behavior,  there  was  no 
significant  change  in  p02  of  thigh  muscles  of  rat  at  any  time  point  after  CA4P.  These 
data  are  consistent  with  dynamic  contrast  enhanced  (DCE)  MRI,  which  revealed  a  ~70% 
decrease  in  perfusion/permeability  (initial  area  under  signal-intensity  curve  (lAUC))  at  2 
h  after  30  mg/kg  CA4P  (p  <  0.001).  lAUC  fully  recovered  in  a  thin  peripheral  region,  but 
not  the  tumor  center  24  h  post-treatment.  Vascular  perfusion  marker  Hoechst  33342 
staining  confirmed  that  the  number  of  perfused  vessels  decreased  significantly  at  2  h,  and 
then  recovered  at  24  h. 

These  results  provide  the  first  insight  into  regional  tumor  oxygen  dynamics  in  response  to 
CA4P  in  a  syngeneic  rat  tumor.  While  dynamic  contrast  enhanced  MRI  will  reveal 
vascular  changes,  dynamic  p02  measurements  are  potentially  more  important  for 
optimizing  combinations  of  VTAs  with  conventional  anti-cancer  therapies  against  breast 
cancers. 


The  U.S.  Army  Medical  Research  and  Materiel  Command  under  DAMDl  7-03-1 -0363 
supported  this  work. _ 
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The  vascular  targeting  agent,  combretastatin  A-4-phosphate  (CA4P)  causes  tumor  vascular  shutdown  inducing 
massive  cell  death.  We  have  recently  shown  acute  hypoxiation  within  90  mins  following  CA4P  administration  to 
rats  bearing  syngeneic  breast  13762NF  tumors  using  MRI.  We  have  now  applied  MRI  and  bioluminescent  imaging 
(BLI)  to  probe  the  acute  effects  of  CA4P  on  human  breast  MDA-MB-23 1  tumors. 

231  cells  were  infected  with  a  lend  virus  expressing  a  lucif erase  reporter  and  highly  expressing  clones  isolated.  10^ 
cells  were  implanted  in  the  flank  of  nude  mice  and  allowed  to  grow  to  ~  6  mm  diameter.  For  BLI  studies,  mice 
were  anesthetized  (isoflurane/02)  and  a  solution  of  D-luciferin  (450  mg/kg)  was  administered  s.c.  in  the  neck 

region  and  light  images  acquired  immediately  using  one  camera  of  our  Light  Emission  Tomography  System 
(LETS).  Serial  images  (30  s  each)  were  acquired  over  a  period  of  20  -  30  mins  and  the  light  intensity-time  curves 
evaluated.  Saline  or  CA4P  in  saline  (120  mg/kg;  OXiGENE,  Inc.  Waltham,  MA)  were  injected  i.p.  immediately 
after  baseline  BLI  and  then  2  h  and  24  h  later  the  BLI  time  course  was  repeated.  We  also  undertook  3D  imaging  by 
acquiring  multiple  images  simultaneously  using  3  cameras  arranged  in  a  circular  gantry  around  the  mouse.  MRI 
studies  were  performed  using  a  4.7  T  Varian  Inova  imaging  system.  The  transverse  relaxation  rate  R2*  was 

measured  using  multigradient  echo  sequence  before  and  2  h  after  i.p.  CA4P  (120  mg/kg).  Dynamic  contrast 
enhanced  (DCE)  MRI  based  on  i.v.  bolus  injection  of  Gd-DTPA-BMA  through  a  tail  vein  catheter  was  also 
acquired  before  and  2  h  after  CA4P,  respectively. 

Control  tumors  showed  intense  light  emission  peaking  within  8  mins  and  generally  decreasing  to  about  50-70% 
after  20  mins.  By  contrast,  CA4P  led  to  a  significantly  lower  light  emission  (peak  ~  2  to  10  times  lower)  and 
delayed  peak  emission  when  animals  were  imaged  after  2  h.  24  h  later  signal  remained  considerably  decreased. 
Traditional  BLI  provides  a  planar  image  only,  but  tumors  are  3D  and  known  to  exhibit  heterogeneity,  particularly, 
after  vascular  targeting  agents.  LETS  successfully  provides  a  3D  representation  of  the  tumors.  In  good  agreement 
with  the  BLI  data,  DCE  MRI  revealed  a  -70%  decrease  in  perfusion/permeability  of  tumors  2  h  after  CA4P  (p  < 
0.001),  while  little  change  was  observed  in  perfusion  of  the  femoral  artery.  R2*  measurement  showed  a  significant 

increase  in  R2*  values  2  h  after  CA4P  treatment  (mean  of  2  h  =  131  s'^  vs  pre  =  1 13  s’^  p  <  0.01),  which  may 

suggest  an  elevated  deoxyhemoglobin  from  hemorrhagic  thrombosis. 

Both  BLI  and  MRI  enabled  accurate  imaging  of  tumor  vascular  shutdown  after  CA4P  treatment.  However,  BLI  is 
much  cheaper  and  offers  a  high  throughput  method  for  evaluating  novel  drugs  and  drug  combinations  and 
scheduling. 
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Dynamic  in  vivo  imaging  of  breast  tumors  enhances  therapeutic  response  to  Combretastatin  A4  phosphate 
Dawen  Zhao,  Karen  Chang,  Edmond  Richer,  Nikolai  Slavine,  Peter  Antich,  Ralph  P.  Mason 
University  of  Texas  Southwestern  Medical  Center,  Dallas,  TX,  United  States 

The  vascular  disrupting  agent  (VDA),  Combretastatin  A-4-phosphate  (CA4P)  causes  tumor  vascular 
shutdown  inducing  massive  cell  death.  We  have  applied  magnetic  resonance  imaging  (MRI)  and 
bioluminescent  imaging  (BLI)  to  evaluate  pathophysiological  response  to  CA4P  in  breast  tumors.  Although 
massive  necrosis  can  be  induced  by  CA4P,  tumors  usually  regrow  from  a  thin  viable  rim.  Thus,  a  combination 
of  VDAs  with  additional  conventional  therapeutic  approaches,  e.g.,  radiation,  will  be  required.  Based  on 
imaging  data  of  temporal  changes  in  tumor  perfusion  and  especially  hypoxiation,  induced  by  CA4P,  we  have 
enhanced  the  combined  treatment  with  radiation  in  a  breast  tumor  model. 

For  BLI  study  of  acute  vascular  perfusion  change,  MDA-MB  23 1  human  breast  tumors  were  infected  to 
stably  express  firefly  luciferase  and  highly  expressing  clones  isolated  and  implanted  subcutaneously  in  the 
thigh  of  a  nude  mouse.  CA4P  (120  mg/kg;  OXiGENE,  Inc.  Waltham,  MA)  was  injected  i.p.  immediately  after 
baseline  BLI  and  then  2  h  and  24  h  later  the  BLI  time  course  was  repeated.  In  the  CA4P  treated  group,  the 
detected  light  emission  decreased  between  50%  and  90%  and  time  to  maximum  was  significantly  delayed. 
Twenty-four  hours  later,  there  was  some  recovery  of  light  emission.  Comparable  vascular  changes  were 
evidenced  by  dynamic  contrast  enhance  (DCE)  'H  MRI,  and  further  validated  by  histological  study.  Tumor 
hypoxiation  induced  by  CA4P  in  syngeneic  13762NF  rat  breast  tumors  was  evaluated  by  using  the  '^F  NMR 
oximetry  approach,  FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen 
Mapping).  Tumor  p02  was  found  to  decline  within  60  min,  become  significantly  lower  at  90  min,  and 
decrease  further  at  2  h  after  CA4P  infusion.  At  this  time  there  was  no  response  to  breathing  O2.  Some  regional 
recovery  was  seen  24  h  later,  but  the  p02  was  still  significantly  lower  than  the  pretreatment  level.  Oxygen 
breathing  at  the  24  h  point  modulated  tumor  p02  significantly,  which  essentially  resulted  in  elimination  of 
tumor  hypoxia.  Thus,  a  combination  treatment  with  radiation  plus  oxygen  24  h  post  CA4P  was  tested  and 
compared  with  alternatives.  As  shown  in  the  figure,  the  optimized  combination  scheme  led  to  a  significantly 
slower  growth  rate  than  any  other  group  (p  <  0.05). 

Both  BLI  and  MRI  revealed  tumor  vascular  shutdown  after  CA4P  treatment.  The  quantitative  p02  values 
measured  by  FREDOM  provide  the  potential  for  exploiting  synergy  with  other  oxygen  dependent  therapies. 
Figure.  Growth  delay  versus  time  curve  for  cohorts  of  13762NF  breast  tumors. 
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ABSTRACT 

We  previously  suggested  that  the  two  time  constants  quantified  from  the  increase  of  tumor  oxyhemoglobin 
concentration,  A[Hb02],  during  hyperoxic  gas  intervention  are  associated  with  two  blood  fiow/perfusion  rates  in  well 
perfused  and  poorly  perfused  regions  of  tumors.  In  this  study,  our  hypothesis  is  that  when  cancer  therapy  is  applied  to  a 
tumor,  changes  in  blood  perfusion  will  occur  and  be  detected  by  the  NIRS.  For  experiments,  systemic  chemotherapy, 
cyclophosphamide  (CTX),  was  applied  to  two  groups  of  rats  bearing  syngeneic  13762NF  mammary  adenocarcinomas: 
one  group  received  a  single  high  dose  i.  p.  (200  mg/kg  CTX)  and  the  other  group  continuous  low  doses  (20  mg/kg  CTX 
i.  p.  for  10  days).  Time  courses  of  changes  in  tumor  A[Hb02]  were  measured  at  four  different  locations  on  the  breast 
tumors  non-invasively  with  an  inhaled  gas  sequence  of  air-oxygen-air  before  and  after  CTX  administration.  Both  rat 
body  weight  and  tumor  volume  decreased  after  administration  of  high  dose  CTX,  but  continuous  low  doses  showed 
decrease  of  tumor  volume  only.  Baselines  (without  any  therapy)  intra-  and  inter-tumor  heterogeneity  of  vascular 
oxygenation  during  oxygen  inhalation  were  similar  to  our  previous  observations.  After  CTX  treatment,  significant 
changes  in  vascular  hemodynamic  response  to  oxygen  inhalation  were  observed  from  both  groups.  By  fitting  the 
increase  of  A[Hb02]  during  oxygen  inhalation,  we  have  obtained  changes  of  vascular  structure  ratio  and  also  of  perfusion 
rate  ratio  before  and  after  chemotherapy.  The  preliminary  results  suggest  that  cyclophosphamide  has  greatest  effect  on 
the  well  perfused  tumor  vasculature.  Overall,  our  study  supports  our  earlier  hypothesis,  proving  that  the  effects  of 
chemotherapy  in  tumor  may  be  monitored  non-invasively  by  using  NIRS  to  detect  changes  of  hemodynamics  induced 
with  respiratory  challenges. 

Keywords:  Breast  Cancer,  Cyclophosphamide,  Hemodynamics,  NIR  Spectroscopy,  Tumor  vascular  oxygenation 


1.  INTRODUCTION 

In  addition  to  surgical  resection,  many  other  types  of  cancer  therapy  are  available  for  patients  including  radiotherapy, 
photodynamic  therapy  and  chemotherapy.  Chemotherapy  plays  an  important  role  to  treat  cancers  even  though  it  has 
some  side  effects.  Currently,  the  effect  of  chemotherapy  is  monitored  by  MRI  or  CT  that  can  measure  the  tumor  volume 
changes  during  cancer  treatment.  However,  it  can  take  up  to  3  weeks  to  detect  such  changes,  and  this  is  considerably  late 
for  clinicians  to  decide  whether  initial  therapeutic  strategy  should  be  continued  or  modified.  This  delay  in  detection  of 
chemotherapy  effect  can  reduce  the  quality  of  a  patient’s  life  and  ineffective  therapy  is  costly.  Therefore,  many 
researchers  are  trying  to  develop  tools  that  can  detect  the  early  response  to  cancer  treatment.  For  example,  Li  et  al  have 
used  nuclear  magnetic  resonance  spectroscopy  (NMRS)  to  measure  the  effectiveness  of  cyclophosphamide  (CTX) 
treatment  in  radiation-induced  fibrosarcoma  (RIF).^  They  found  that  the  ratio  of  inorganic  phosphate  to  other  phosphate 
metabolites  in  CTX  treated  group  was  significantly  decreased  during  the  tumor  growth  delay  period  compared  to  age- 
matched  controls.  Poptani  et  al  studied  the  effects  of  CTX  treatment  in  RIF-1  tumors  in  terms  of  tumor  oxygenation 
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and  glycolytic  rate  changes  by  utilizing  MRS,  Eppendorf  electrode,  and  Redox  scanning.^  They  observed  that  CTX 
treatment  caused  reduction  in  glycolytic  rate,  a  significant  decrease  in  tumor  tissue  p02,  and  also  an  increase  of  NADH 
levels  24  hours  after  the  treatment  while  tumor  volume  did  not  show  any  significant  difference  between  the  CTX-treated 
and  control  groups.  Zhao  et  al.  have  reported  significant  changes  in  rat  breast  tumor  perfusion  following  either  single 
dose  CTX  or  continuous  low  dose  “metronomic”  therapy.^ 

In  the  last  decade,  near  infrared  spectroscopy  (MRS)  has  been  developed  to  examine  tissue  oxygenation  and  has  been 
widely  applied  to  investigate  hemoglobin  oxygenations  of  muscles,"^’^’^  the  brain,^’^’^  and  animal  tumors. Since 
tumors  have  higher  vascular  density  and  also  higher  metabolism  than  normal  tissues,  total  hemoglobin  concentration 
([Hbtotai]),  oxyhemoglobin  concentration  ([Hb02]),  and  reduced  scattering  coefficient  (pis’)  were  used  as  markers  to 
identify  tumors  from  the  human  breast  by  using  MRS.^"^’^^’^^  In  addition,  an  MR  spectrometer  is  a  low-cost,  portable, 
and  real-time  display  instrument.  Therefore,  MRS  has  a  good  potential  to  be  used  as  a  monitoring  tool  for  tumor 
treatment  planning  and  tumor  prognosis. 

We  have  previously  studied  breast  tumor  oxygenation  under  gas  intervention  using  MRS  and  found  that  oxyhemoglobin 
concentration  changes  (A[Hb02])  during  gas  intervention  can  be  fitted  by  a  two-exponential  equation  containing  two 
time  constants.  Based  on  the  model,  we  formed  a  hypothesis  that  changes  in  oxygenated  hemoglobin  concentration 
result  from  well  perfused  and  poorly  perfused  regions  of  an  animal  tumor  to  explain  why  there  are  two  different  time 
constants  in  the  A[Hb02]  data.  The  model  further  allows  us  to  associate  the  signal  amplitudes  and  time  constants  to  the 
ratio  of  vascular  density  and  the  ratio  of  the  perfusion  rates  in  the  two  different  regions,  respectively. 

In  this  study,  we  applied  the  MRS  system  to  monitor  the  tumor  oxygenation  changes  during  oxygen  intervention  before 
and  after  CTX  administration.  The  purpose  for  this  study  is  to  explore  the  MRS  as  a  possible  tool  for  monitoring  tumor 
responses  to  chemotherapy.  This  work  is  based  on  the  following  hypothesis:  when  tumor  is  treated  with  chemotherapy, 
changes  in  blood  perfusion  and  vascular  density  in  the  tumor  will  occur  and  will  be  seen  as  changes  of  the  two  fitted 
parameters  from  the  MRS  measurements.  In  addition,  by  developing  a  non-in vasive  tool  for  monitoring  cancer  therapy, 
we  are  not  only  monitoring  the  reduction  of  tumor  size,  but  also  detecting  the  changes  of  tumor  physiological  conditions, 
which  are  essential  for  tumor  treatment  planning  and  tumor  prognosis. 


2.  MATERIALS  AND  METHODS 


2.1  Tumor  Model  and  Experimental  Procedure 

Rats  were  divided  into  three  groups  for  this  study.  Two  groups  were  treated  with  CTX  at  different  doses,  and  the  other 
group  was  administered  saline  instead  of  CTX  as  a  control  group.  Cyclophosphamide  was  chosen  as  a  chemotherapeutic 
agent  for  this  study  since  our  tumor  line  is  highly  responsive  to  alkylating  agents  and  platinum  chemotherapeutic 
agents. CTX  is  an  antineoplastic  alkylating  agent,  and  it  has  been  used  to  treat  lymphomas,  cancers  of  the  ovary,  breast 
and  bladder,  and  chronic  lymphocytic  leukaemia.  The  tumor  line  was  rat  mammary  adenocarcinomas  13762M^ 
(cells  originally  provided  by  the  Division  of  Cancer  Therapeutics,  NCI),  and  the  tumors  were  implanted  in  the  hind  limb 
of  adult  female  Fisher  344  rats  (-200  g). 

The  rats  were  anesthetized  with  0.2  ml  ketamine  HCl  (100  mg/ml;  Aveco,  Fort  Dodge,  lA)  when  the  tumors  reached 
approximately  1  cm  in  diameter  and  maintained  under  general  gaseous  anesthesia  using  a  small  animal  anesthesia  unit 
with  air  (1  dmVmin)  and  1%  isoflurane  through  a  mask  placed  over  the  mouth  and  nose.  During  the  experiments,  the  rat 
was  placed  on  a  warm  blanket  to  maintain  body  temperature,  which  was  monitored  with  a  rectally  inserted  thermal  probe 
connected  to  a  digital  thermometer  (Digi-Sense,  model  91100-50,  Cole-Parmer  Instrument  Company,  Vernon  Hills,  IF). 
Tumors  were  shaved  before  measurements  to  improve  optical  contact  for  transmitting  light.  A  pulse  oximeter  (model: 
8600V,  Nonin,  Inc.)  was  placed  on  the  hind  foot  to  monitor  arterial  oxygenation  (Sa02)  and  heart  rate.  For  the  single 
high  dose  group  (n=5),  a  light  source  and  four  detectors  from  a  multi-channel,  CW  (continuous  wave)  MRS  (MM,  Inc, 
Philadelphia,  PA)  were  attached  to  the  tumor  using  posts  and  swivel  post  clamps  (see  Figure  1(a)).  For  the  multi  low 
dose  (n=3)  and  control  (n=3)  groups,  we  have  used  four-channel,  frequency  domain  (FD),  MRS  (ISS,  Champaign,  IF). 
In  the  latter  case,  the  four  sets  of  light  sources  replaced  four  detectors  shown  in  Figure  1(b),  and  one  detector  fiber  was 
placed  on  the  top  center  of  the  tumor  to  obtain  signals  from  four  different  regions  of  tumor. 
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Top  View 

(b) 


Fig.  1.  Schematic  diagram  of  experimental  setup  for  the  multi¬ 
channel  NIR  spectroscopy  system,  (a)  A  temperature  probe  was 
located  rectally,  and  pulse  oximeter  was  placed  on  the  hind  paw  to 
monitor  the  arterial  oxygen  saturation,  (b)  A  light  source  probe  was 
placed  on  the  top  center  of  the  tumor,  and  four  detectors  were 
attached  on  four  sides  of  the  tumor  to  detect  the  regional  signals 
during  oxygen  intervention,  (multi-channel  CW  NIRS).  In  the 
setup  of  4  channel  FD  NIRS  system,  the  four  sets  of  light  sources 
replaced  four  detectors  and  one  detector  fiber  was  placed  on  the 
top  center  of  the  tumor  to  obtain  signals  from  four  different  regions 
of  tumor. 


All  the  measurements  were  performed  in  a  dark  room,  and  the  measurements  were  initiated  while  the  rats  breathed  air  for 
10  minutes  to  get  a  stable  baseline.  After  10  minutes  of  baseline  measurement,  the  inhaled  gas  was  switched  to  oxygen 
for  15  minutes,  and  then  back  to  air  for  15  minutes.  This  experimental  procedure  was  repeated  before  and  after  the  CTX 
treatment,  and  the  four  detectors  were  intended  to  be  located  at  the  same  positions  for  each  measurement  on  different 
days.  Using  an  ellipsoidal  approximation,  tumor  volume  V  (cm^)  was  estimated  as  V  =  (7r/6)  L  W  H  where  L,  W,  and  H 
are  the  three  respective  orthogonal  dimensions.  Raw  amplitude  data  from  four  locations  were  recorded  simultaneously 
during  the  experiments  and  processed  after  the  experiments  to  obtain  the  A[Hb02]  and  A[Hbtotai]-  The  amplitude  and 
time  constant  of  A[Hb02]  were  calculated  by  fitting  the  two- exponential  model  to  the  data  using  Kaleidagraph  (Synergy 
Software,  Reading,  PA).  The  corresponding  ratios  of  vasculature  coefficients  and  perfusion  rates,  i.e.,  71/72  dindfi/f2, 
were  also  calculated  to  show  static  and  dynamic  heterogeneities  of  the  tumors  at  different  locations. 


2.2  Measurement  System 

As  mentioned  above,  we  have  used  two  NIRS  systems  for  this  study:  a  multi-channel,  CW,  NIRS  system  was  used  in  the 
single  high  dose  study,  while  an  ISS,  4-channel,  frequency-domain,  NIRS  system  was  applied  in  the  continuous  low  dose 
study.  The  original  CW  NIRS  system  has  the  ability  to  measure  light  signals  from  eight  different  locations,  but  due  to  the 
finite  tumor  size,  we  used  only  four  detectors  to  monitor  tumor  vascular  oxygenation  dynamics  during  respiratory 
challenges.  In  common  with  our  previous  work^^’^®,  we  assume  that  oxyhemoglobin  and  deoxyhemoglobin  are  the  only 
significant  absorbing  materials  in  the  blood-perfused  tumor  tissue  in  the  NIR  range.  The  absorption  coefficients 
comprise  the  extinction  coefficients  for  deoxyhemoglobin  and  oxyhemoglobin  multiplied  by  their  respective 
concentrations  (Eqs.  1  and  2)  at  730  nm  and  850  nm,  that  have  been  employed  in  our  multi-channel  NIRS  system. 


Fa 

|.la' 


730 

850 


=  8Hb^^lHb]  +  8Hb02^^1Hb02], 
=  eHb''°[Hb]  +  eHb02''°[Hb02] . 


(1) 

(2) 
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Based  on  modified  Beer-Lambert’s  law,  the  data  presented  in  this  paper  were  analyzed  using  amplitude  values  to  find  the 
changes  in  absorption  (Eq.  3).  By  manipulating  Equations  1-3,  changes  in  oxygenated  hemoglobin,  deoxygenated 
hemoglobin  and  total  hemoglobin  concentrations  were  calculated  from  the  transmitted  amplitude  of  the  light  through  the 
tumor  (Eqs.  4,  5  and  6). 


|.laB  -  ^aT  =  log  (Ag/A^)  /  L,  (3) 

A[Hb02]  =  [-0.674*log  +1.1 17*log  (Ag/A^)*^®]  /  L,  (4) 

A[Hb]  =  [0.994*log  (Ag/A^)’^°  -  0.376*log  (Ag/A^)*^°]  /  L,  (5) 

A[Hb,^^^,]  =  A[Hb]  +  A[HbO2]=[0.32*log  (Ag/A^)''°  +  0.741  *log  (Ag/A^)''°]  /  L,  (6) 


where  Ab  =  baseline  amplitude;  At  =  transition  amplitude;  L  =  optical  pathlength  between  source/detector.  The 
constants  contained  in  these  equations  were  computed  with  the  extinction  coefficients  for  oxy  and  deoxyhemoglobin  at 
the  two  wavelengths  used.^^  Notice  that  these  coefficients  have  accounted  for  four  hemes  per  hemoglobin  molecule.  In 
principle,  L  should  be  equal  to  the  source-detector  separation,  d,  multiplied  by  a  differential  pathlength  factor  (DPE), 
i.e.,  L=d*DPE.  Little  is  known  about  DPP  for  tumors,  although  a  DPP  value  of  2.5  has  been  used  by  others. Since  our 
focus  is  on  dynamic  changes  and  relative  values  of  tumor  [Hb02]  in  response  to  oxygen  intervention,  we  have  taken  the 
approach  of  including  the  DPP  in  the  unit,  and  eq.  (4)  becomes  as  follows: 

A[Hb02]  =  [-0.674*log  (Ag/A^)’^°  +1.1 17*log  (Ag/A^)*^°]  /  d,  (7) 

where  d  is  the  direct  source-detector  separation  in  cm,  and  the  unit  of  A[Hb02]  in  Eq.  (7)  is  mM/DPP. 

Since  the  wavelengths  of  light  sources  from  the  ISS,  frequency-domain  system  were  750  nm  and  830  nm,  the 
corresponding  equations  for  A[Hb02],  A[Hb],  and  A[Hbtotai]  are  modified  as  follows. 

750  830 

A[Hb02]  =  [-0.709*log  (A^/A^)  -r  1.404*log  (A^/A^)  ]  /  d,  (8) 

750  830 

A[Hb]  =  [0.9546*log  (A^/A^)  -  0.4992*log  (A^/A^)  ]  /  d,  (9) 

750  830 

=  A[Hb]  -r  A[HbO2]=[0.2456*log  (A^/A^)  -r  0.9048*log  (AJA^)  ]  /  d.  (10) 

2.3  Bi-exponential  Model  of  Tumor  Vascular  Oxygenation 

In  our  previous  report, we  followed  an  approach  used  to  measure  regional  cerebral  blood  flow  (rCBP)  with  diffusible 
radiotracers,  as  originally  developed  by  Kety^^  in  the  1950’s.  By  applying  Pick’s  principle  and  defining  j^as  the  ratio  of 
Hb02  concentration  changes  in  the  vascular  bed  to  that  in  veins,  we  arrived  at  Eq.  (11): 

AHbO2™*“'“”70  =  ^[l-exp(-^/Y)]  =  Ai[l-exp(-?/x)]  (11) 

where  ^is  the  vasculature  coefficient  of  the  tumor.  Ho  is  the  arterial  oxygenation  input  and/is  the  blood  perfusion  rate. 

If  a  tumor  has  two  distinct  perfusion  regions,  and  the  measured  signal  results  from  both  of  the  regions  (Pigure  2),  then  it 
is  reasonable  to  include  two  different  blood  perfusion  rates,/i  and/2,  and  two  different  vasculature  coefficients,  and 
in  the  model.  Therefore,  Eq.  (11)  can  be  modified  to  count  for  the  double  exponential  feature  observed  in  the 
experiments: 

AHbO2™*“'“”70  =  riHo[l-exp(-/i?/ri)]  +  )iHo[l-exp(-/2l/)i)] 

=  Ai[l-exp(-?/xi)]  +  A2[l-exp(-?/x2)]  (12) 

where  /  and  are  the  blood  perfusion  rate  and  vasculature  coefficient  in  region  1  for  the  well  perfused  region, 
respectively;  /2  and  have  the  same  respective  meanings  in  region  2  for  the  poorly  perfused  region,  and  Ai  =  A2  = 
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^Ho,  Xi  =  Yilfi,  X2  =  JI//2.  Then,  if  Ai,  A2,  Xi,  and  X2  are  determined  by  fitting  the  measurements  with  the  model,  we  can 
obtain  the  ratios  of  two  vasculature  coefficients  and  the  two  blood  perfusion  rates: 

_  A  ^  /i  _  A /A 

Yi  4  ’  fi  ^1/^2 

With  these  two  ratios,  we  are  able  to  understand  more  about  tumor  physiology,  such  as  tumor  vasculature  and  blood 
perfusion. 


3.  RESULTS 


3.1  Body  weight  and  tumor  volume  changes  during  chemotherapy 

Body  weight  and  tumor  volume  were  monitored  before  and  after  the  CTX  treatments  to  see  the  tumor  responses  and  side 
effects  from  chemotherapy.  In  the  single  high  dose  treatment  group,  body  weight  decreased  until  6  days  after  the 
treatment,  but  later  increased  for  the  rest  of  days  of  observation.  Two  rats  among  five  in  this  group  failed  to  survive  at 
day  6  due  to  the  toxicity  from  the  high  dose  CTX  treatment.  Therefore,  the  data  shown  at  day  8  and  10  represent  the 
smaller  group  of  rats  which  survived  during  the  high  dose  treatment.  Tumor  volume  did  not  further  decreased  after  day 
4.  (Fig.  2(b))  In  comparison,  rats  in  the  continuous  low  dose  group  initially  lost  weight  after  a  low  dose  of  CTX 
administration,  but  gradually  gained  the  weight  during  the  treatment,  presenting  low  toxicity  from  the  treatment.  This 
group  also  showed  a  significant  reduction  in  tumor  volume  during  the  treatment.  For  a  control  group,  saline  was  injected 
into  the  rats  instead  of  CTX,  and  a  gradual  decrease  of  body  weight  was  observed,  while  their  tumor  volumes  increased 
exponentially.  Changes  in  rat  body  weight  and  tumor  volume  were  normalized  to  day  0  (before  CTX  or  saline 
administration).  (Figures  2(a)  and  (b))  Solid  circles  represent  the  data  from  a  control  group,  and  open  squares  and  open 
diamonds  represent  the  continuous  low  dose  group  and  single  high  dose  group,  respectively. 


Fig.  2.  Normalized  changes  in  rat  body  weight  (a)  and  tumor  volume  (b)  during  the  saline  and  CTX  treatments. 


3.2  Intratumoral  heterogeneity  of  vascular  oxygenation  observed  by  the  multi-channel  NIRS 

In  our  experiment,  we  have  utilized  either  one  light  source  and  four  detectors  and  from  the  CW  system  (NIM,  Inc)  or 
four  light  sources  and  one  detector  from  the  FD  system  (ISS,  Inc.)  for  the  rat  breast  tumor  measurements,  and  the  setups 
(  Fig.  1).  After  10  minutes  of  baseline  measurement  with  air  breathing,  gas  was  switched  to  pure  oxygen,  causing  a  rapid 
increase  in  tumor  [Hb02].  These  changes  were  measured  simultaneously  from  four  locations  of  the  tumor.  Figure  3 
shows  a  representative  set  of  data  before  the  CTX  treatment,  with  the  DC  NIRS  system.  Open  circles  show  the  raw  data 
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measured  by  multi-channel  NIRS,  and  the  solid  black  lines  represent  the  fitted  curves  using  our  bi-exponential  model  for 


hemodynamics  during  oxygen  intervention, 
similar  trends  of  A[Hb02]. 


It  is  apparent  that  the  data  from  each  location  differs  though  there  are 
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Fig.  3.  Dynamic  changes  of  tumor  [Hb02]  from  four  locations  in  a  rat  breast  tumor.  The  rising  parts  of  A[Hb02]  from  the  four 
locations  were  fitted  using  a  double-exponential  expression.  Figures  3(a)-3(d)  were  taken  from  locations  #l-#4,  respectively.  In  this 
case,  the  tumor  was  not  treated  yet. 

Table  1.  Summary  of  vascular  oxygen  dynamics  determined  at  the  four  detectors  from  tumor  shown  in  the  Fig.  3. 


Location 

Ai 

A2 

Xi 

X2 

A1/A2  (=71/72) 

T1/X2 

Ml 

1  (Fig.  3a) 

0.013 

0.027 

0.48 

24.3 

0.48 

0.02 

24 

2  (Fig.  3b) 

0.029 

0.026 

0.09 

29.8 

1.12 

0.003 

373 

3  (Fig.  3c) 

0.008 

0.036 

0.15 

10.6 

0.22 

0.014 

16 

4  (Fig.  3d) 

0.021 

0.037 

0.12 

7.25 

0.57 

0.017 

34 

To  compare  the  data  taken  from  four  locations  of  the  tumor  more  clearly,  the  time  constants  and  amplitudes  from  the 
four  fitted  curves  are  summarized  in  Table  1.  The  ratios  of  71/72  dindfi/f2  characterize  tumor  vascular  structure  and  blood 
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perfusion  within  the  volume  of  tumor  interrogated  by  light.  In  principle,  when  71/72  is  close  to  1,  it  implies  that  the 
measured  optical  signal  results  equally  from  both  regions  1  {Le.,  well  perfused  region)  and  2  (i.e.,  poorly  perfused 
region);  if  71/72  <  1,  the  measured  signal  results  more  from  region  2  than  region  1  [Figures  3(a),  3(c)  and  3(d)].  As  Table 
1  demonstrates,  only  location  #2  has  a  ratio  of  71/72  slightly  higher  than  1,  and  the  readings  from  locations  #1,  #3  and  #4 
have  the  ratios  of  71/72  less  than  1.  This  may  suggest  that  the  tumor  volume  that  was  optically  interrogated  from  location 
#2  was  dominated  by  well  perfused  regions,  while  most  of  other  tumor  volumes  detected  from  locations  #1,  #3  and  #4 
are  composed  of  more  poorly  perfused  regions.  Furthermore,  all  the  ratios  of  /1//2  from  four  locations  of  the  tumor 
shown  in  Fig.  3  are  much  greater  than  1,  indicating  that  the  blood  perfusion  rate  in  well  perfused  region  is  much  greater 
than  that  in  poorly  perfused  region.  Especially,  fi/f2  from  location  #2  is  10  to  20  times  higher  than  those  from  locations 
#1,  #3,  and  #4,  showing  a  high  level  of  intratumoral  heterogeneity  in  dynamic  vascular  structure. 

3.3  Monitoring  vascular  hemodynamics  of  breast  tumors  before  and  after  chemotherapy 

The  tumor  hemodynamics  during  oxygen  intervention  were  measured  before  and  after  administration  of  CTX  and  saline. 
The  representative  data  from  the  control  group  and  continuous  low  dose  group  are  shown  in  Figs.  4(a)  and  4(b), 
respectively.  Similar  to  Figure  3,  open  symbols  are  the  raw  data  from  measurements,  and  solid  lines  are  the  fitted  curves 
using  our  double  exponential  model.  As  mentioned  before,  there  were  4  light  sources  placed  on  the  surface  of  tumor. 
Figure  4  shows  the  acute  and  then  gradual  changes  of  [Hb02]  after  switching  the  breathing  gas  from  air  to  oxygen,  and 
the  data  were  observed  at  the  same  (or  nearly  the  same)  location  of  tumor  from  day  0  to  day  6.  From  Figure  4(a),  we  can 
see  different  tumor  hemodynamics  at  different  days,  but  having  similar  trends  and  maximum  A[Hb02]  for  the  control 
group.  However,  the  data  taken  from  the  continuous  low  dose  group  show  quite  different  hemodynamics  throughout  the 
treatment  days  (Fig.  4(b)).  Especially,  we  notice  that  the  fast  increase  part  became  much  smaller  at  Days  2  and  6 
compared  to  Day  0,  implying  a  significant  decrease  in  signal  from  the  well  perfused  region. 


Fig.  4.  Dynamic  changes  of  A??[Hb02]  taken  at  location  #1  from  a  rat  breast  tumor  before  and  after  administration  of  (a)  saline  and 
(b)  continuous  low  dose  of  CTX  (20  mg/kg  for  10  days).  The  rising  part  of  A[Hb02]  from  location  #1  was  fitted  using  the  double¬ 
exponential  expression. 

Figure  5  also  shows  the  changes  of  tumor  hemodynamics  during  oxygen  intervention,  before  and  after  a  single  high  dose 
of  CTX  treatment,  measured  by  the  CW  NIRS.  This  figure  clearly  demonstrates  that  we  can  observe  significant  changes 
in  tumor  hemodynamics  after  chemotherapy  by  using  respiratory  challenge  as  a  mediator.  The  fitted  parameters  from 
our  mathematical  model  are  summarized  in  Table  2  to  compare  the  changes  in  hemodynamic  parameters  before  and  after 
administration  of  CTX.  The  rising  part  of  A[Hb02]  from  location  #1  was  fitted  using  either  a  single-  or  double¬ 
exponential  expression. 
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Fig.  5.  Dynamic  changes  of  [Hb02]  taken  at  location  #1  from  a  rat  breast  tumor  before  and  after  a  single  high  dose  of  CTX  treatments 
(200mg/kg).  The  rising  part  of  A[Hb02]  from  location  #1  was  fitted  using  either  a  single-  or  double-exponential  expression. 


Table  2.  Summary  of  vascular  oxygen  dynamics  determined  at  location  #1  from  the  tumor  shown  in  Fig.  5 

before  and  after  CTX  treatment. 


Day 

Ai 

A2 

Xi 

X2 

A1/A2  (=71/72) 

X1/X2 

Ml 

0 

0.044 

0.031 

0.23 

25.21 

1.42 

0.0091 

156 

1 

0.032 

0.033 

0.13 

11.36 

0.97 

0.0114 

85 

3 

0.0087 

0.014 

0.089 

8.36 

0.62 

0.0106 

58 

5 

0.0087 

1.27 

4.  DISCUSSION 

Changes  of  rat  body  weight  and  tumor  volume  clearly  show  that  CTX  treatment  is  effective  for  the  tumor  type  that  we 
have  used  in  this  study.  For  the  control  group,  the  average  rat  body  weight  gradually  decreased  during  the  entire  course 
of  treatment,  which  implies  the  sickness  of  rats  possibly  due  to  the  tumor  growth  (cachexia).  (At  Day  10,  the  tumor 
volume  was  ~5  times  larger  than  that  at  Day  0.)  It  is  clear  that  there  is  a  different  effectiveness  of  CTX  treatment 
between  the  single  high  dose  group  (200mg/kg)  and  continuous  low  dose  group  (20mg/kg  for  10  days).  Both  of  the  CTX 
treatments  delayed  the  tumor  growth  and  even  further  reduced  the  tumor  volume.  However,  a  single  high  dose  of  CTX 
treatment  caused  the  death  of  two  rats,  and  the  tumor  volume  was  not  decreased  further  4  days  after  the  treatment,  while 
the  continuous  low  dose  CTX  treatment  continued  to  provide  tumor  regression  without  causing  severe  sickness.  From 
this  observation,  it  is  obvious  that  continuous  low  dose  of  CTX  treatment  is  working  much  better  than  a  single  high  dose 
of  CTX  treatment  for  a  rat  mammary  adenocarcinomas  13762NF  tumor. 

MRS  is  a  portable,  low  cost,  and  real  time  measurement  system  that  can  monitor  changes  of  vascular  oxygen  levels  by 
using  two  wavelengths.  We  have  previously  used  a  single-channel  MRS  system  with  one  light  source  and  one  detector 
for  global  measurements  of  A[Hb02]  and  A[Hbtotai]  in  tumors  during  respiratory  challenges. Through  those 
experiments,  we  have  found  that  most  tumors  have  a  bi-phasic  behavior  in  A[Hb02]  increase  (i.e.,  a  rapid  increase 
followed  by  a  slow  and  gradual  increase)  after  switching  the  gas  from  air  to  carbogen/oxygen.  To  explain  this  bi-phasic 
behavior,  we  developed  a  mathematical  model  and  formed  a  hypothesis  that  the  bi-phasic  behavior  of  A[Hb02]  during 
carbogen/oxygen  inhalation  results  from  two  different  vascular  regions  in  tumor  with  two  blood  perfusion  rates  and 
vascular  structures. 
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By  giving  an  oxygen  intervention,  tumor  blood  vessels  are  acutely  subject  to  an  increase  of  [Hb02]  due  to  higher  supply 
of  oxygenated  blood  from  artery  compared  to  that  from  air  breathing.  However,  due  to  the  irregular  vascular  structure  in 
tumor,  the  well  perfused  regions  in  tumor  may  have  an  increase  in  [Hb02]  much  faster  than  other  parts  of  tumor  that  are 
poorly  perfused.  Therefore,  two  time  constants  obtained  from  tumor  hemodynamic  measurements  during  oxygen 
intervention  are  able  to  reveal  two  blood  flow/perfusion  rates  in  tumor,  more  precisely,  two  speeds  of  blood  flow  within 
the  tumor  blood  vessels.  More  recently,  we  have  shown  that  the  bi-phasic  increase  in  optical  density  changes  occurs 
when  there  exist  two  different  flow  rates  in  tumor  vascular  phantom.^^  The  bi-phasic  model  is  a  basis  of  our  current 
study  where  we  wish  to  detect  any  changes  in  vascular  structures,  hemodynamic  features,  or  perfusion  rate  within  a 
tumor  after  CTX  treatment. 

The  amplitude  and  time  constants  obtained  from  A[Hb02]  increase  (Fig.  5)  are  summarized  in  Table  2.  At  day  0,  we  can 
see  that  71/72  is  higher  than  1,  indicating  that  the  measured  signal  results  more  from  the  well  perfused  region  than  poorly 
perfused  region.  However,  this  ratio  becomes  less  than  1  after  injection  of  cyclophosphamide  (Day  1  and  3).  This  may 
be  explained  by  destruction  of  vascular  structure  in  tumor  after  chemotherapy.  We  expect  that  after  a  single  high  dose 
administration  of  CTX,  the  drug  circulates  in  the  blood  vessels  and  is  delivered  to  the  tumor  cells  more  in  the  well 
perfused  region  than  in  the  poorly  perfused  region.  This  will  lead  to  death  of  tumor  cells  in  the  well  perfused  region 
more  effectively  than  that  in  the  poorly  perfused  region,  eventually  resulting  in  decreases  in  tumor  volume  in  the  well 
perfused  region  more  than  in  the  poorly  perfused  region.  Then,  the  tumor  volume  containing  the  well  perfused  regions 
will  consequently  decrease,  so  will  the  contribution  of  detected  NIR  signals  from  the  well  perfused  region.  In  other 
words,  a  decrease  in  71/72  may  indicate  deceases  in  well  perfused  regions  in  tumor  volume,  after  the  administration  of 
CTX. 

As  shown  in  Table  2,  moreover,  the  perfusion  rate  ratio,  /1//2,  was  also  decreased  after  a  single  high  dose  of  CTX 
administration.  At  Day  0,  /1//2  was  very  high,  meaning  that  there  was  a  big  difference  of  perfusion  rate  between  the  well 
perfused  and  poorly  perfused  region  in  tumor.  However,  this  ratio  significantly  decreased  at  Day  1  and  3  after  CTX 
treatment,  representing  that  the  perfusion  rate  gap  between  the  well  perfused  region  and  poorly  perfused  region  became 
much  smaller  than  that  at  day  0.  At  Day  5,  changes  in  [Hb02]  during  oxygen  intervention  do  not  show  any  bi-phasic 
behavior  anymore,  and  it  was  fitted  by  a  single-exponential  model.  This  may  indicate  that  most  of  tumor  cells  and/or 
tumor  vasculature  in  the  well  perfused  region  are  possibly  destroyed  by  the  effect  of  CTX,  resulting  in  quite  different 
hemodynamic  behavior. 


5.  CONCLUSION 

In  conclusion,  we  have  conducted  this  study  to  show  the  possibility  of  using  NIRS  to  monitor  tumor  hemodynamics  in 
response  to  chemotherapy  by  comparing  the  changes  in  tumor  vascular  oxygenation  before  and  after  CTX  treatment. 
The  heterogeneity  of  tumor  vasculature  was  easily  observed  by  quantifying  the  blood  perfusion  rate  and  vascular 
coefficients  at  four  different  locations  of  the  tumor.  Tumor  hemodynamics  has  been  significantly  changed  before  and 
after  CTX  treatment  compared  to  the  saline-treated  control  group,  showing  high  possibility  of  the  NIRS  system  to  be 
used  as  a  monitoring  tool  for  cancer  treatments.  Our  future  studies  will  include  the  development  of  NIR  imaging  systems 
to  obtain  a  map  of  tumor  hemodynamic  changes  from  whole  tumor,  allowing  us  to  predict  the  efficacy  of  tumor 
treatment. 
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results  exhibit  a  linear  correlation  between  AHbO  and  Ap02  of  the  tumors 
under  hyperoxic  gas  intervention,  suggesting  that  the  NIRS  approach  could 
have  a  good  potential  value  in  the  clinic.  Finally,  the  newly  developed 
tumor  hemodynamic  model  allows  us  to  reveal  tumor  heterogeneities  at 
different  tumor  locations  based  on  the  multichannel  NIRS  results.  Through 
this  chapter,  we  lay  a  foundation  for  an  NIR  imaging  technique  to  be  fur¬ 
ther  developed  to  facilitate  investigations  of  tumor  heterogeneity  and  vas¬ 
cular  perfusion.  Such  a  noninvasive  imaging  approach  can  enhance  our 
understanding  of  the  dynamics  of  tumor  oxygenation  and  the  mechanism 
of  tumor  physiology  under  baseline  and  perturbed  conditions. 
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Introduction 

Significance  of  pO 2  in  Oncology 

It  has  long  been  appreciated  that  hypoxic  tumor  cells  are  more  resistant 
to  radiotherapy.^  Indeed,  a  3-fold  increase  in  radio  resistance  may  occur 
when  cells  are  irradiated  under  hypoxic  conditions  compared  with  oxygen 
pressure  p02  >  15  torr  for  a  given  single  radiation  dose.  However,  recent 
modeling  has  indicated  that  the  proportion  of  cells  in  the  range  0-20  torr 
may  be  most  significant  in  terms  of  surviving  a  course  of  fractionated  radio¬ 
therapy.^  Certain  chemotherapeutic  drugs  also  present  differential  efficacy, 
depending  on  hypoxia.^’"^  Increasingly,  there  is  evidence  that  hypoxia  also 


^  L.  Gray,  A.  Conger  et  al.,  Br.  J.  Radiol  26,  638  (1953). 

^  B.  G.  Wouters  and  J.  M.  Brown,  Radial  Res.  147,  514  (1997). 
^  B.  Teicher,  J.  Lazo  et  al,  Cancer  Res.  41,  73  (1981). 

^  A.  C.  Sartorelli,  Cancer  Res.  48,  775  (1988). 
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influences  such  critical  characteristics  as  angiogenesis,  tumor  invasion,  and 
metastasis.^”^  Moreover,  repeated  bouts  of  intermittent  hypoxic  stress  may 
be  important  in  stimulating  tumor  progression.^  Thus  the  ability  to  mea¬ 
sure  p02  noninvasively  and  repeatedly,  with  respect  to  acute  or  chronic 
interventions,  becomes  increasingly  important. 

Early  work  examined  cells  in  vitro,  where  ambient  oxygen  concentra¬ 
tions  are  readily  controlled.  In  vivo,  hypoxia  may  be  achieved  by  clamping 
the  blood  supply  to  a  tumor,^^  but  other  levels  of  oxygenation  reflect  the 
interplay  of  supply  and  consumption.^^’^^  Robust  fine-needle  polarographic 
electrodes  opened  the  possibility  of  measuring  p02  in  tumors  in  situ  and 
in  vivo  to  define  local  p02  under  baseline  conditions  or  with  respect  to 
interventions.  In  early  work.  Cater  and  Silver^^  showed  the  ability  to  moni¬ 
tor  p02  at  individual  locations  in  patients’  tumors  with  respect  to  breathing 
oxygen.  Later,  Gatenby  et  al}"^  showed  that  p02  in  a  tumor  was  correlated 
with  clinical  outcome.  Tumor  oximetry  received  its  greatest  boost  with  the 
development  of  the  Eppendorf  Histograph  polarographic  needle  electrode 
system.  This  computer-controlled  device  equipped  with  a  stepper  motor 
can  reveal  distributions  of  tumor  oxygenation  and  has  been  applied  exten¬ 
sively  to  clinical  trials.  Many  reports  have  now  shown  that  tumors  are 
highly  heterogeneous  and  have  extensive  hypoxia;  furthermore,  strong 
correlations  have  been  shown  in  cervix  and  head  and  neck  tumors  between 
median  p02  or  hypoxic  fraction  and  survival  or  disease-free  survival.^’^^”^^ 
Extensive  hypoxia  also  has  been  found  in  tumors  of  the  prostate  and 
breast.^^“^^  Thus  tumor  oxygenation  is  now  recognized  as  a  strong 
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prognostic  indicator,  and  this  device  has  laid  a  convincing  foundation 
for  the  value  of  measuring  p02  in  patients.  However,  the  Histograph  is 
highly  invasive,  and  it  is  not  possible  to  make  repeated  measurements  at  in¬ 
dividual  locations,  precluding  dynamic  studies  to  assess  the  influence  of 
interventions  on  tumor  p02. 

Given  that  hypoxic  tumors  are  more  resistant  to  certain  therapies,  it  be¬ 
comes  important  to  assess  tumor  oxygenation  as  part  of  therapeutic  plan¬ 
ning.  Patients  could  be  stratified  according  to  baseline  hypoxia  to  receive 
adjuvant  interventions  designed  to  modulate  p02,  or  more  intense  therapy 
as  facilitated  by  intensity  modulated  radiation  therapy  (IMRT).  Tumors, 
which  do  not  respond  to  interventions,  may  be  ideal  candidates  for  hyp¬ 
oxia-selective  cytotoxins  (e.g.,  tirapazamine^"^).  Noting  that  any  therapy 
and  intervention  may  have  side  effects  or  simply  add  to  clinical  costs,  it  is 
vital  that  efficacy  be  established  and  therapy  be  optimized  for  an  individual 
patient.  Whether  initially  hypoxic  regions  of  a  tumor  can  be  modified  to 
become  better  oxygenated  has  long  been  considered  a  key  to  improving 
outcome  of  irradiation.  However,  many  attempts  to  improve  therapeutic 
outcome  by  manipulation  of  tumor  oxygenation  have  shown  only  modest 
success  in  the  clinic,^^  and  it  is  thought  that  lack  of  success  may  have 
resulted  from  inability  to  identify  those  patients  who  would  benefit  from 
adjuvant  interventions. 

Although  p02  determinations  could  be  of  great  clinical  value,  they  are 
also  vital  to  many  laboratory  investigations  of  new  drugs  and  studies  of 
tumor  development.  Given  the  potential  importance  of  measuring  p02, 
many  diverse  techniques  have  been  developed,  as  reviewed  by  others 
previously and  here,  in  the  next  section. 

Methods  of  Measuring  Tumor  Oxygenation 

Table  I  lists  various  techniques  that  have  been  reported  to  provide 
quantitative  estimates  of  p02.  Historically,  polarographic  needle  oxygen 
electrodes  have  been  considered  a  “gold  standard,”  and  they  have  been  ap¬ 
plied  in  the  clinic  since  the  1950s.  One  or  more  electrodes  may  be  placed  in 
a  tumor,  facilitating  measurement  of  baseline  p02  and  dynamic  response  to 
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interventions.  ’  “  Initially,  the  focus  was  on  generating  finer  needles, 
which  would  be  less  invasive,  and  tips  as  fine  as  a  few  microns  have  been 
applied  to  animal  tissues.^^  However,  such  needles  are  progressively  brittle 
and  generate  such  small  current  that  stray  electromagnetic  fields  can  inter¬ 
fere.  Stationary  electrodes  sample  limited  volumes,  and  recognizing  tumor 
heterogeneity,  the  Eppendorf  Histograph  was  developed  to  generate  mul¬ 
tiple  measurements  along  tracks  in  tumors. Following  extensive  stud¬ 
ies  in  animals,  the  Histograph  has  found  widespread  application  in  the 
clinical  setting  and  has  unequivocally  revealed  hypoxia  in  many  tumor 
types,  for  example,  head  and  neck,^^’^^  cervix,^^’^^  breast,^^’^^  and  pros- 
tate.^^  Moreover,  p02  distributions  have  been  found  to  have  prognostic 
value.  Disease-free  survival  is  significantly  worse  for  patients  with  hypoxic 
tumors,  though  the  optimal  prognostic  parameter  has  variously  been 
median  p02  or  percent  measurements  <5  torr  (HF5). 

Although  the  Eppendorf  Histograph  uses  a  large  invasive  needle  (size  = 
26  G  or  about  0.35  mm),  it  has  provided  great  impetus  for  further  investiga¬ 
tions.  One  aspect  is  the  application  of  less  invasive  probes.  Fiber-optic 
probes  are  typically  finer  and  do  not  consume  oxygen  during  measurement. 
Typically,  only  two  or  four  locations  are  sampled  simultaneously,  but  as  with 
the  earlier  electrodes,  these  optical  probes  facilitate  observation  of  dynamic 
changes  in  p02  in  response  to  interventions. Both  the  current  commer¬ 
cial  systems,  the  OxyLite  (http://www.oxford-optronix.com/tissmon/oxylite/ 
oxylite.htm)  and  FOXY  (http://www.oceanoptics.com/Products/foxyfaqs. 
asp),  exploit  the  fluorescent  quenching  by  oxygen  of  a  ruthenium  complex 
coating.  OxyLite  measures  fluorescent  lifetime,  whereas  FOXY  uses  a 
simple  intensity  integration  and  is  correspondingly  much  cheaper.  Fibers 
are  fragile,  and  coatings  have  a  limited  lifetime. 
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TABLE  I 

Tumor  Oximetry  Methods 


Characteristic  resolution 

Parameter  - 

Technique  Reporter  measured  Invasiveness  Spatial  Temporal  References 


FREDOM 

HFB 

Ri 

Minimal 

32-G 

needle 

Map  multiple 
locations 
each  8  mm^ 

6.5  min 

Hunjan,  Zhao  et 

Zhao,  Constantinescu  et  alP\ 
Song,  Constantinescu  et 

Zhao,  Constantinescu  et  al?^\ 
Kim,  Zhao  et 

Zhao,  Constantinescu  et 

Zhao,  Ran  et  alP^ 

MRI 

PFC 

Ri 

IV 

Perfused  regions 

min 

Hees  and  Sotak^^^; 

Dardzinski  and  Sotak^^^; 
McIntyre,  McCoy  et 

Fan,  River  et  alP; 

Wang,  Su  et  alP 

^‘^F  MRS 

PFC 

Ri 

IV 

Perfused  regions 

s  to  min 

Hees  and  Sotak^^^; 

Mason,  Antich  et 

Baldwin  and  Ng^^^; 

McIntyre,  McCoy  et 
van  der  Sanden,  Heerschap 
et  alP^ 

DCE  MRI 

Gd-DTPA 

Contrast 

kinetics 

IV 

Maps 

min 

Cooper,  Carrington  et  alP\ 

Lyng,  Vorren  et  alf"', 

Wang,  Su  et  alP 

ESR/EPR 

Charcoal, 

phthalocyanine 

Linewidth 

Needle 

Single  location 

23  G  IT 

seconds 

O’Hara,  Goda  et  alP, 

Goda,  Bacic  et  alf; 

O’Hara,  Goda  et  alP; 

Gallez,  Jordan  et  alp 
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ESR/EPR 

Nitroxides 

Linewidth 

IV 

OMRI 

Free  radical 

Overhauser 

enhancement 

IV 

Needle 

electrode 

Oxygen 

Current 

Needle 
26  G 

(Histograph) 

Oxygen 

Current 

Needle 
26  G 

Global  or  map 
mm^  maps 

Single  location 

Multiple  tracks 


Jordan,  Misson  et 
He,  Begliein  et  al^^; 

Jiang,  Beghei  et 
O’Hara,  Blumenthal  et  alf; 
Baudelet  and  Gallez^"^; 

Dunn,  O’Hara  et  alf; 

Jordan,  Gregoire  et 
Matiy,  De  Bast  et  alP 
s  to  min  Bias,  Williams  et  alP 
10  min  Krishna,  English  et  alP 

s  Cater  and  Silver^^; 

Evans  and  Naylor^^; 

Hasegawa,  Rhee  et  al.^; 
Gatenby,  Kessler  et  alP; 

Song,  Shakil  et  alp 
Zhao,  Constantinescu  et  alP 
1  s  per  Eble,  Wenz  et  al!\ 
location  Falk,  Laurence  et  alJ; 

Vaupel,  Kelleher  et  alP; 

Brizel,  Scully  et  alP; 

Hockel,  Schlenger  et  alP, 
Nozue,  Lee  et  alP; 

Fyles,  Milosevic  et  aP; 
Siemann,  Johansen  et  al!\ 
Mason,  Constantinescu  et  alP^; 
Aquino-Parsons,  Green  et  alP; 
Jenkins,  Evans  et  alP\ 

Hockel  and  Vaupel^ 
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TABLE  I  (continued) 


Parameter 

measured 

Characteristic  resolution 

Technique 

Reporter 

Invasiveness 

Spatial 

Temporal 

References 

Optical  probe 
(OxyLite, 
FOXY) 

Rh  complex 

Fluorescent 

lifetime 

Needle 
(26  G) 

2-4  locations 

Real  time 

Griffiths^^;  Bussink, 

Kaanders  et  al.^^\ 

Braun,  Lanzen  et  al.^^; 

Zhao,  Constantinescu  et  al.^^\ 
Gu,  Bourke  et  al.^^; 

Jordan,  Beghein  et  al.^^ 

Phosphorescence 

PD  complex 

Lifetime 

IV 

Maps 

<  1  min 

Wilson"^; 

Vinogradov,  Lo  et  a/.”; 

Dewhirst,  Ong  et  al.^; 

Wilson,  Vinogradov  et  al.^^; 
Erickson,  Braun  et  al.^ 

Fluorescence 

EF5 

Fluorescent 

intensity 

IV  +  biopsy 

Maps 

microscopic 

Once 

Koch^^ 

Mass 

spectrometry 

Oxygen 

Atoms 

Needle 

Single 

location 

Potapov,  Sirovskii  et  al.^ 
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Reporter  molecules  have  been  developed  for  use  with  electron  spin 
resonance  (ESR  or  EPR),  where  the  line  width  is  highly  sensitive  to 
oxygen.^^’'^'^“^^  Two  primary  approaches  are  used:  (1)  direct  intratumoral 
(IT)  injection  of  char  crystals, phthalocyanine,^^  or  India  ink^^  into  a 
tissue  or  (2)  intravenous  (IV)  infusion  of  water-soluble  agents  which  dis¬ 
seminate  throughout  the  tumor  vasculature. Direct  IT  injection  is  inva¬ 
sive  and  has  generally  been  applied  as  a  spectroscopic  approach  to  report 
p02  at  single  locations  only.  Nonetheless,  significant  data  have  been 
achieved  demonstrating  hypoxiation  and  reoxygenation  with  respect  to  ir¬ 
radiation,  and  the  importance  of  timing  successive  radiation  doses  to  coin¬ 
cide  with  reoxygenation.^^  Char  particles  may  be  stable  in  tissue  for  weeks 
to  years,  allowing  measurements  of  chronic  changes  in  tissues  (e.g.,  accom¬ 
panying  tumor  growth).^^  The  IV  approach  is  noninvasive,  but  reporter 
molecules  may  predominately  distribute  in  the  well-perfused  vasculature, 
potentially  biasing  measurements  toward  the  well-oxygenated  tumor 
regions.  Progressive  uptake  and  clearance  of  agents  produces  variable  con¬ 
centrations,  and  some  agents  degrade  in  tissue  requiring  appropriate  cor¬ 
rection  factors.^^  Nonetheless,  images  of  tumor  oxygen  distribution  have 
been  reported,  including  three-dimensional  representations.^^  Spin  radicals 
also  may  be  applied  to  a  combined  ESR-NMR  (nuclear  magnetic  res¬ 
onance)  approach,  Overhauser-enhanced  magnetic  resonance  imaging 
(OMRI),  exploiting  the  Overhauser  enhancement  in  the  tissue  water 
proton  MRI  signal  that  occurs  by  polarization  transfer  from  free  radicals 
upon  electromagnetic  irradiation.^^ 
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Vascular  oxygenation  has  been  probed  by  fluorescence  or  phosphores¬ 
cence  imaging  based  on  reporter  complexes  delivered  Historically, 

the  approach  was  limited  to  superflcial  tissues  due  to  limited  light  penetra¬ 
tion.  The  latest  molecules  are  active  in  the  near-infrared,  permitting 
greater  depth  of  signal  penetration.^^ 

NMR  facilitates  interrogation  of  deep  tissues  noninvasively,  and 
NMR  approaches  will  be  reviewed  in  detail  in  the  following  section.  The 
methods  discussed  earlier  provide  direct  quantitative  measurements  of 
p02  based  on  various  physiochemical  parameters,  such  as  electric  current, 
fluorescent  lifetime,  magnetic  resonance  linewidth,  or  relaxation.  Other 
approaches  are  less  direct,  but  can  reveal  hypoxia  or  correlates  of  p02. 

Speciflc  classes  of  reporter  molecules  have  been  developed  to  reveal 
hypoxia^*’’*’^  (e.g.,  pimonidazole,*’^’*’®  EFS,*’’’®®  CC1-103F,®  Cu-ATSM'^“’'^^ 
galactopyranoside  lAZA^^).  Following  IV  infusion,  these  agents  become 
reduced  in  tissues  and  are  trapped.  However,  in  the  presence  of  oxygen 
they  are  reoxidized  and  ultimately  clear  from  the  body.  Histologic  assess¬ 
ment  of  the  distribution  of  these  agents  provides  microscopic  indications 
of  local  hypoxia.  EF5,  pimonidazole,  and  Cu-ATSM  are  currently  being 
tested  in  clinical  trials,  and  correlations  have  been  reported  with  clinical 
outcome.^^’^^’^^  Many  variants  have  been  proposed  over  the  past  20  years, 
and  incorporation  of  radionuclides  has  facilitated  noninvasive  investiga¬ 
tions  using  positron  emission  tomography  (PET)  or  single  photon  emission 
computed  tomography  (SPECT),  while  ^^F  labels  permitted  NMR  spec¬ 
troscopy.  Generally,  only  a  single  time  point  is  investigated,  but  dy¬ 
namic  variations  in  hypoxia  may  be  assessed,  even  in  biopsy  specimens, 
by  applying  pairs  of  hypoxia  reporters  in  a  pulse-chase  fashion  with  respect 
to  an  intervention,  as  shown  by  Ljungkvist  et  al.^^ 
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Several  studies  have  shown  a  lack  of  correlation  between  hypoxic 
marker  binding  and  p02  assessed  using  the  Eppendorf  Histograph,  which 
may  be  related  to  chronic  versus  acute  hypoxia,  or  the  extent  of  necro- 
sis7^“^^  The  ultimate  value  of  the  techniques  is  evidenced  by  correlations 
between  uptake  and  outcome.^^’^^’^^’^^  Recent  data  also  indicate  that  EF5 
fluorescence  may  be  correlated  with  p02.^^ 

The  techniques  discussed  so  far  all  depend  on  exogenous  reporter  mol¬ 
ecules  or  probes.  Ideally,  oxygenation  could  be  related  to  endogenous 
characteristics.  Because  many  biochemical  pathways  are  under  oxygen 
regulation,  they  can  provide  an  elegant  window  on  hypoxia,  for  example, 
induction  of  hypoxia-inducible  factor  1  (HIF-1)  and  glucose  transporter  1 
(Glut-1)  together  with  secondary  responses,  such  as  increased  production 
of  vascular  endothelial  growth  factor  (VEGF),  NIPS  and  tumor-associated 
macrophage  activity.^  Such  molecules  indicate  hypoxia,  though  they  may 
be  induced  by  other  factors.  Intrinsic  radiation  sensitivity  also  may  be 
assessed  using  the  Comet  assay.^^  These  assays  each  require  biopsy.  Other 
markers  potentially  associated  with  hypoxia  may  be  found  in  the  plasma  or 
urine  and  have  been  correlated  with  clinical  outcome.^^  An  attractive  alter¬ 
native  is  the  introduction  of  transgenes  with  hypoxic  response  elements 
(HREs)  as  promoter  sequences  coupled  to  reporter  genes  such  as  GFP 
(green  fluorescent  protein)^^’^^  or  luciferase.^"^’^^  GFP  synthesis  is  an  ener¬ 
getic  process,  which  could  be  hindered  under  hypoxia  conditions.  Likewise, 
bioluminescence  accompanying  action  of  luciferase  on  luciferin  requires 
adenosine  triphosphate  (ATP)  and  O2,  but  reports  suggest  that  even  under 
exceedingly  low  p02,  sufficient  oxygen  remains  to  reveal  hypoxia. 

Many  practical  considerations  govern  clinical  application  of  oximetry 
methods.  Proton  MRI  is  routinely  applied  for  anatomic  evaluation  of 
tumors  and  would  provide  an  ideal  conduit  for  prognostic  investigations. 
Application  of  contrast  agents  may  reveal  tumor  boundaries  to  enhance 
detectability,  and  the  dynamic  contrast  enhancement  (DCE)  changes 
provide  insight  into  vascular  perfusion  and  surface  permeability  area.^^ 
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Specific  studies  have  shown  a  correlation  between  DCE  and  p02,^^  and 
indeed,  a  theoretical  underpinning  has  been  provided  based  on  the  Krogh 
cylinder  model.^^  However,  the  correlation  is  unlikely  to  be  widely  applic¬ 
able,  since  DCE  is  sensitive  to  vascular  flow,  perfusion,  and  permeability, 
where  p02  depends  on  oxygen  consumption  as  well  as  delivery. 

Blood  oxygen  level  dependent  (BOLD)  contrast  proton  NMR  facili¬ 
tates  rapid  interrogation  of  vascular  oxygenation  and  is  particularly 
appropriate  for  examining  dynamic  responses  to  interventions."^^’^^”^^ 
Deoxyhemoglobin  is  paramagnetic  and  induces  signal  loss  in  T2*-weighted 
images.  However,  BOLD  does  not  provide  absolute  p02  values  and  is  con¬ 
founded  by  the  influence  of  blood  flow,  as  investigated  extensively  by 
Howe  et  al.^^  who  termed  the  expression  FLOOD  (flow  and  oxygen  level 
dependent)  contrast.  In  addition,  variation  in  vascular  volume  can  intro¬ 
duce  signal  perturbation.^^  Nonetheless,  some  studies  have  indicated  a 
correlation  with  relative  p02,  but  poor  indication  of  absolute 

Near-infrared  spectroscopy  (NIRS)  offers  an  alternative  approach 
based  on  the  differential  light  absorption  of  the  strong  chromophores  oxy¬ 
hemoglobin  and  deoxyhemoglobin.  NIRS  provides  a  noninvasive  means  to 
monitor  global  tumor  vascular  oxygenation  in  real  time  based  on  endoge¬ 
nous  molecules.  Although  many  NIRS  investigations  have  been  conducted 
in  the  brain  and  breast  in  both  laboratory  and  clinical  settings  over  the 
past  decade,  there  have  been  relatively  few  reports  regarding  solid 
tumors."^^’^^”^^^  Most  studies  to  date  have  used  reflectance  mode.  By 
contrast,  we  have  favored  transmission  mode,  so  as  to  interrogate  deep 
tumor  regions,  and  we  have  presented  preliminary  studies  in  rat  breast 
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and  prostate  tumors  with  respect  to  various  interventions."^^’^^^’^^^  NIR 
approaches  are  presented  in  detail  in  Chapter  17  of  this  volume. 

Each  technique  has  specific  virtues  and  drawbacks,  which  must  be  con¬ 
sidered  for  any  given  application,  particularly  the  degree  of  invasiveness, 
the  ability  to  generate  maps  of  heterogeneity,  and  the  ability  to  assess  dy¬ 
namic  changes.  In  addition,  the  location  of  a  measurement  (e.g.,  vascular 
versus  tissue  compartments,  the  precision  of  measurements,  and  spatial 
and  temporal  resolution)  must  be  considered.  For  further  details  of  the 
techniques  described  earlier,  the  reader  is  referred  to  the  references.  In 
the  next  section,  we  present  NMR  approaches  in  greater  detail. 

NMR  Approaches  to  Measuring  p02 

Nuclear  magnetic  resonance  (NMR)  is  attractive  because  it  is  inherently 
noninvasive.  Liquid-state  NMR  is  characterized  by  several  parameters,  in¬ 
cluding  signal  amplitude,  chemical  shift  {6),  spin-spin  relaxation  (T2),  and 
spin-lattice  relaxation  (Ti).  Oxygen  could  be  quantified  using  NMR, 
but  this  is  rather  esoteric. Alternatively,  it  has  long  been  recognized  that 
the  oxygen  molecule  (O2)  is  paramagnetic,  causing  increased  spin-lattice  re¬ 
laxation  rates  (Ri  =  1/Ti).  Indeed,  physical  and  theoretical  chemists  must  go 
to  great  lengths  to  rigorously  remove  oxygen  from  solutions  (using  freeze- 
thaw  procedures)  to  achieve  inherent  relaxation  rates  for  studying  nuclear 
interactions. Proton  NMR  studies  have  reported  changes  in  the  water  re¬ 
laxation  rate  as  a  result  of  tissue  oxygenation,^^^  but  many  other  processes 
(metal  ions,  cellularity,  pH,  ionic  strength)  also  cause  relaxation,  and  thus 
it  is  not  suitable  for  detecting  p02,  except  under  rare  circumstances,  such 
as  with  the  eye.^^^  There  is  also  a  substantial  temperature  response,  whereas 
the  relaxivity  due  to  oxygen  is  only  0.0002  s“^/torr.^^^ 

However,  several  investigators  showed  that  the  ^^F  NMR  spin-lattice 
relaxation  rates  for  fluorocarbons  are  much  more  sensitive  to  p02.^^^’^^^ 
Thomas  et  al}^^  pioneered  the  application  of  ^^F  NMR  relaxometry  to 
measure  p02  in  tissues,  in  vivo,  including  lung,  liver,  and  spleen;  several 
other  investigators  demonstrated  feasibility  and  applications,^^"^”^^^  as 


X.  H.  Zhu,  H.  Merkle  et  al,  Magn.  Reson.  Med.  45,  543  (2001). 

M.  A.  Hamza,  G.  Serratice  et  al,  J.  Magn.  Reson.  42,  227  (1981). 

E.  Tadamura,  H.  Hatabu  et  al,  J.  Magn.  Reson.  Imaging  7,  220  (1997). 

A.  Berkowitz,  C.  McDonald  et  al,  Magn.  Reson.  Med.  46,  412  (2001). 

J.-J.  Delpuech,  M.  A.  Hamza  et  al,  J.  Chem.  Rhys.  70,  2680  (1979). 

A.  Hamza,  G.  Serratice  et  al,  J.  Am.  Chem.  Soc.  103,  3733  (1981). 

S.  R.  Thomas,  in  “Magnetic  Resonance  Imaging”  (C.  L.  Partain,  R.  R.  Price,  J.  A.  Patton, 
M.  V.  Kulkarni,  and  A.  E.  J.  James,  eds.),  Vol.  2,  p.  1536.  W.  B.  Saunders,  London,  1988. 
J.  E.  Fishman,  P.  M.  Joseph  et  al,  Magn.  Reson.  Imaging  5,  279  (1987). 


[18] 


DYNAMIC  ASSESSMENT  OF  TUMOR  OXYGENATION 


391 


reviewed  some  years  ago  by  Mason. The  NMR  Ri  of  perfluorocar- 
bons  (PFCs)  varies  linearly  with  and  each  resonance  is  sensitive 

to  p02,  temperature,  and  magnetic  field,  but  importantly,  is  essentially  un¬ 
responsive  to  pH,  CO2,  charged  paramagnetic  ions,  mixing  with  blood,  or 
emulsification.^^^-^^'^ 

A  particular  PFC  molecule  may  have  multiple  resonances,  and  each 
resonance  has  a  characteristic  Ri  response  to  p02.  This  is  attributed  to 
steric  effects  of  O2,  as  it  approaches  the  molecule,^^^  which  implies  that 
perfluorinated  groups,  which  are  both  geometrically  and  magnetically  com¬ 
parable,  should  have  similar  Ri  responses  to  oxygen  tension.  At  a  fixed 
temperature  and  magnetic-field  strength,  the  Ri  response  to  p02  of  any 
single  resonance  obeys  the  simple  formula 

Ri  =  Ria  +  (RipX)  (1) 

where  X  is  the  mole  fraction  of  O2  dissolved  in  the  PFC,  is  the  anoxic 
relaxation  rate,  and  Rjj^  is  the  relaxation  rate  due  to  the  paramagnetic  con¬ 
tribution  of  oxygen.  According  to  Henry’s  law,  the  dissolved  mole  fraction 
is  related  directly  to  the  partial  pressure  of  oxygen, 

PO2  =  KX  (2) 

where  K  represents  Henry’s  constant  for  a  given  solution  of  gas  at  a  speci¬ 
fied  temperature.  By  substitution, 

Ri  =  Ria  +  (Rip/K)p02  (3) 

The  slope  {RiplK)  indicates  the  response  of  a  particular  resonance  to  p02. 

PFCs  essentially  act  as  molecular  amplifiers,  since  the  solubility  of 
oxygen  is  greater  than  in  water,  but  thermodynamics  require  that  the  p02 
in  the  PFC  will  rapidly  equilibrate  with  the  surrounding  medium,  and 
estimates  of  diffusion  suggest  the  equilibration  can  occur  within  seconds. 
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Because  relaxation  is  proportional  to  oxygen  concentration,  the  effect 
will  be  greater  at  a  given  p02  than  for  water.  Importantly,  ions  do  not  enter 
the  hydrophobic  PFC  phase,  and  thus  do  not  affect  the  bulk  relaxation. 
Indeed,  PFCs  are  typically  exceedingly  hydrophobic  and  do  not  mix 
with  the  aqueous  phases,  but  rather  form  droplets  or  emulsions.  Based  on 
these  principles,  PFCs  have  been  applied  to  in  vivo  p02  measurements. 
Characteristics  of  many  diverse  PFCs  are  summarized  in  Table  II. 

At  any  given  magnetic  field  (Bo)  and  temperature  (T),  sensitivity  to 
changes  in  p02  is  given  by  Ri  =  a  +  bp02.  Thus  a  greater  slope  is  impor¬ 
tant,  and  the  ratio  r]  =  b/a  has  been  proposed  as  a  sensitivity  index. Gen¬ 
erally,  a  small  “a”  value  (intercept)  represents  greater  sensitivity,  but  it 
also  generates  longer  Ti  values  under  hypoxic  conditions,  potentially  in¬ 
creasing  data  acquisition  times.  Indeed,  the  Ti  of  hexafiuorobenzene 
(HFB)  at  4.7  T  may  reach  12  s,  potentially  creating  long  imaging  cycles, 
but  this  is  readily  overcome  by  applying  single-shot  (echo  planar)  imaging 
techniques,  as  presented  in  a  later  section. 

Many  PFCs,  such  as  perfluorotributylamine  (PFTB),  perfiubron  (for¬ 
merly  referred  to  as  perfiuorooctyl  bromide;  PFOB),  and  Therox  (F44- 
E),  have  several  ^^F  NMR  resonances,  which  can  be  exploited  to  provide 
additional  information  in  spectroscopic  studies,  but  seriously  hamper  ef¬ 
fective  imaging.  Multiple  resonances  can  lead  to  chemical  shift  artifacts 
in  images,  which  compromise  the  integrity  of  relaxation  time  measure¬ 
ments,  though  they  can  be  avoided  by  selective  excitation,  or  detection, 
chemical  shift  imaging,  deconvolution,  or  sophisticated  tricks  of  NMR 
spin  physics.^^^’^^^’^^^"^^^  These  approaches  add  to  experimental  com¬ 
plexity  and  are  generally  associated  with  lost  signal  to  noise  ratio  (SNR). 
Thus  we  strongly  favor  PFCs  with  a  single  resonance,  and  we  will  describe 
the  use  of  though  some  research  groups  favor 

5-crown-5-ether 
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TABLE  II 

NMR  Characteristics  and  Applications  of  PFCs  for  Tissue  Oximetry 


PFC 

Sensitivity 
to  PO2'* 

Temp. 

sensitivity 

(torr/°) 

Magnetic 
field  Bo(T) 

Application/ 

comments 

References 

Hexafluorobenzene 

(HFB) 

A  =  0.0835 

B  =  0.001876 

0.13 

4.7 

Rat  breast  tumor, 
prostate  tumor, 
human  lymphoma 
xenograft 

Hunjan,  Zhao  et 

Zhao,  Constantinescu  et 

Mason,  Ran  et  al}^\ 

Song,  Constantinescu  et  alP^; 

Zhao,  Ran  et  alP^ 

HFB 

A  =  0.074 

B  =  0.00158 

4.7 

Rat  prostate  tumor 

Hunjan,  Mason  et 

Mason,  Constantinescu  et  al}^^ 

HFB 

A  =  0.093 

B  =  0.103 

1.40 

7 

Phantom 

Mason,  Rodbumrung  et  al}^^ 

Perfluoro  - 1 5  -  Crown- 
5-ether  (15C5) 

A  =  0.345 

B  =  0.0034 

2.94 

2.0 

Tumor  cells 

Helmer,  Han  et 

15C5 

A  =  0.333 

B  =  0.0033 

3.98 

2.0 

Mouse  tumor, 
spleen,  liver 

Dardzinski  and  Sotak^^^ 

15C5 

A  =  0.44 

B  =  0.0028 

4.3 

Human  glioma 
tumor  in  mice 

van  der  Sanden,  Heerschap  et 
van  der  Sanden,  Heerschap  et  alP^ 

15C5 

A  =  0.375 

B  =  0.00198 

4.7 

Rat  breast  tumor 

Fan,  River  et  alP 

15C5 

A  =  0.362 

B  =  0.1239 

4.7 

Rat  brain 

Duong,  Ladecola  et 

Perfluorotributyl-amine 
(FC-43)  PFTB 

PFTB 

A  =  1.09 

B  =  0.00623 

4.43 

0.14 

0.14 

Pig  liver,  spleen, 
lung 

Rat  liver,  spleen, 
lung 

Thomas,  Pratt  et  al}^^ 

Pratt,  Zheng  et 

Thomas,  Gradon  et  al}^^ 

(continued) 
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PFC 

Sensitivity 
to  p02^ 

Temp. 

sensitivity 

(torr/°) 

Magnetic 
field  Bo(T) 

Application/ 

comments 

References 

PFTB 

A 

=  0.684 

4.7 

Rabbit  eye 

Berkowitz,  Wilson  et 

B 

=  0.00305 

Wilson,  Berkowitz  et  al.^ 

PFTB 

A 

=  0.9072 

1.5 

Human  eye 

Wilson,  Berkowitz  et  al}^^ 

B 

=  0.004486 

PFTB 

A 

=  0.8848 

8.17 

7 

Mouse  Meth-A  tumor 

Mason,  Shukla  et  al}^^ 

B 

=  0.1307 

and  heart 

Perfluorotripropylamine 

A 

=  0.314 

1.9 

Rat  subcutaneous 

Fishman,  Joseph  et  al.^ 

(FTPA) 

B 

=  0.002760"^ 

tumor 

FTPA 

A 

=  0.301 

1.4 

Rat  spleen,  lung, 

Fishman,  Joseph  et  al}^^ 

B 

=  0.00312 

liver 

FTPA 

A 

=  0.4052 

2.0 

Rat  liver,  spleen 

Holland,  Kennan  et  al}^^ 

B 

=  0.0023 

FTPA 

4.4 

Cells 

Taylor  and  Deutsch"^ 

Bis-perfluoro-butylethylene 

A 

=  0.3421 

7.05 

Rat,  alginate  capsules 

Noth,  Grohn  et  al}^^ 

(F-44E) 

B 

=  0.11172 

F-44E 

A 

=  0.342 

7.05 

Rat  spleen,  liver. 

Noth,  Morrissey  et  al}^^ 

B 

=  0.1201 

abdominal  aorta. 

vena  cava 

F-44E 

A 

=  0.2525 

0.59 

2.0 

Mouse  tumor 

Hees  and  Sotak^^^ 

B 

=  0.16527 

Perfluoroctyl-bromide 

A 

=  0.517 

9.4 

Rat  heart 

Shukla,  Mason  et  alf 

(perflubron)  (PFOB) 

B 

=  0.0038 

PFOB 

A 

=  0.2677 

1.26 

4.7 

Prostate  tumor  in  rat 

Antich  et  al}^^ 

B 

=  0.12259 
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PFOB 

A 

=  0.328 

2.85 

7 

Phantom 

Mason,  Shukla  et  al.^^^ 

B 

=  0.12137 

PFOB 

2.0 

Rat  tumor 

Sostman,  Rockwell  et  al.^^ 

PFOB 

A 

=  0.085 

1.5 

Rabbit  liver 

Tran,  Guo  et  al.^^^ 

B 

=  0.0033 

PFOB 

8.5 

Pig  liver,  lung,  spleen 

Millard  and  McGoron^ 

PFOB 

1.45 

Rat  lung.  Mouse  lung 

Thomas,  Clark,  Jr.  et  al.^ 

PFOB 

1.5 

Pig  lung 

Laukemper-Ostendorf, 

Scholz  et  al.^^^ 

Perfluoro-2,2,2',2'- 

A 

=  0.50104 

2.0 

Phantom 

Sotak,  Hees  et  al.^ 

tetr  amethy  1-4,4' - 

bis(l,3-fioxolane) 

(PTBD) 

B 

=  0.1672 

Some  original  papers  presented  calibration  curves  in  other  forms  (e.g.,  including  coefficients  for  temperature  dependence).  In  those  cases, 
equations  have  been  derived  assuming  37°.  Where  a  PFC  has  more  than  one  resonance,  the  equation  presented  is  either  for  the  most  sensitive 
signal,  or  the  equation  used  where  the  signal  may  have  been  unresolved. 

(s“^)  =  A  +  B  X  p02  (torr). 

^C.  A.  Wilson,  B.  A.  Berkowitz  et  ah,  Exp.  Eye  Res.  55,  119  (1992). 

""  J.  E.  Fishman,  P.  M.  Joseph  et  al..  Invest.  Radiol.  24,  65  (1989). 

^J.  Taylor  and  C.  J.  Deutsch,  Biophys.  J.  53,  227  (1988). 

^H.  P.  Shukla,  R.  P.  Mason  et  al.,  Magn.  Reson.  Med.  35,  827  (1996). 

^R.  W.  Millard  and  A.  J.  McGoron,  Artif.  Cells  Blood  Subst.  Immobil.  Biotechnol.  22,  1251  (1994). 

^S.  R.  Thomas,  L.  C.  Clark,  Jr.  et  al,  J.  Comput.  Assist.  Tomogr.  10,  1  (1986). 

^C.  H.  Sotak,  P.  S.  Hees  et  al.,  Magn.  Reson.  Med.  29,  188  (1993). 
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Ri  is  sensitive  to  temperature,  although  the  response  varies  greatly  be¬ 
tween  PFCs  and  between  individual  resonances  of  each  individual  PFC. 
Over  small  temperature  ranges,  a  linear  correction  to  calibration  curves 
is  appropriate,  but  over  larger  temperature  ranges,  the  response  can  be 
complex,  as  investigated  extensively  by  Shukla  et  al}^^  for  several  PFCs. 
Differential  sensitivity  of  pairs  of  resonances  to  p02  and  temperature 
allowed  Mason  et  al}^^  to  simultaneously  determine  both  parameters  by 
solving  simultaneous  equations.  However,  generally  it  is  preferable  for  a 
p02  sensor  to  exhibit  minimal  response  to  temperature,  since  this  is 
not  always  known  precisely  in  vivo  and  temperature  gradients  may  occur 
across  tumors.  As  shown  in  Table  II,  even  a  relatively  small  error  in  tem¬ 
perature  estimate  can  introduce  a  sizable  discrepancy  into  the  apparent 
PO2;  for  example,  the  relative  error  introduced  into  a  p02  determination 
by  a  1°  error  in  temperature  estimate  ranges  from  8  torr/°  for  FFTB^"^^ 
to  3  torr/°  for  PFOB  (perfiubron)^'^^  or  15C5^^^  and  0.1  torr/°  for  HFB,^"^^ 
when  PO2  is  actually  5  torr.  It  must  be  noted  that  error  depends  on  actual 
PO2  and  the  error  varies  with  magnetic  field  and  temperature.  Ri  response 
does  depend  on  magnetic  field,  necessitating  calibration  curves  for  each 
type  of  magnet  system  (e.g.,  1.5,  4.7,  or  7  T).  Thus  comparison  of  PFC 
utility  for  p02  measurements  is  complicated  by  the  field  used  for  specific 
published  investigations,  and  in  Table  II,  we  consider  sensitivity  as 
presented. 

Choice  of  PFC  may  be  governed  by  practical  considerations,  such  as 
cost  and  availability,  since  several  products,  particularly  proprietary  emul¬ 
sions,  may  be  difficult  to  obtain.  HFB  and  15C5  offer  the  immediate  advan¬ 
tage  of  a  high  symmetry  and  a  single  ^^F  NMR  resonance.  This  offers 
maximum  SNR  and  simplifies  imaging,  which  may  otherwise  require  fre¬ 
quency  selective  excitation,  deconvolution,  or  other  NMR  tricks  to  avoid 
chemical  shift  artifacts. 

Route  of  Administration.  The  most  popular  route  for  the  delivery  of 
PFCs  is  as  emulsions  injected  intravenously.  Given  the  extremely  hydro- 
phobic  nature  of  PFCs,  they  do  not  dissolve  in  blood  directly,  but  may  be 
formulated  as  biocompatible  emulsions.  Much  effort  has  been  applied 
to  formulate  stable  homogenous  emulsions,  as  reviewed  else  where. 
Following  IV  infusion,  the  emulsion  circulates  in  the  vasculature  with  a 


P.  Shukla,  R.  P.  Mason  et  at,  J.  Magn.  Reson.  B  106,  131  (1995). 

R.  P.  Mason,  H.  P.  Shukla  et  al,  Magn.  Reson.  Med.  29,  296  (1993). 

R.  P.  Mason,  H.  P.  Shukla  et  al,  Biomater.  Artif.  Cells  Immobilization  Biotechnol  20,  929 
(1992). 

R.  P.  Mason,  W.  Rodbumrung  et  al,  NMR  Biomed.  9,  125  (1996). 

J.  G.  Riess,  Biomater.  Artif.  Cells  Immobilization  Biotechnol  20,  183  (1992). 
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typical  half-life  of  12  h,  (depending  on  the  nature  of  the  emulsion)  provid¬ 
ing  substantial  clearance  within  2  days.^"^^’^"^^  Primary  clearance  is  by 
macrophage  activity,  leading  to  extensive  accumulation  in  the  liver,  spleen, 
and  bone  marrow.^^^’^"^^  Indeed,  this  is  a  major  shortcoming  of  IV  delivery, 
since  animals  may  exhibit  extensive  hepatomegaly  or  splenomegaly. The 
emulsions  are  not  toxic,  and  other  than  causing  swelling,  appear  not  to 
cause  health  problems.  PFC  clearance  occurs  from  the  liver  with  a  typical 
half-life  of  60  days  for  perfluorotripropylamine  and  3  days  for  perflubron, 
with  primary  clearance  by  migration  to  the  lungs  and  exhalation. 

Some  investigators  have  examined  p02  of  tissues,  while  PFC  remained 
in  the  blood,  providing  a  vascular  Flow  can  generate  artifacts, 

and  correction  algorithms  have  been  proposed. Many  investigators  have 
measured  p02  in  liver,  spleen,  and  tumors  following  clearance  from  the 
blood,  thus  providing  measurement  of  tissue 

Both  spectroscopic  and  imaging  approaches  have  been  applied  to  tissue 
p02  measurements  depending  on  the  available  SNR.  It  appears  that  uptake 
and  distribution  efficiency  vary  with  tumor  type,  but  in  general,  maximum 
signal  is  detected  from  the  tumor  periphery  corresponding  with  regions  of 
greater  perfusion.^^^’^^^’^^^’^^"^’^^^’^^^  Several  reports  have  examined  changes 
in  tumor  p02  in  response  to  acute  interventions  such  as  vasoactive  drugs 
and  hyperoxic  Spectroscopic  time  resolution  has 

ranged  from  seconds  to  minutes,^^^’^^^  whereas  imaging  often  takes  longer.^^^ 
Long  tissue  retention  facilitates  chronic  studies  during  tumor  develop¬ 
ment,  and  progressive  tumor  hypoxiation  has  been  observed  over  extended 
time  periods  of  many  days.^^"^’^^^  Correlated  ^^F  and  proton  MRI  suggest 


M.  C.  Malet-Martino,  D.  Betbeder  et  al,  J.  Pharm.  Pharmacol  36,  556  (1984). 

T.  F.  Zuck  and  J.  G.  Riess,  Crit.  Rev.  Clin.  Lab.  ScL  31,  295  (1994). 

R.  P.  Mason,  P.  P.  Antich  et  al,  Magn.  Reson.  Imaging  7,  475  (1989). 

W.  I.  Rosenblum,  M.  G.  Hadfield  et  al.  Arch.  Pathol  Lab.  Med.  100,  213  (1976). 

R.  F.  Mattrey  and  D.  C.  Long,  Invest.  Radiol  23,  s298  (1988). 

Eidelberg,  G.  Johnson  et  al,  J.  Cereb.  Blood  Flow  P.  Metab.  8,  276  (1988). 

T.  Higuchi,  S.  Naruse  et  al,  in  “7th  SMRM,”  p.  435,  1988. 

R.  P.  Mason,  R.  L.  Nunnally  et  al,  Magn.  Reson.  Med.  18,  71  (1991). 

P.  S.  Hees  and  C.  H.  Sotak,  Magn.  Reson.  Med.  29,  303  and  erratum  29,  716  (1993). 

R.  P.  Mason,  P.  P.  Antich  et  al,  Int.  J.  Radiat.  Oncol  Biol  Phys.  29,  95  (1994). 

S.  R.  Thomas,  R.  W.  Millard  et  al,  Artif.  Cells  Blood  Subst.  Immobil  Biotechnol  22, 1029 
(1994). 

N.  J.  Baldwin  and  T.  C.  Ng,  Magn.  Reson.  Imaging  14,  541  (1996). 

K.  G.  Helmer,  S.  Han  et  al,  NMR  Biomed.  11,  120  (1998). 

P.  P.  Antich,  R.  P.  Mason,  A.  Constantinescu  et  al,  Proc.  Soc.  Nucl  Med.  35,  216P  (1994). 
D.  J.  O.  McIntyre,  C.  L.  McCoy  et  al,  Curr.  Scl  76,  753  (1999). 
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that  PFC  does  not  redistribute,  but  remains  associated  with  specific  tissues, 
analogous  to  tree  rings. Thus  in  principle,  a  whole  tumor  can  be  in¬ 
vestigated  by  administering  successive  doses  of  PFC  emulsion  during 
growth. 

PFC  emulsions  may  also  be  administered  intraperitoneally  (IP), 
resulting  in  similar  distribution  to  IV  administration  (unpublished  observa¬ 
tions).  Given  the  volatile  nature  of  many  PFCs,  they  could  be  inhaled,  but 
although  this  is  a  popular  route  for  delivery  of  anesthetics  and  blood  flow 
tracers,  it  does  not  appear  to  have  been  widely  exploited  for  oximetry. 
Nonetheless,  aerosols  have  been  delivered  to  the  lungs  by  inhalation  to 
facilitate  p02  measurements. 

Two  approaches  have  been  applied  to  circumvent  reticuloendothelial 
uptake.  PFC  has  been  incorporated  in  polyalginate  beads  for  direct  im¬ 
plantation  at  a  site  of  interest.^^"^’^^^  We  favor  direct  IT  injection  of  neat 
PFC,  allowing  any  region  of  interest  in  a  tumor  to  be  interrogated  immedi¬ 
ately.  Use  of  a  fine  needle  ensures  minimal  tissue  damage,  as  described  in 
detail  in  a  later  section.  Others  have  used  direct  injection  of  emulsions  into 
tumors,  but  this  increases  the  volume  considerably,  making  it  more  inva- 
sive.^^  Investigators  have  suggested  that  emulsification  improves  retention 
at  the  site  of  injection.  Direct  injection  of  neat  PFC  also  has  been  used  to 
investigate  retinal  oxygenation^^^“^^^  and  cerebral  oxygenation  in  the 
interstitial  and  ventricular  spaces. 

As  described  in  the  following  section  we  favor  direct  intratumoral  injec¬ 
tion  of  neat  HFB  followed  by  echo  planar  imaging  to  generate  p02  maps  in 
tumors. 


FREDOM  (Fluorocarbon  Relaxometry  using  Echo  Planar  Imaging  for 
Dynamic  Oxygen  Mapping) 

Recognizing  that  tumors  are  heterogeneous  and  that  p02  may  fluctuate, 
we  developed  a  procedure,  which  allows  repeated  quantitative  maps  of 
regional  p02  to  be  achieved  with  multiple  individual  locations  simulta¬ 
neously  in  6.5  min  with  a  precision  of  1-3  torr,  when  p02  is  in  the 


S.  R.  Thomas,  L.  Gradon  et  al,  Invest.  Radiol.  32,  29  (1997). 

Noth,  P.  Grohn  et  al.,  Magn.  Reson.  Med.  42,  1039  (1999). 

U.  Zimmermann,  U.  Noth  et  al.,  Artif.  Cells  Blood  Subst.  Immobil.  Biotechnol.  28,  129 
(2000). 

A.  Berkowitz,  C.  A.  Wilson  et  al..  Invest.  Ophthalmol.  Vis.  Sci.  32,  2382  (1991). 

C.  Wilson,  B.  Berkowitz  et  al..  Arch.  Ophthalmol.  110,  1098  (1992). 

W.  Zhang,  Y.  Ito  et  al..  Invest.  Ophthalmol.  Visual  Sci.  44,  3119  (2003). 
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Fig.  1.  Hexafluorobenzene  (HFB)  is  a  perfluorocarbon  (PFC)  exhibiting  extensive 
symmetry. 


range  0-10  torr.^^^  We  have  applied  FREDOM  to  diverse  tumor  types  and 
interventions,  as  reviewed  in  a  later  section. 

MRI  is  attractive  because  it  is  readily  available  at  many  institutions.  For 
small  animal  work,  ^^F  NMR  is  widely  available  at  4.7,  7,  and  9.4  T  by 
minor  adaptation  of  routine  instrumentation,  [e.g.,  retuning  proton  radio¬ 
frequency  (RF)  coils].  Within  the  recent  past  ^^F  MRI  is  also  becoming 
available  on  clinical  systems,  facilitating  translation  of  these  techniques  to 
patients.  ^^F  NMR  is  particularly  facile  because  there  is  essentially  no 
background  signal  in  tissues  to  interfere  with  measurements,  yet  the  reso¬ 
nance  frequency  and  sensitivity  approach  that  of  proton  NMR.  The  pio¬ 
neering  work  of  Thomas^^^  showed  that  tissue  p02  could  be  imaged  in 
various  organs  based  on  the  ^^F  NMR  spin-lattice  relaxation  rate  (Ri)  of 
PFC  reporter  molecules  following  IV  infusion.  Prompted  by  these  studies, 
we  surveyed  a  number  of  PFCs  and  identified  that  HFB  (Fig.  1)  has  many 
virtues  as  a  p02  reporter. Symmetry  provides  a  single  narrow  ^^F  NMR 
signal,  and  the  spin-lattice  relaxation  rate  is  highly  sensitive  to  changes 
in  p02,  yet  minimally  responsive  to  temperature  (Fig.  2)}^^  HFB  also  has 
a  long  spin-spin  relaxation  time  (T2),  which  is  particularly  important 
for  imaging  investigations.  From  a  practical  perspective,  HFB  is 
cheap  (<$2/g)  and  readily  available  commercially  in  high  purity  (>99%). 
We  obtain  supplies  from  Lancaster  Synthesis  (Windham,  NH),  though 
many  other  fine  chemical  supply  houses  also  offer  HFB.  We  do  favor 
bottles  over  sealed  ampoules,  since  they  are  easier  to  handle.  HFB  is  well 
characterized  in  terms  of  lack  of  toxicity, exhibiting  no  mutage- 
nicity,^^^  teratogenicity,  or  fetotoxicity,^^^  and  the  manufacturer’s  material 


S.  Gorsman  and  T.  A.  Kapitonenko,  Izv.  Estestvennonauchu.  Inst.  Pevinsk.  15,  155 
(1973). 

I.  M.  C.  M.  Rietjens,  A.  Steensma  et  al.,  Eur.  J.  Pharmacol.  293,  292  (1995). 

K.  E.  Mortelmans  and  V.  F.  Simmon,  Gov.  Rep.  Announce.  Index  (US)  81,  2555  (1981). 
K.  D.  Courtney  and  J.  E.  Andrews,  /.  Environ.  Sci.  Health  B  19,  83  (1984). 
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A  B 


p02  (torr)  Temperature  (°) 

Fig.  2.  (A)  The  NMR  spin-lattice  relaxation  rate  (Ri)  of  HFB  shows  a  linear  response 
to  PO2:  at  37°  and  4.7  T,  Ri  =  0.0835  +  0.001876  x  p02.  (B)  The  ^‘^F  NMR  Ri  of  HFB  shows 
minimal  response  to  temperature  in  the  physiologic  range;  separate  curves  are  shown  at  four 
different  p02s:  0,  1,  9.8,  and  21%  O2. 


data  safety  sheet  indicates  LD50  >  25  g/kg  (oral — rat)  and  LC50  95  g/ni^/2  h 
(inhalation — mouse). 

HFB  had  been  proposed  as  a  veterinary  anesthetic  and  has  been  used  in 
many  species,  including  ponies,  sheep,  cats,  dogs,  rats,  and  mice,  but  was 
abandoned  because  of  its  high  volatility  (boiling  point  81°)  and  low  flash 
point  (10°).^^^  This  presented  unacceptable  dangers  in  veterinary  suites 
for  inhalation  anesthetics.  It  is  not  a  problem  in  our  studies,  where  small 
quantities  of  liquid  (typically,  50  /jI)  are  injected  directly  into  the  tumor. 
HFB  requires  no  special  storage,  other  than  a  sealed  bottle  to  prevent 
evaporation.  Melting  point  is  4-6°,  and  density  is  1.62  g/ml. 

Methodology 

Tumor  Preparation.  We  have  applied  this  technique  to  tumors  in  rats 
and  mice,  but  our  methodology  will  focus  on  our  standard  application  to 
rats.  As  for  electrode  techniques,  the  tumor  must  be  accessible  (e.g.,  sub¬ 
cutaneous  in  flank  of  thigh),  though  we  favor  the  pedicle  model,  which  pro¬ 
vides  a  tumor  remote  from  the  body,  analogous  to  an  additional  limb.^^"^ 
This  model  is  optimal  for  selective  therapy  such  as  local  hyperthermia,  ir¬ 
radiation,  or  excision.  Noting  the  increasing  interest  in  orthotopic  tumors, 
the  approach  is  also  facile  in  breast  tumors  in  the  mammary  fat  pad  and 
applicable  in  the  prostate  with  a  little  practice. 

L.  W.  Hall,  S.  R.  K.  Jackson  et  al,  in  “Recent  Progress  in  Anaesthesiology  and 
Resuscitation”  (A.  Arias,  R.  Llaurado,  M.  A.  Nalda,  and  J.  N.  Lunn,  eds.),  p.  201.  Excerpta 
Medica,  Oxford,  UK,  1975. 

W.  Hahn,  P.  Peschke  et  al,  Magn.  Reson.  Imaging  11,  1007  (1993). 
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Rats  are  preanesthetized  with  ketamine  hydrochloride  (100  mg/ml)  as 
a  relaxant  and  maintained  under  general  gaseous  anesthesia  with  air 
(1  dm^/min)  and  1.2%  isoflurane  [we  do  note  that  the  appropriate  concen¬ 
tration  of  anesthesia  may  depend  on  strain;  for  example,  Copenhagen  rats 
require  (and  tolerate)  higher  concentrations  than  Fisher  rats].  In  our  latest 
refinement  of  the  procedure,  isoflurane  may  be  stopped  for  a  period  of 
minutes  during  HFB  administration,  since  rats  occasionally  exhibit  respira¬ 
tory  distress,  which  may  be  caused  by  the  anesthetic  properties  of  HFB 
interacting  with  isoflurane.  Within  a  few  minutes,  rats  are  stable  and  rou¬ 
tine  isoflurane  anesthesia  is  maintained.  In  our  earlier  work  we  used 
0.5%  methoxyflurane  in  33%  oxygen  with  66%  but  isoflurane 

appears  to  be  a  less  stressful  anesthetic. 

HFB  is  deoxygenated  by  bubbling  nitrogen  for  5  min  before  use.  We 
have  previously  shown  in  hypoxic  tumor  biopsies  that  use  of  aerated 
HFB  may  introduce  a  systematic  apparent  elevation  in  p02  (~2-3  torr). 
HFB  is  injected  directly  into  tumors  using  a  gas  tight  syringe  with  a 
custom-made  fine  sharp  needle  (32G  #7803-04;  Hamilton,  Reno,  NV). 
Generally,  HFB  is  administered  along  three  to  five  tracks  in  the  form  of 
a  fan  in  a  single  central  plane  of  the  tumor  coronal  to  the  rat’s  body 
(Fig.  3A).  The  needle  is  inserted  manually  to  penetrate  across  the  whole 
tumor  and  withdrawn  mm  to  reduce  pressure,  and  3  //I  HFB  is 
deposited.  The  needle  is  repeatedly  withdrawn  a  further  2-3  mm,  and 
additional  HFB  is  deposited.  Typically,  HFB  is  deliberately  deposited  at 
about  16  individual  locations  per  tumor,  in  both  the  central  and  pe¬ 
ripheral  regions  of  the  tumors,  to  ensure  that  the  interrogated  regions  are 
representative  of  the  whole  tumor. 

The  animal  is  placed  on  its  side  in  a  cradle  with  a  thermal  blanket  to 
maintain  body  temperature.  A  fiber-optic  probe  is  inserted  rectally  to 
monitor  core  temperature.  Temperature  measurement  is  optional,  since 
the  Ri  is  essentially  invariant  with  temperature  and  does  not  need  correc¬ 
tion  for  p02  estimates.  Of  course,  temperature  regulation  and  measure¬ 
ment  is  important  to  ensure  stable  tumor  physiology.  An  MR-compatible 
pulse  oximeter  equipped  with  “rat”  software  (8600V;  Nonin,  Inc.,  Ply¬ 
mouth,  MN)  may  be  applied  to  a  hind  foot  to  monitor  arterial  oxygenation 
(Sa02)  and  heart  rate  (optional,  but  provides  additional  useful  data 
regarding  animal  health  and  physiology). 

Most  of  our  MR  experiments  were  performed  using  an  Omega  CSI  4.7 
horizontal  bore  magnet  system  with  actively  shielded  gradients  (GE 
systems,  acquired  by  Bruker  Instrument,  Inc.,  Fremont,  CA).  Recently, 
our  MR  system  has  been  upgraded  to  a  Varian  Unity  INOVA,  providing 

R.  P.  Mason,  A.  Constantinescu  et  al.,  Radial  Res.  152,  239  (1999). 
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Fig.  3.  (A)  Recommended  pattern  of  IT  injection.  Several  tracks  in  the  form  of  a  fan  in  a 
single  plane.  (B)  Four  contiguous  MR  images  (5  mm  thick)  of  a  representative  Dunning 
prostate  R3327-AT1  tumor  (volume  =  1.6  cm^).  Top:  MRI;  middle:  corresponding  ^^F; 

bottom:  overlay  to  show  interrogated  regions.  B  marks  point  of  attachment  of  the  tumor 
to  back  of  rat. 
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MBEST 


RF  pulses 
G,(Slice) 


Fig.  4.  MRI  sequence  used  to  achieve  relaxation  curves. 


enhanced  data  acquisition  and  processing  software  and  stronger  imaging 
gradients.  A  tunable  (^H/^^F)  MR  coil,  2  or  3  cm  in  diameter  matched  to 
the  tumor  size  (constructed  from  a  cylindrical  copper  tube  about  2  cm  deep 
and  acting  as  a  single-turn  solenoid),  is  placed  around  the  tumor-bearing 
pedicle.  Shimming  is  performed  on  the  proton  water  signal  to  a  typical  line- 
width  of  60-100  Hz.  Proton  images  are  obtained  for  anatomic  reference 
using  a  three-dimensional  (3D)  spin-echo  sequence.  The  coil  is  then  re¬ 
tuned  in  place  to  188.27  MHz,  and  corresponding  MR  images  are 
obtained.  Overlaying  the  ^^F  MR  images  on  the  corresponding  proton 
images  reveals  the  distribution  of  HFB  (Fig.  3B).  Typically,  ^^F  NMR 
signal  is  obtained  from  6-10%  of  the  total  tumor  voxels. 

p02  maps  are  obtained  using  our  standard  MR  oximetry  proto- 
This  applies  pulse  burst  saturation  recovery  (PBSR)  echo 
planar  imaging  (EPI)  relaxometry  using  the  ARDVARC  (alternated  relax¬ 
ation  (ielays  with  variable  acquisitions  to  reduce  clearance  effects)  acquisi¬ 
tion  protocol  to  map  the  tumors  (Fig.  4).i'75,i77  ^pj  ^  single  spin-echo 
with  “Flipped”  phase  encoding  (modulus-blipped  echo-planar  single-pulse 
technique  [MBEST]),  although  other  EPI  sequences  should  be  equally  ap¬ 
plicable.  We  chose  the  PBSR  approach  to  Ti  relaxation  measurements  for 
historical  reasons;  our  earliest  work  used  ^^F  NMR  spectroscopy  of  tumors 
with  excitation  and  detection  based  on  surface  coils.^^^  The  PBSR  ap¬ 
proach  is  suitable  for  use  with  the  nonuniform  excitation  typical  of  surface 

B.  R.  Barker,  R.  P.  Mason  et  al,  Magn.  Reson.  Imaging  11,  1165  (1993). 

D.  Le,  R.  P.  Mason  et  al,  Magn.  Reson.  Imaging  15,  971  (1997). 
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Fig.  5.  Typical  relaxation  curve  showing  order  of  data  acquisition,  y  =  A  [1  —  (1  +  W) 
exp(-Ri  X  r)]  and  Ti  =  10.34  ±  0.16  s,  A  =  99.5,  W  =  0.04,  and  hence  p02  =  7.0  ±  0.7  torr. 


coils  as  presented  by  Evelhoch  and  Ackerman. For  EPI,  we  do  apply 
carefully  calibrated  7r/2  pulses  to  ensure  accurate  refocusing  of  images, 
but  we  continue  to  favor  the  PBSR  approach,  since  it  provides  particularly 
rapid  estimation  of  Ti  by  avoiding  the  need  for  extended  relaxation  recov¬ 
ery  times  between  acquisitions.  Saturation  is  achieved  by  a  series  of  20  non- 
spatially  selective  7r/2  pulses  with  50-ms  spacing,  which  is  sometimes  called 
a  “Comb”  format. 

Typical  FREDOM  parameters  use  32  x  32  data  points  across  a  field  of 
view  of  40  X  40  mm,  providing  1.25  mm  in-plane  resolution.  Recently,  we 
have  applied  a  slice  selection  gradient,  providing  5-mm-thick  slices,  and 
hence,  8-/il  voxels.  However,  the  deposition  of  HFB  is  designed  to  occur 
in  a  plane,  and  thus  the  discrete  distribution  of  the  reporter  molecule  itself 
can  be  used  to  define  the  slice.  We  apply  a  PBSR  preparation  sequence  be¬ 
cause  it  is  ideally  suited  for  measuring  diverse,  long  Ti  values.  Unlike  more 
traditional  inversion  recovery  sequences,  there  is  no  need  to  wait  >5  Ti  be¬ 
tween  successive  images.  As  shown  in  Fig.  5,  we  use  14  delays  in  the  order 
90  s,  200  ms,  60  s,  400  ms,  40  s,  600  ms,  20  s,  800  ms,  16  s,  1  s,  8  s,  1.5  s,  4  s, 
and  2  s,  selected  to  cover  the  whole  range  of  TiS  (viz  p02  values).  Many 
papers  have  been  published  on  optimizing  relaxation  curves  and  choosing 
parameters  to  enhance  precision. Our  experience  suggests  that  our 

J.  L.  Evelhoch  and  J.  J.  H.  Ackerman,  /.  Magn.  Resort.  53,  52  (1983). 

A.  P.  Crawley  and  R.  M.  Henkelman,  Magn.  Reson.  Med.  7,  23  (1988). 

S.  J.  Doran,  J.  J.  Attard  et  al,  J.  Magn.  Reson.  100,  101  (1992). 

F.  Franconi,  F.  Seguin  et  al,  Medical  Phys.  22,  1763  (1995). 
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parameters  are  appropriate  for  interrogating  a  broad  range  of  p02  values, 
though  they  have  not  been  rigorously  or  theoretically  optimized.  More 
delays  could  improve  curve  fitting,  but  would  increase  experiment  time. 
We  believe  fewer  delays  would  degrade  quality  of  the  curve  fitting. 

Traditional  Ti  measurement  sequences  acquire  data  with  delays  in 
monotonic  order,  whereas  we  alternate  longer  and  shorter  delays  to  mini¬ 
mize  any  systematic  errors,  which  would  be  introduced  if  the  signal  ampli¬ 
tude  varies  during  the  measurement.  We  have  found  that  HFB  clears  from 
tissue  with  a  typical  half-life  of  600  min,^^^  which  would  introduce  errors 
into  the  amplitude.  Although  the  total  acquisition  time  for  a  Ti  map  is 
6.5  min,  we  reduce  the  time  between  first  and  last  acquisitions  further 
by  applying  the  longest  delay  first.  Our  experience  shows  that  it  is  import¬ 
ant  to  measure  at  least  two  data  points  on  the  plateau  of  the  curve  (i.e., 
>5  X  Ti,  and  we  choose  60  and  90  s,  respectively.  The  data  points  with 
longer  recovery  times  have  greater  SNR,  and  we  find  that  a  minimum 
SNR  =  10  is  required  to  produce  satisfactory  Ti  curves  and  p02  estimates. 
The  variation  in  amplitude  between  longest  and  shortest  delays  approaches 
100-fold.  The  poorest  SNR  data  points  (short  delay;  r)  could  compromise 
the  quality  of  the  Ti  curves,  and  thus  we  obtain  multiple  acquisitions  at 
these  times  to  provide  enhanced  SNR  by  signal  averaging.  Because  the  r 
values  are  short,  it  adds  little  to  the  overall  experimental  time.  We  use 
NA  =  12  for  r  <  1  s;  NA  =  8  for  r  =  1  s;  NA  =  4  for  r  =  1.5,  2,  and  4  s; 
NA  =  2  for  r  =  8  s;  and  16  s  and  NA  =  1  for  the  longer  delays.  Signal  am¬ 
plitudes  are  corrected  for  the  additional  acquisitions,  and  a  2-d  Fourier 
transform  (FT)  is  applied  to  each  image.  Curve  fitting  is  then  applied  on 
a  voxel-by-voxel  basis  to  each  image  set.  Data  quality  could  probably  be 
enhanced  by  application  of  apodization  and  filtering  functions,  such  as 
Fermi  and  Hanning  filters,  though  we  have  not  yet  implemented  these 
approaches.  The  acquisition  protocol  is  not  part  of  the  standard  software 
supplied  by  the  manufacturer,  but  we  can  assist  interested  investigators. 

Data  are  transferred  to  a  personal  computer  (PC)  for  further  analysis 
using  a  program  created  in  our  laboratory  (Timap  [Fig.  6],  PASCAL  by 
Dennis  Le).  Timap  recognizes  issues  in  data  analysis  and  provides  filters, 
clustering  algorithms,  and  temporospatial  correlation  to  assist  in  effective 
data  reduction.  At  top  right  of  Fig.  6  is  a  map,  which  may  be  toggled  to  dis¬ 
play  signal  amplitude,  Ti,  Ti  error  (Ti  err),  or  p02  (Fig.  7).  Blank  pixels  in¬ 
dicate  that  no  curve  fit  could  be  achieved.  Gray  pixels  indicate  a  successful 
curve  fit  and  potential  p02  value,  but  with  errors  beyond  the  specified 
range.  In  this  example,  thresholds  were  set  as  Ti  err  <  2.5  s  and  the  ratio 
Ti  err/Ti  <  50%.  Tightening  these  thresholds  can  produce  higher-quality 
data,  but  eliminates  more  data.  Colored  pixels  report  measurements 
satisfying  the  threshold  criteria,  and  these  are  tabulated  in  the  form  of  a 
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Fig.  6.  Example  Timap  used  for  data  filtering,  clustering,  and  analysis. 
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histogram  (center),  which  can  be  plotted  with  selected  data  bins  and  ranges 
of  interest.  At  bottom  right  of  Fig.  6  are  data  for  each  pixel,  including  co¬ 
ordinate  and  fitted  parameters.  In  the  box  on  the  left  are  shown  the  number 
of  fitted  data  points,  number  of  accepted  data  points,  and  statistics  such  as 
mean  Ti,  p02,  and  corresponding  amplitude-weighted  values. 

p02  is  calculated  on  the  basis  of  a  curve  fit  to  exponential  data.  Any 
curve  fitting  is  associated  with  uncertainty,  and  indeed  our  procedure  deter¬ 
mines  Ti  err  values.  Although  the  concept  of  negative  p02  values  may  appear 
impossible,  it  is  perfectly  reasonable,  provided  that  error  bars  and  uncer¬ 
tainties  are  considered.  Thus  a  value  of  -2  torr  ±  3  torr  legitimately  indi¬ 
cates  an  actual  p02  very  close  to  zero.  Some  oximetry  approaches  ignore 
negative  values  or  bin  them  all  as  zero.  We  accept  any  value  providing  that 
the  error  estimated  for  Ti  is  within  a  specified  range. 

The  most  powerful  aspect  of  FREDOM  is  the  ability  to  follow  the  fate 
of  individual  voxels,  and  thus  we  usually  acquire  at  least  three  baseline  p02 
maps,  followed  by  further  maps  accompanying  interventions,  such  as  hy- 
peroxic  gas  breathing.  Even  under  baseline  conditions,  fluctuations  in  Ti 
are  apparent.  These  may  arise  from  uncertainty  in  Ti,  which  may  be 
reflected  in  Ti  err  or  transient  fluctuations  in  p02. 

Ti  map  allows  us  to  investigate  sequential  maps  and  select  only  those 
pixels  that  consistently  show  small  Ti  errS,  for  example,  white  regions  of 
interest  (ROI).  In  subsequent  maps  we  may  include  only  these  regions, 
which  satisfy  the  inclusion  criteria  for  every  map.  This  allows  us  to  follow 
population  dynamics  based  on  specific  tumor  regions,  avoiding  potential 
anomalies  due  to  varying  numbers  of  “good  quality”  pixels  in  sequential 
images.  Typically,  50-150  pixels  provide  high-quality  data  in  any  given 
map,  but  generally  20-80  pixels  may  be  followed  for  a  series  of  23  maps 
associated  with  interventions. 

The  ROI  tool  also  facilitates  clustering,  for  example,  selection  of  only 
those  tumor  regions,  which  consistently  have  p02  <  10  torr  throughout 
the  baseline  period  (these  may  be  considered  as  chronically  hypoxic,  as 
opposed  to  regions,  which  fluctuate  and  may  only  be  <10  torr  in  some 


Fig.  7.  Maps  showing  (A)  NMR  signal  intensity,  (B)  spin-lattice  relaxation  times  (Ti), 
(C)  estimated  errors  in  Ti,  and  (D)  p02  values.  Black  squares  show  regions  where  signal  to 
noise  or  data  quality  was  so  poor  as  to  not  provide  a  Ti  curve  fit.  Gray  voxels  show  regions  that 
provided  a  curve  fit,  but  where  the  uncertainty  in  the  data  exceeded  the  threshold  criteria  (i.e., 
Ti  err  >  2.5  s  or  Ti  err/Ti  >  50%).  All  colorcd  voxels  (n  =  100)  provided  a  curve  fit  within  the 
acceptance  criteria  in  this  map,  but  only  those  56  voxels  within  the  white  regions  of  interest 
(ROI)  provided  consistently  high-quality  data  throughout  the  sequence  of  eight  maps  acquired 
during  baseline  and  with  oxygen  breathing.  Data  were  obtained  using  the  FREDOM  approach 
from  a  representative  Dunning  prostate  R3327-AT1  tumor  (2  cm^). 
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baseline  maps).  The  ability  to  follow  groups  of  pixels  with  particular  base¬ 
line  characteristics  has  revealed  heterogeneity  in  response  to  many  inter¬ 
ventions — often  those  regions  initially  well  oxygenated  show  rapid  and 
large  response  to  hyperoxic  gas  breathing,  whereas  those  that  are  initially 
poorly  oxygenated  show  little  and  sluggish  response  (e.g.,  Fig.  8).  This  ap¬ 
proach  can  also  reveal  phenomena  such  as  the  “steal”  effect,  whereby  ini¬ 
tially  well-oxygenated  regions  decrease  in  p02,  while  others  increase, 
which  might  appear  as  “no  change,”  when  histogram-based  population 
statistics  are  used. 

Figure  9  shows  that  the  data  quality  is  strongly  related  to  signal  ampli¬ 
tude.  Although  these  data  represent  56  voxels  from  a  single  tumor,  we  have 
previously  shown  similar  data  for  the  alternate  PFC  perflubron.^^^  Given  the 
high  solubility  of  O2  in  PFCs,  there  could  be  a  concern  that  PFCs  act  as 
reservoirs,  perturbing  local  p02.  Fig.  9B  shows  that  there  is  no  correlation 
(r^  <  0.05)  between  signal  amplitude  (viz.  HFB  concentration)  and  p02. 

As  with  any  measurement,  sampling  is  a  critical  issue.  FREDOM  is 
analogous  to  the  Eppendorf  Histograph  in  that  it  samples  multiple  loca¬ 
tions,  which  appear  to  reflect  interstitial  p02.  Histograph  data  suggest  that 
a  minimum  of  100-140  data  points^^^  along  five  tracks^^^  are  required  to  ac¬ 
curately  represent  the  p02  distribution  of  a  tumor,  though  such  criteria  do 
depend  on  tumor  size  and  heterogeneity.  We  applied  Monte  Carlo  simula¬ 
tion  to  assess  the  data  requirements  for  FREDOM.  A  data  set  was 
selected,  which  provided  120  high-quality  data  points,  and  these  data  points 
were  accessed  in  random  order  with  continuous  calculation  of  mean  p02  to 
assess  the  asymptotic  trend  lines  toward  the  actual  mean  p02  (Fig.  10).  It 
appears  that  about  50  data  points  are  required  to  well  represent  the  p02 
distribution  of  the  tumor,  which  is  generally  achievable  in  any  individual 
PO2  map.  However,  the  unique  ability  to  observe  dynamic  changes  in 
PO2  at  multiple  locations  simultaneously  is  the  greatest  strength  of  FRE¬ 
DOM.  Detection  of  heterogeneous  responses  is  useful  even  if  fewer  data 
points  are  examined,  since  each  location  serves  as  its  own  control. 

Validation  of  Measurements.  The  spin-lattice  relaxation  rate  (Ri)  is  de¬ 
termined  for  each  voxel  using  a  three-parameter  fit  of  signal  intensities  to 

yi  =  A[1  -  (1  +  W)  exp(-Ri  x  r)]  (4) 

using  the  Levenberg-Marquardt  least  squares  fitting  protocol.^^^  Typically, 
^100-300  voxels  provide  an  Ri  fit  and  potential  p02  value.  Because  noise 
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A 


Fig.  9.  (A)  Relationship  of  Ri  err  with  signal  amplitude  showing  strong  correlation  (r^  > 
0.69).  (B)  Relationship  of  p02  and  signal  amplitude  showing  lack  of  correlation  (P  <  0.05). 


Fig.  8.  PO2  maps  obtained  using  the  FREDOM  approach  from  the  tumor  shown  in  Fig.  3. 
(A)  Under  baseline  conditions,  those  voxels  (n  =  56)  within  the  white  region  provided 
consistently  reliable  data  during  repeated  measurements.  (B)  Baseline  map  (breathing  air: 
FO2  =  21%):  mean  p02  =  7.2  ±  2.6  (SE)  torr,  median  p02  =  1.3  torr  (range  -7-88  torr). 
(C)  Breathing  oxygen  (FO2  =  100%):  fourth  map  obtained  24-32  min  after  switching  from  air: 
mean  p02  =  47.2  ±  10.7  torr  (p  <  0.0001  compared  with  baseline),  median  p02  =  8.3  torr 
(range  -3-204  torr).  (D)  Breathing  carbogen  (FO2  =  95%):  fifth  map  after  switching  to 
carbogen:  mean  p02  =  43.1  ±  9.2  torr  (p  <  0.0001  compared  with  baseline),  median  p02  =  5.6 
torr  (range  -10-216  torr). 
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Fig.  10.  Monte  Carlo  simulation  of  tumor  oxygenation.  Asymptotic  behavior  shows  that 
tumor  may  be  characterized  by  about  50  data  points.  For  this  large  ATI  tumor  in  a  rat 
breathing  oxygen,  all  estimates  converge  on  the  mean  p02  =  31.7  +  1.0  (SE)  torr. 


itself  may  give  an  apparent  relaxation  curve  (Ri)  fit,  data  are  selected 
within  a  region  of  interest,  having  Ti  err  <  2.5  s  and  the  ratio  Ti  err/Ti  < 
50%.  With  respect  to  respiratory  interventions,  only  those  voxels  that  pro¬ 
vided  consistently  reliable  data  throughout  the  measurements  are  included 
for  further  analysis.  At  37°  and  4.7  T 

p02(torr)  =  [Ri(s-i)  -  0.0835]/0.001876  (5) 

Equation  (5)  provides  p02  in  units  of  torr.  The  literature  can  be  compli¬ 
cated  by  use  of  various  units,  and  it  may  be  instructive  to  provide  conver¬ 
sion  factors.  We  favor  torr  (1  torr  =  1  mmHg),  since  radiobiologic  hypoxia 
develops  in  the  range  0-15  torr  and  this  is  the  traditional  unit  favored 
by  radiation  biologists  and  oncologists.  Further,  760  torr  =  1  standard  at¬ 
mosphere  (atm.),  and  gases  are  often  quoted  in  %atm.  For  SI  units,  101, 
325  Pa  (or  N/m^)  =  1.01325  x  10^  dynes/cm^  =  1  atm.  Some  investigators 
quote  oxygen  concentrations  in  //M,  often  assuming  that  the  solubility  of 
oxygen  in  water  at  37°  is  1.35  //M/torr.^^  Quoting  concentrations  can  be 
confusing  because  the  solubility  of  oxygen  is  highly  variable  in  solutions. 
FREDOM  calibration  is  based  on  p02  values,  and  as  shown  by  Eqs.  1-3, 
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these  are  directly  related  to  the  solubility  of  oxygen  in  HFB.  However, 
partitioning  of  oxygen  between  the  aqueous  and  PFC  phases  depends  on 
p02,  not  concentrations.  Thus  the  p02  determined  by  FREDOM  will  ac¬ 
curately  reflect  the  ambient  tissue  p02  and  if  the  solubility  of  oxygen  can 
be  estimated  in  the  particular  milieu,  then  [O2]  could  be  calculated. 

Both  systematic  and  random  errors  may  interfere.  Random  errors  may 
be  diminished  by  performing  multiple  repeat  measurements  (provided  that 
the  biologic  system  is  stable).  Systematic  errors  are  more  complex  and 
could  arise  inter  alia  from  erroneous  calibrations,  inappropriate  curve  fit¬ 
ting,  and  temperature  changes.  Appropriate  curve  fitting  may  be  the 
greatest  problem  in  relaxation  analysis.  Provided  signal/noise  >  10  for  the 
most  intense  signals,  we  generally  obtain  excellent  curve  fits.  Each  voxel 
(or  ROI)  will  comprise  HFB  at  a  range  of  p02  values,  creating  a  multiex¬ 
ponential  curve.  However,  modeling  shows  that  the  fit  provides  an  “aver¬ 
age”  value.  For  population  data,  we  sometimes  determine  error-weighted 
means  [E(x/cr2)]/[E(l/cr2)]  in  order  to  exploit  as  much  available  data  as 
possible. 

In  terms  of  absolute  p02  values,  the  most  critical  aspect  is  effective  cali¬ 
bration  curves.  Over  the  years  we  have  achieved  p02-dependent  ^^F  NMR 
relaxation  curves  for  many  and  encountered  potential 

pitfalls.  The  calibration  curve  we  recommend  at  4.7  T  and  37°  is  given  by 
Eq.  (5)  which  was  originally  presented  by  Hunjan  et  al}^^  Briefly,  125  //I 
HFB  was  added  to  each  of  four  gas-tight  NMR  tubes  together  with 
0.5  ml  water  and  saturated  at  37°  by  bubbling  with  carbon  dioxide,  1% 
O2  (balance  N2),  9.8%  O2  (balance  N2),  or  air,  respectively.  Tubes  were 
sealed,  and  the  phantom  was  maintained  at  37°  in  a  water  bath  within  a  coil 
in  the  magnet.  FREDOM  was  applied  using  the  parameters  described 
earlier,  and  the  spin-lattice  relaxation  rates  were  estimated  on  a  voxel- 
by-voxel  basis  using  three-parameter  fit.  Equation  (5)  was  established  using 
linear  regression  analysis  of  amplitude  squared  weighted  mean  values  for 
each  gas.^^^ 

Although  the  relationship  Ri  =  f  (PO2)  is  theoretically  expected  to  be 
linear  and  empirically  found  to  be  so,  we  believe  that  it  is  important  to  use 
calibration  gases  in  the  range  of  physiologic  p02.  We  recommend  purchase 
of  rigorously  calibrated  gases.  We  bubble  gases  for  30  min  and  use  gas-tight 
Wilmad  NMR  tubes,  which  may  be  sealed  with  ground  glass  joints. 
Samples  are  used  within  hours  of  saturation.  As  desired,  multiple  tubes 
may  be  prepared  at  given  p02,  but  since  a  linear  relationship  is  expected, 
it  may  be  equally  appropriate  to  use  additional  gases.  Ti  may  be  deter¬ 
mined  multiple  times  for  each  sample,  and  we  recommend  using  the  pulse 
sequence  to  be  applied  for  in  vivo  investigations.  In  1997,  we  published  a 
calibration  curve 
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Ri(s-i)  =  0.074(±0.003)  +  0.012(±0.0002)p02(%atm.)  (6) 

for  HFB  at  37°  and  4.7  Considering  temperature  sensitivity  gave 

Ri(s-i)  =  0.77(±0.03)  -  (0.00009  ±  0.001)T(°C) 

+  (0.018  ±  0.003)P(%atm.)  -(0.00017  ±  0.00008)TP(°C%  atm.) 

(7) 

More  recently,  we  achieved  Eq.  (5)  as  a  new  calibration  curve  at  37°  and 
4.7  T  using  the  ARDVARC  approach.^^^  This  calibration  curve  is  based 
on  higher-quality  data  and  provides  superior  p02  estimates.  In  particular, 
when  tissue  is  expected  to  be  hypoxic  (excised  tissue),  we  find  fewer 
apparently  negative  p02  values. 

A  major  strength  of  the  FREDOM  approach  is  that  calibration  curves 
remain  valid  between  samples  and  across  experimental  platforms.  Calibra¬ 
tion  curves  obtained  in  vitro  are  valid  in  vivo.  Thus  we  believe  other  inves¬ 
tigators  can  apply  Eq.  (5)  without  the  need  to  establish  their  own 
calibrations,  provided  that  studies  are  undertaken  at  about  4.7  T  and  in 
the  range  30-45°.  As  discussed  in  an  earlier  section  and  Table  II,  a  1°  error 
in  temperature  estimates  will  only  introduce  a  0.13-torr  error  in  p02  esti¬ 
mate,  when  p02  is  about  5  torr.  We  also  note  that  the  calibration  curves 
for  HFB  are  relatively  insensitive  to  magnetic  field. 

We  have  previously  investigated  microbiodistribution  of  HFB  based  on 
Oil  Red  O  stain,  which  indicated  that  HFB  occurs  as  microscopic  droplets 
(1-20  fim)  widely  distributed  across  tumor  tissue.^^^  Significantly,  there 
was  no  evidence  for  formation  of  films,  which  could  act  as  conduction  con¬ 
duits,  causing  oxygen  equilibration.  Occasionally,  an  animal  (viz.  tumor) 
will  move  slightly  during  a  long  time  series  of  measurements.  In  this  case 
the  requirement  of  consistently  high-quality  curve  fits  throughout  the  data 
set  for  individual  voxels  fails  because  the  tissues  have  moved  relative  to  the 
voxel  grid.  Such  motion  is  immediately  apparent  when  examining  the 
images.  There  is  a  choice  of  eliminating  such  a  data  set  or  relaxing  the  ac¬ 
ceptance  criteria  and  examining  p02  population  distribution  without 
spatial  continuity.  It  is  important  to  recognize  that  all  methods  include  a 
degree  of  sampling.  In  some  approaches,  this  involves  the  selective  place¬ 
ment  and  tracking  of  electrodes;  in  others,  the  choice  of  locations  for 
biopsy  and  the  number  of  microscopic  fields  of  view.  Sampling  can  be 
avoided  by  obtaining  global  measurements,  as  commonly  acquired  using 
near-infrared  approaches,  but  this  may  itself  mask  the  fundamental  tumor 
heterogeneity^^^  (see  Chapter  17  in  this  volume). 

We  have  previously  shown  that  HFB  shows  little  macroscopic  redistri¬ 
bution  over  a  period  of  hours. It  does  clear  from  the  tumors  with  a 
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typical  half-life  of  HFB  of  about  600  min,  though  some  tumors  show  essen¬ 
tially  no  detectable  clearance  over  a  period  of  6  Clearance  of  HFB 

precludes  long-term  studies  of  chronic  oxygenation,  unless  further  doses  of 
HFB  are  administered. 

Comparison  of  p02  distributions  using  FREDOM  or  the  Eppendorf 
Histograph  has  shown  close  similarity  in  both  small  and  large  tumors. 
Dynamic  studies  in  several  tumor  types  have  shown  equivalent  behav¬ 
ior  when  assessed  using  polarographic  oxygen  electrodes  or  OxyLite  or 
FOXY  optical  probes.^^’"^^’^^"^  Relative  hypoxia  has  been  compared  with 
the  histologic  reporter  pimonidazole,  revealing  similar  trends  across  tumor 
types. 

Applications  of  FREDOM 

FREDOM  has  been  applied  to  investigations  of  diverse  tumor  types 
(syngeneic  rat  prostate  and  breast  tumors  and  xenograft  human  lym¬ 
phomas)  with  respect  to  growth  and  acute  interventions.^^’^^’^^"^”^^^’^^^  Per¬ 
haps  the  greatest  strength  is  the  ability  to  investigate  regional  dynamic 
changes  in  p02  accompanying  acute  interventions.  Figure  8  shows  changes 
in  regional  p02  in  an  ATI  tumor  accompanying  hyperoxic  gas  breathing. 
As  published  previously,^^^’^^^’^^^  the  areas  that  were  initially  better  oxy¬ 
genated  responded  with  elevated  p02,  whereas  those  relatively  hypoxic 
regions  showed  little  response.  Data  also  may  be  presented  as  histograms 
(Fig.  11),  revealing  significant  differences  between  mean  and  median  p02 
and  hypoxic  fractions  between  small  and  large  tumors,  and  between  the 
slow-  and  fast-growing  sublines  H  and  ATI  of  the  Dunning  prostate 
R3327.  Figure  12  shows  variation  in  global  mean  p02  accompanying  re¬ 
spiratory  challenge  with  oxygen  or  carbogen  and  return  to  baseline 
following  the  intervention.  Because  individual  tumor  regions  may  be  ob¬ 
served,  local  response  may  be  compared  at  identical  locations,  efficiently 
comparing  the  efficacy  of  interventions  (Fig.  13).  Figure  14  shows  the  dif¬ 
ferential  behavior  of  regions  in  undifferentiated  Dunning  prostate  R3327- 
ATl  tumors  versus  highly  differentiated  H  tumors.  In  each  tumor  type 
there  are  both  well  and  poorly  oxygenated  regions  under  baseline  condi¬ 
tions  (see  Fig.  11).  In  response  to  hyperoxic  gas  breathing,  well-oxygenated 
regions  in  both  tumor  types  show  a  rapid  and  significant  elevation  in  p02, 
which  is  reversible  when  inhaled  gas  is  returned  to  air  (Fig.  14).  Poorly  oxy¬ 
genated  regions  in  H  tumors  respond  slowly,  but  ultimately  rise  above  the 
range  of  radiobiologic  hypoxia,  whereas  corresponding  regions  in  ATI 
tumors  do  not.  We  also  have  investigated  vasoactive  drugs"^^  and  vascular 
targeting  agents,^^  and  the  short-term  changes  in  p02  following  irradiation 
also  have  been  examined.^^"^  Perhaps  the  most  significant  results  to  date 
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Fig.  11.  Histograms  of  tumor  oxygenation  pooled  for  six  H  (A)  and  seven  ATI  (B)  small 
Dunning  prostate  R3327  rat  tumors.  The  H  tumors  (219  voxels)  showed  significantly  higher 
baseline  p02  (mean  =  31  ±  2  torr,  median  =  27  torr)  than  the  size-matched  ATI  tumors  (338 
voxels;  mean  =  14  ±  1  torr,  median  =  11  torr;  p  <  0.0001).  With  respect  to  oxygen  challenge 
(bottom),  mean  p02  in  both  groups  increased  significantly  (mean  =  122  ±  7  in  the  H  versus, 
mean  =  58  ±  4  in  the  ATI;  p  <  0.0001).  Arrows  indicate  mean  (x)  and  median  (m)  p02, 
respectively. 


show  that  PO2  measurements  and  detection  of  changes  in  p02  accompany¬ 
ing  interventions  correlate  with  the  efficacy  of  tumor  irradiation. 


Future 

Ultimately,  the  value  of  a  technique  depends  on  its  robustness,  ease  of 
use,  and  widespread  implementation.  To  date,  few  laboratories  had 
adopted  the  FREDOM  approach  because  efficient  investigation  of  HFB 
requires  an  unusual  NMR  pulse  sequence.  With  the  recent  upgrade  of 
our  own  instrumentation  to  the  Varian  Unity  INOVA,  the  software  is 
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Fig.  12.  Dynamic  p02  (mean  ±  SE)  obtained  from  sequential  maps  of  the  ATI  tumor 
shown  in  Fig.  8  with  respect  to  respiratory  challenge.  *,  p  <  0.001;  p  <  0.0001  versus 
baseline. 


Fig.  13.  Correlation  between  maximum  p02  detected  in  each  of  56  voxels  from  the  ATI 
tumor  shown  in  Fig.  8,  when  the  rat  breathed  oxygen  versus  carbogen  (r^  >  0.85). 


now  available  on  this  popular  platform,  facilitating  ready  implementation 
elsewhere.  In  terms  of  research  applications,  it  is  known  that  tumor  tissue 
p02  varies  rapidly  in  response  to  many  acute  interventions,  ranging  from 
irradiation  to  photodynamic  therapy,  and  various  chemotherapies.  We 
foresee  FREDOM  as  a  valuable  tool  for  assessing  the  dynamic  time  course 
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Fig.  14.  Variation  in  p02  with  respiratory  challenge  for  individual  regions  chosen  as 
initially  well  oxygenated  (p02  >  10  torr)  or  hypoxic  (p02  <  5  torr)  from  an  ATI  (dotted  lines) 
and  H  (solid  lines)  tumor,  respectively.  All  the  well-oxygenated  regions  in  the  two  sublines 
increased  significantly  in  p02  in  response  to  oxygen  breathing  {p  <  0.01).  The  hypoxic  regions 
from  the  H  tumor  increased,  whereas  those  in  the  ATI  tumor  did  not. 


of  such  interventions  to  provide  clear  insight  into  the  mode  of  action  of 
therapeutic  approaches  and  aid  in  the  high-throughput  screening  of  new 
drugs,  such  as  vascular  targeting  and  antiangiogenic  agents. 
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TUMOR  PHYSIOLOGIC  RESPONSE  TO  COM BRETASTATIN  A4  PHOSPHATE 

ASSESSED  BY  MRI 


Dawen  Zhao,  M  .D.,  Ph.D.,  Lan  Jiang,  M  .S.,  Eric  W.  Hahn,  Ph.D.,  and  Ralph  P.  Mason,  Ph.D. 

Department  of  Radiology,  University  of  Texas  Southwestern  Medical  Center,  Dallas,  TX 

Purpose:  To  evaluate  the  effect  of  the  vascular  targeting  agent,  combretastatin  A4  phosphate,  on  tumor 
oxygenation  compared  with  vascular  perfusion/permeability. 

Methods  and  Materials:  MRI  oximetry  and  dynamic  contrast-enhanced  (DCE)-MRI  were  used  to  monitor 

tumor  oxygenation  and  perfusion/permeability  in  syngeneic  13762NF  rat  breast  carcinoma. 

Results:  A  significant  drop  was  found  in  the  mean  tumor  pOz  (23  to  9  mm  Hg,/;  <0.05)  within  90  min  after 
treatment  (30  mg/kg  of  combretastatin  A4  phosphate)  and  a  further  decrease  was  observed  at  2  h  (mean  2  mm 
Hg;p  <0.01).  The  initial  changes  in  PO2  in  the  central  and  peripheral  regions  were  parallel,  but  by  24  h  after 
treatment,  a  significant  difference  was  apparent;  the  PO2  in  the  periphery  had  improved  significantly,  and  the 
center  remained  hypoxic.  These  data  are  consistent  with  DCE-MRI,  which  revealed  an  ~70%  decrease  in 
perfusion/permeability  (initial  area  under  signal-intensity  curve)  at  2  h  (p  <0.001).  The  initial  area  under 
signal-intensity  curve  recovered  fully  after  24  h  in  a  thin  peripheral  region,  but  not  in  the  tumor  center. 
Conclusion;  The  response  observed  by  DCE-M  Rl,  indicating  vascular  shutdown,  paralleled  the  PO2  measure- 
ments  as  expected,  but  quantitative  PO2  measurements  are  potentially  important  for  optimizing  the  therapeutic 
combination  of  vascular  targeting  agents  with  radiotherapy.  ©  2005  Elsevier  Inc. 

MRI,  Vascular  targeting  agent.  Vasculature,  Oxygenation,  Breast  tumor. 


INTRODUCTION 

T umor  growth,  survival,  and  metastasis  depend  critically  on 
the  development  of  new  blood  vessels  (1).  Therefore,  ex¬ 
tensive  research  has  focused  on  developing  strategies  to 
attack  the  tumor  vasculature  (1,  2).  T ubulin-binding  agents 
(e.g.,  combretastatin  A4  phosphate  [CA4P]  and  ZD6126) 
represent  one  kind  of  vascular  targeting  agent  (VTA)  (3,  4). 
Promising  preclinical  studies  have  shown  that  such  agents 
selectively  cause  tumor  vascular  shutdown  and  subse¬ 
quently  trigger  a  cascade  of  tumor  cell  death  in  experimen¬ 
tal  tumors  (4, 5).  A  Ithough  massive  necrosis  can  be  induced, 
tumors  usually  regrow  from  a  thin  viable  rim.  Thus,  a 
combination  of  VTAs  with  additional  conventional  thera¬ 
peutic  approaches  will  be  required  (6,  7).  Several  studies 
involving  the  combination  of  VTAs  with  radiotherapy  (8- 
11)  or  chemotherapeutic  agents  (12)  have  shown  enhanced 
tumor  response. 

To  better  understand  the  mode  of  action,  and  hence, 
optimize  such  combinations,  in  vivo  imaging  approaches 
have  been  initiated  to  monitor  the  physiologic  changes 


resulting  from  VTA  administration  (13-15).  Dynamic  con¬ 
trast-enhanced  (DCE)-MRI  based  on  the  transport  proper¬ 
ties  of  gadolinium-diethylenetri amine  pentaacetic  acid  (Gd- 
DTPA)  is  the  most  commonly  used  imaging  approach  to 
study  tumor  vascular  perfusion  and  permeability.  DCE-M  Rl 
was  included  as  part  of  the  Phase  I  clinical  trials  of  CA4P 
(16,  17).  The  results  of  preclinical  and  clinical  DCE-MRI 
studies  have  shown  a  reversible  change  in  vascular  perfu¬ 
sion  in  the  tumor  periphery  after  a  single  dose  of  VTA 
(18-21).  For  combination  with  radiotherapy,  measurement 
of  tumor  oxygen  dynamics  will  be  especially  important, 
because  reduced  perfusion  can  induce  hypoxia,  potentially 
modulating  the  radiation  response.  A  number  of  studies 
have  reported  an  improved  response  when  administering 
VTAs  after  radiotherapy,  with  the  enhancement  reduced  or 
lost,  if  VTAs  were  administered  before  radiotherapy,  im¬ 
plying  increased  hypoxia  induced  by  VTAs  (8,  10).  Direct 
measurements  of  tissue  PO2  using  the  Eppendorf  electrode, 
conducted  by  Horsman  et  ai.  (22,  23),  found  increased 
hypoxia  3  h  after  CA  4P  or  ZD 6126  administration.  W  e  have 
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developed  a  method  for  measuring  tumor  oxygenation  and 
dynamics  based  on  nuclear  magnetic  resonance  echo 
planar  imaging  after  direct  intratumoral  injection  of  the 
reporter  molecule  hexafluorobenzene  (HFB)  called  "fluoro¬ 
carbon  relaxometry  using  echo  planar  imaging  for  dynamic 
oxygen  mapping"  (FREDOM  )  (24,  25).  This  technique  pro¬ 
vides  PO2  measurements  at  multiple  specific  locations  si¬ 
multaneously  within  a  tumor  and  reveals  the  dynamic 
changes  at  individual  locations  with  respect  to  interven¬ 
tions.  We  have  previously  evaluated  the  tumor  oxygen 
response  to  varying  interventions  such  as  hyperoxic  gas 
breathing  (26,  27).  We  also  had  an  anecdotal  example  of 
PO2  response  to  a  tumor-selective  infarcting  agent  (28).  We 
have  now  applied  both  DCE-M  Rl  and  FREDOM  to  evalu¬ 
ate  tumor  perfusion/permeability  and  oxygen  dynamics  in 
response  to  CA4P  in  conjunction  with  confirmatory  histo¬ 
logic  examination. 

METHODS  AND  MATERIALS 

Tumor  model 

Rat  mammary  carcinoma  13762NF  was  implanted  syngenei- 
cally  in  a  skin  pedicle  surgically  created  on  the  fore  back  of  Fisher 
344  adult  female  rats  (n  =  25,  -150  g,  Harlan),  as  previously 
described  in  detail  (29).  Of  the  25  rats,  9  were  used  for  DCE-M  Rl, 
10  for  FREDOM,  and  6  for  histologic  study.  The  Institutional 
Animal  Care  and  Use  Committee  approved  the  investigations. 

Drug  preparation  and  dosing 

CA4P  was  provided  by  OXiGENE  (Waltham,  M  A).  CA4P  was 
dissolved  in  0.9%  saline  at  a  concentration  of  30  mg/mL  before 
each  experiment.  A  single  dose  of  30  mg/kg  CA4P  was  chosen  for 
this  study,  because  it  is  considered  a  clinically  relevant  dose  (18). 

MRI  experiments 

When  tumors  reached  -1  cm  diameter  (-0.6  cm^),  M  Rl  was 
performed  using  a  4.7  T  horizontal  bore  magnet  with  a  Varian 
Unity  Inova  system.  Each  rat  was  given  intraperitoneal  ketamine 
hydrochloride  (120  julL;  100  mg/mL,  Aveco,  Fort  Dodge,  lA)  as  a 
relaxant  and  maintained  under  general  anesthesia  (air  and  1% 
isoflurane,  Baxter  International,  Deerfield,  IL).  A  27-gauge  butter¬ 
fly  catheter  (Abbott  Laboratories,  Abbott  Park,  IL)  was  placed 
intraperitoneal  I  y  for  infusion  of  CA4P  or  saline  alone.  For  DCE- 
M  R I,  a  tail  vein  wascatheterized  using  a  second  27-gauge  butterfly 
catheter  for  contrast  agent  administration.  For  oximetry,  hexaflu¬ 
orobenzene  (50  /jii,  Lancaster,  Gainesville,  FL)  was  injected  di¬ 
rectly  into  the  tumor  along  two  or  three  tracks  in  a  single  central 
plane  of  the  tumor,  coronal  to  the  rat's  body  using  a  Hamilton 
syringe  (Reno,  NV)  with  a  custom-made,  fine,  sharp  needle  (32- 
gauge),  as  previously  described  in  detail  (25).  A  tunable  (^H/^®F) 
volume  radiofrequency  (RF)  coil  was  placed  around  the  tumor¬ 
bearing  pedicle.  Each  animal  was  placed  on  its  side  in  the  magnet 
with  no  change  in  position  during  the  whole  study,  so  that  indi¬ 
vidual  regions  could  be  tracked.  A  thermal  blanket  was  used  to 
maintain  body  temperature. 

DCE-MRI 

Nine  tumor-bearing  rats  were  studied  before  CA4P  injection  (n 
=  6)  or  saline  alone  (n  =  3)  and  2  and  24  h  after  treatment.  On 


each  occasion,  a  series  of  T^-weighted  spin  echo  images  (TR  160 
ms,  TE  16  ms,  field  of  view  40  x  40  mm,  matrix  128  x  128,  voxel 
size  2.0  X  0.3  x  0.3  mm,  total  timefor3slices23  s)  was  acquired 
before  and  after  a  bolus  injection  of  Gd-DTPA-BMA  (injection 
within  1  s;  0.1  mmol/kg.  Omniscan,  Amersham  Health,  Princeton, 
NJ)  on  three  2-mm-thick  cross-sections  parallel  to  the  animal.  Data 
were  processed  on  a  voxel-by-voxel  basis  using  software  written 
by  us  using  Interactive  Data  Language  (IDL),  version  5. 3/5. 4 
(Research  Systems,  Boulder,  CO).  For  each  slice,  the  tumor  was 
separated  into  central  and  peripheral  regions.  The  tumor  periphery 
was  taken  to  be  a  1-2-mm-thick  rim  aligned  around  the  whole 
tumor.  Signal  intensity  vs.  time  curves  were  plotted  and  relative 
signal  intensity  changes  (AS/)  of  each  tumor  voxel  were  analyzed 
using  the  equation:  (AS/)  =  (S/^  -  51^^)131^^,  where S/^  refers  to  the 
enhanced  signal  intensity  in  the  voxel  and  S/^  is  defined  as  the 
average  of  the  baseline  images.  The  area  under  the  normalized 
signal  intensity-time  curve  (lAUC)  for  the  first  1.5  min  after 
Gd-DTPA-BMA  injection  was  integrated. 


tumor  oximetry- FREDOM 

A  separate  cohort  of  10  tumors  (7  treated  and  3  controls)  was 
used  for  PO2  measurement.  A  single  4-mm  slice  parallel  to  the  rat 
body  containing  the  strongest  fluorine  signal  was  chosen  for  the 
^^F  M  Rl  PO2  studies.  ^H  and  ^®F  M  R  images  were  acquired  using 
a  spin-echo  sequence.  Overlaying  the  ^^F  M  R  image  on  the  cor¬ 
responding  ^H  image  revealed  the  distribution  of  HFB.  After 
conventional  M  Rl,  tumor  oxygenation  was  estimated  on  the  basis 
of  ^®F  pulse  burst  saturation  recovery  echo  planar  imaging  relax¬ 
ometry  of  the  HFB,  as  previously  described  (24).  This  approach 
provided  PO2  maps  with  1.25  mm  in  plane  resolution  and  6-/xL 
voxel  size  in  6.5  min.  The  spin-lattice  relaxation  rate  [Rl  (s“^)  = 
1/Tl]  was  estimated  on  a  voxel-by-voxel  basis  using  a  three- 
parameter  monoexponential  function.  PO2  was  estimated  using  the 
relationship  PO2  mm  Hg  =  (Rl  -  0.0835)/0.001876  (24).  The 
data  are  presented  in  bins  of  5  mm  Hg,  except  for  the  highest  and 
lowest  bins,  which  were  open  ended.  Before  CA4P  or  saline 
injection,  a  series  of  PO2  maps  was  acquired  with  respect  to 
respiratory  challenge  with  oxygen:  typically,  two  baseline  mea¬ 
surements,  three  with  oxygen  and  four  on  return  to  air.  Immedi¬ 
ately  after  the  last  (fourth)  air  measurement,  CA4P  (30  mg/kg)  or 
saline  (0.15  mL)  was  injected  intraperitoneal ly.  An  additional 
series  of  PO2  maps  was  acquired  after  10,  30,  60,  90,  and  120  min, 
and  finally  another  three  maps  while  breathing  oxygen.  The  24  h 
follow-up  study  comprised  two  measurements  with  air  and  four 
with  oxygen.  The  oxygen  challenge  was  included  to  evaluate 
vascular  function. 


Markers  of  vascular  perfusion  and  endothelium 
Six  animals  were  used  to  study  the  total  blood  vessels  and 
perfused  vessels  before  (n  =  2),  2  h  (n  =  2),  and  24  h  (n  =  2)  after 
CA4P.  The  blue  fluorescent  dye  Hoechst  33342  (Molecular 
Probes,  Eugene,  OR)  was  injected  into  the  tail  vein  of  the  anes¬ 
thetized  rats  at  a  concentration  of  10  mg/kg  in  0.9%  saline  (0.1 
mL),  and  the  tumors  were  excised  1  min  later.  Tumor  specimens 
were  immediately  immersed  in  liquid  nitrogen  and  then  stored  at 
-80°C.  Immediately  after  cryostat  sectioning  (6  jam  thick),  the 
slices  were  imaged  for  Hoechst  33342  under  ultraviolet  wave¬ 
length  (330-380  nm).  Perfused  vessels  were  determined  by  count¬ 
ing  the  total  number  of  structures  stained  by  Hoechst  33342  in  four 
fields  per  section  selected  to  show  high  perfusion  and  calculating 
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the  mean  number  of  vessels  per  millimeter  squared.  On  the  fol¬ 
lowing  day,  the  same  slices,  as  well  as  their  adjacent  slices,  were 
immunostained  for  the  endothelial  marker,  CD 31.  A  primary 
mouse  anti-rat  CD31  monoclonal  antibody  (1:20  dilution,  Serotec, 
Raleigh,  NC)  was  added  and  incubated  for  2  h  at  37°C  in  a  humid 
box.  The  slides  were  incubated  with  Cy3-conjugated  goat  anti¬ 
mouse  secondary  antibody  (1:100  dilution,  Jackson  Immunore- 
search  Laboratories,  West  Grove,  PA)  for  1  h  at  37°C.  After 
mounting  with  Vectorshield  medium  (Vector  Laboratories,  Bur¬ 
lingame,  CA),  the  slides  were  observed  under  red  fluorescence 
(530-550  nm  excitation)  to  detect  anti-CD31  and  then  the  corre¬ 
sponding  Hoechst  33342  again  under  ultraviolet  light.  Image  anal¬ 
ysis  was  performed  using  M  etaview  software  (Universal  Imaging, 
West  Chester,  PA). 


Statistical  analysis 

Statistical  significance  was  assessed  using  analysis  of  variance 
on  the  basis  of  Fisher's  Protected  Least  Significant  Difference 
(StatView,  SAS  Institute,  Cary,  NC)  or  Student's  t  tests. 
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RESULTS 

DC E -MR I  findings 

Baseline  MRI  showed  stable  signal  throughout  the 
tumor  (not  shown).  Immediately  after  a  bolus  injection  of 
Gd-DTPA-BMA,  the  signal  intensity  increased  signifi¬ 
cantly,  reaching  a  peak  after  23-69  s,  and  then  gradually 
decreasing  toward  baseline.  The  maximal  signal  enhance¬ 
ment  averaged  over  the  whole  tumor  in  all  three  slices 
ranged  from  ~52%  to  85%  for  the  nine  tumors. 

A  single  dose  of  CA4P  30  mg/kg  caused  a  dramatic 
decrease  in  signal  enhancement  at  the  2  h  point,  and  a  24  h 
follow-up  image  showed  signal  recovery  in  a  thin  tumor 
rim,  as  shown  for  a  representative  tumor  (Fig.  1).  Distribu¬ 
tion  of  the  lAUC  values  showed  a  distinct  shift  to  the  left  at 
both  2  and  24  h  after  treatment  (Fig.  lb).  The  highly 
enhancing  fraction  (lAUC  >1.5)  at  24  h  after  CA4P  reached 
the  pretreatment  level,  indicating  complete  recovery  in  the 
tumor  periphery  (Fig.  lb).  For  all  six  tumors  treated  with 


Fig.  1.  Dynamic  contrast-enhanced-magnetic  resonance  imaging  (DCE-MRI)  performed  before,  2  h,  and  24  h  after 
treatment  with  CA4P  (30  mg/kg).  Series  of  T^-weighted  contrast-enhanced  images  acquired  at  these  times  with  respect 
to  infusion  of  bolus  injection  of  contrast  agent,  (a)  Normalized  images  of  signal  enhancement  23  s  after  contrast  agent 
administration  shown  for  representative  tumor.  Significantly  less  signal  enhancement  was  observed  for  whole  tumor 
region  2  h  after  treatment.  Full  recovery  is  apparent  in  tumor  rim  24  h  after  treatment,  (b)  Based  on  serial  images 
obtained  with  respect  to  contrast  agent  infusion.  Initial  area  under  signal-intensity  curve  (lAUC)  was  calculated  on 
voxel-by-voxel  basis.  Frequency  distribution  of  lAUC  values  is  shown.  Distributions  at  2  and  24  h  showed  significant 
increases  in  number  of  voxels  with  no  signal  enhancement  (lAUC  <0.05).  However,  frequency  of  highly  enhancing 
voxels  (lAUC  >1.5)  at  24  h  recovered  to  pretreatment  level. 
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CA4P,  a  significant  decrease  in  lAUC  was  observed  for  the 
whole  tumor  region  at  2  h:  the  tumor  peripheral  and  central 
region  was  reduced  by  65%  and  70%,  respectively  (p 
<0.05;  Table  1).  A  complete  recovery  in  I A  DC  for  the 
tumor  periphery  was  seen  24  h  after  treatment,  but  the  tumor 
center  remained  low  (Table  1).  lAUC  frequency  data  from 
these  six  tumors  also  showed  a  significant  Increase  In  the 
number  of  nonenhancing  voxels  (lAUC  <0.05),  and  the 
percentage  of  voxels  with  lAUC  >0.5  decreased  from  55% 
to  5%  2  h  after  administration  of  CA4P  (p  <0.05,  Table  2). 
In  control  tumors  treated  with  saline,  the  lAUC  was  similar 
at  all  points  (Tables  1  and  2). 

Tumor  oximetry— F  RE  DOM 

Overlay  of  on  Image  (Fig.  2a)  confirmed  thatHFB 
was  distributed  In  both  peripheral  and  central  regions  of  a 
central  plane  of  the  tumor.  In  the  series  of  echo  planar 
imaging  relaxation  data  sets,  typically  —40-100  voxels 
provided  an  R1  fit,  and  potential  PO2  value.  Because  noise 
may  give  an  apparent  relaxation  curve  (Rl)  fit,  data  were 
selected  by  applying  thresholds  of  T1  error  <2.5  s  and  T1 
error/Tl  ratio  <50%.  Only  those  voxels,  which  provided 
consistently  reliable  data  throughout  the  time  course,  were 
included  for  additional  analysis.  As  an  example  (Fig.  2b),  25 
voxels  qualified  for  all  these  prerequisites  and  could  be 
followed  through  the  17  measurements  from  baseline  to  2  h 
after  CA4P.  At  24  h,  31  voxels  were  traceable  through  the 
six  measurements  with  respiratory  challenge.  As  expected, 
on  the  basis  of  the  distribution  of  HFB  (Fig.  2a),  each  PO2 
map  (Fig.  2b)  comprised  two  distinct  groups  of  voxels, 
representing  the  peripheral  and  central  location  of  FIFB, 
respectively.  Variations  in  the  mean  PO2  with  respect  to 
intervention  for  these  two  groups  of  voxels  are  presented 
separately  in  Fig.  2c.  Before  CA4P  administration,  the  mean 
PO2  was  45  mm  Hg  in  the  tumor  periphery  and  36  mm  Hg 
in  the  center  during  air  breathing.  Both  regions  responded 
significantly  to  oxygen  breathing  (p  <0.05),  and  this  was 
reversed  on  return  to  air  breathing.  A  significant  decrease  in 
PO2  was  detected  as  early  as  30  min  after  CA4P  adminis¬ 
tration  in  the  tumor  periphery  and  by  60  min  in  the  center  (p 


<0.05).  The  decline  continued  with  the  lowest  PO2  values 
(13  mm  Hg  and  10  mm  Hg)  at  the  2-h  point.  At  this  time, 
oxygen  challenge  no  longer  produced  an  increase  in  PO2  for 
either  region.  At  24  h  later,  the  peripheral  PO2  had  increased 
significantly  (24  mm  Hg),  but  was  still  significantly  lower 
than  the  pretreatment  value  (p  <0.05).  After  24  h,  PO2  again 
responded  significantly  to  breathing  oxygen  in  the  tumor 
periphery  (p  <0.01).  However,  PO2  in  the  central  regions  (4 
mm  Hg)  was  even  lower  at  24  h  than  at  2  h  and  did  not 
respond  to  oxygen  breathing.  Histograms  of  pooled  voxels 
assessed  for  the  seven  tumors  showed  significantly  de¬ 
creased  PO2  at  1,  2,  and  24  h  after  treatment  (p  <0.01).  The 
hypoxic  fraction  <5  mm  Hg  (HF5)  increased  significantly 
from  a  pretreatment  value  of  21%  to  36%  at  1  h  and  68%  at 
2  h,  but  was  only  34%  at  24  h  (Fig.  3).  For  the  seven  treated 
tumors,  two  tumors  showed  a  significant  decrease  in  PO2  as 
early  as  30  min  after  treatment  (p  <0.05).  A  significant 
decrease  in  the  mean  PO2  for  the  group  was  found  at  90  min 
after  CA4P  (23  ±  5  vs.  9  ±  3  mm  Hg;  p  <0.05)  and  a 
further  decrease  at  2  h  (2  ±  2  mm  Hg;  Table  3).  The  mean 
PO2  increased  significantly  to  15  ±  4  mm  Hg  at  24  h 
compared  with  2  h  (p  <0.05).  For  the  control  tumors,  no 
significant  differences  in  PO2  were  found  among  any  of  the 
measurements  (before  and  2  or  24  h  after). 

Histologic  examination  and  immunohistochemistry 

On  the  basis  of  the  perfusion  assessed  by  distribution  of 
Hoechst  33342,  a  significant  decrease  in  the  number  of 
perfused  vessels  was  evident  at  2  h  (21  ±  6/mm^  vs.  145  ± 
23/mm^  pre-CA4P,  p  <0.01),  followed  by  recovered  per¬ 
fusion  at  24  h  (114  ±  19/mm^).  Comparisons  of  perfused 
vessels  marked  by  Hoechst  33342  and  with  the  total  vessels 
labeled  by  anti -CD 31  before  and  2  and  24  h  after  CA4P  are 
shown  in  Fig.  4. 

DISCUSSION 

In  common  with  previous  reports  in  the  literature  (15, 
18),  CA4P  caused  a  significant  reduction  in  tumor  perfusion 
within  2  h.  DCE-M  Rl  showed  that  the  lAUC  in  the  tumor 


Table  1.  Normalized  lAUC  by  DCE-M  Rl 


Mean  lAUC 

Group 

Baseline 

2  h 

24  h 

CA4P  (n  =  6)  (30  mg/kg) 

Periphery 

0.88  ±  0.06 

0.31  ±  0.08* 

0.86 

±  0.06 

Center 

0.44  ±  0.11+ 

0.13  ±  0.08* 

0.12 

±  0.05*+* 

Control  (n  =  3)  (saline) 

Periphery 

1.09  ±  0.04 

0.93  ±  0.13 

0.86 

±  0.05 

Center 

0.52  ±  0.19+ 

0.46  ±  0.05+ 

0.44 

±  0.03+ 

Abbreviations:  lAUC  =  initial  area  under  signal-intensity  curve;  DCE  =  dynamic  contrast 
enhanced:  CA4P  =  combretastatin  A4  phosphate. 

*  p  <0.05  from  baseline. 

^  p  <0.05  from  periphery. 

*  p  <0.05  from  control. 
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Table  2.  Comparison  of  weakly  (lAUC  <0.05)  and  strongly 
(lAUC  >0.05)  responding  voxels  with  respect  to  DCE-MRI  and 
drug  administration 


Total  voxels  (%) 

Group 

lAUC 

Baseline 

2  h 

24  h 

CA4P  (n  =  6)  (30  mg/kg) 

<0.05 

16 

54*t 

41*t 

>0.5 

55 

32*t 

Control  (n  =  3)  (saline) 

<0.05 

14 

15 

19 

>0.5 

62 

58 

51 

Individual  voxels  (total  =  10,508  in  CA4P  vs.  5376  in  control) 
categorized  according  to  lAUC. 

Abbreviations  as  in  Table  1. 

*  p  <0.01  from  baseline, 
p  <0.05  from  control. 

periphery  was  reduced  70%  2  h  after  a  single  dose  of  CA4P 
(30  mg/kg),  followed  by  a  full  recovery  by  24  h.  However, 
the  drop  in  the  lA  UC  of  the  tumor  center  was  not  reversible 
(Table  1).  The  lAUC  reflects  tumor  blood  flow,  vascular 
permeability,  and  a  fraction  of  interstitial  space.  Histologic 
data  using  the  perfusion  marker,  Hoechst  33342,  confirmed 
a  significant  decrease  in  perfused  vessels  2  h  after  CA4P, 
with  recovery  evident  at  24  h  (Fig.  4). 

Vascular  shutdown  induced  by  VTA  is  expected  to  de¬ 
crease  tumor  PO2.  Horsman  etal.  (23)  recently  reported  that 
tumor  PO2,  measured  by  Eppendorf  electrode,  was  signifi¬ 
cantly  decreased  3  h  after  ZD6126.  In  line  with  that  study, 
M  Rl  oximetry  showed  that  CA4P  caused  a  significant 
decrease  in  PO2  within  90  min  and  a  further  reduction 
within  2  h.  As  for  perfusion,  PO2  recovered  toward  pretreat¬ 
ment  levels  after  24  h  in  the  tumor  periphery  (Fig.  2C). 

Although  VTAs  cause  vascular  shutdown  and  lead  to 
extensive  hemorrhagic  tumor  necrosis,  rapid  tumor  re¬ 
growth  occurs  from  the  surviving  rim.  Many  preclinical 
studies  have  shown  that  a  single  dose  of  VTA  is  associated 
with  very  little  tumor  growth  delay  (4,  5,  23).  Thus,  a 
combination  with  other  therapeutic  modalities  (e.g.,  radio¬ 
therapy  and  conventional  chemotherapy)  will  be  needed  to 
achieve  complete  control.  The  rationale  for  such  combina¬ 
tions  Is  that  VTAs  target  the  poorly  perfused  and  relatively 
hypoxic  regions  of  the  central  tumor,  and  conventional 
therapy  is  applied  to  attack  the  relatively  better  oxygenated, 
highly  proliferating  rim.  A  critical  issue  regarding  such 
combinations  is  the  timing  and  sequencing  of  the  agents. 
Thus,  a  primary  goal  of  this  study  was  to  monitor  longitu¬ 
dinal  PO2  changes  after  CA4P  to  facilitate  and  potentially 
optimize  such  combinations.  Others  have  demonstrated  en¬ 
hanced  therapeutic  efficacy  by  combining  CA4P,  ZD6126, 
or  5, 6-di methyl xanthenone-4-acetic  acid  (DMXAA)  with 
radiotherapy  (8-11),  radioimmunotherapy  (30),  or  chemo¬ 
therapy  (12).  In  the  study  of  combined  ZD6126  with  radio¬ 
therapy,  Siemann  and  Rojiani  (8)  reported  increased  tumor 
cell  kill  by  administering  drug  immediately  (30  min  to  1  h) 
after  a  single  radiation  dose  compared  with  administering 
the  drug  24  h  or  1  h  before.  This  observation  coincides  with 
our  present  findings  that  tumor  PO2  declined  1  h  after  drug 
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administration  and  had  not  recovered  fully  at  24  h  (Figs.  2C 
and  3). 

Oxygen  challenge  before  and  after  treatment  was  used  to 
compare  vascular  function  In  this  study.  In  common  with 
our  previous  observations  (31,  32),  oxygen  breathing  sig¬ 
nificantly  increased  tumor  oxygenation  in  mammary  carci¬ 
noma  13762NF  under  control  conditions  (Fig.  2  and  Table 
3).  The  PO2  response  to  CA4P  was  equally  effective  in  the 
peripheral  and  central  regions  through  2  h  after  administra¬ 
tion  (Fig.  2C).  PO2  decreased  in  both  regions  within  30  min 
and  continued  to  decrease  for  2  h.  Differential  behaviorwith 
respect  to  oxygen  challenge  between  the  peripheral  and 
central  tumor  regions  was  observed  only  after  24  h.  At  24  h, 
PO2  in  the  tumor  center  was  even  lower  than  the  pretreat¬ 
ment  baseline  or  the  2-h  level  and  showed  no  response  to 
oxygen  inhalation  (Fig.  2).  However,  oxygen  breathing  pro¬ 
duced  a  significant  increase  in  PO2  of  the  peripheral  region, 
even  though  the  PO2  was  significantly  lower  than  before 
treatment  (Fig.  2).  These  observations  provide  further  evi¬ 
dence  that  CA4P  at  30  mg/kg  induces  irreversible  vascular 
damage  in  the  tumor  center,  and  peripheral  vessels  survived 
and  were  functional  by  24  h.  Onewould  expect  that  hypoxia 
modifiers  (e.g.,  oxygen  breathing),  if  given  24  h  after  CA4P, 
might  improve  tumor  radiosensitivity. 

Hypoxia  in  solid  tumors  has  been  widely  recognized  as  a 
potent  factor  that  leads  to  resistance  to  radiotherapy  and 
some  anticancer  drugs  (33).  Recently,  Increasing  evidence 
has  shown  that  tumor  malignant  progression  may  be  asso¬ 
ciated  with  a  hypoxic  microenvironment  (34,  35).  Given  the 
importance  of  oxygen,  many  techniques  for  monitoring  PO2 
have  been  developed  (25,  36).  Although  each  method  has 
specific  attributes,  many  are  highly  invasive  and  impractical 
for  longitudinal  studies  of  specific  regions  of  interest.  The 
Eppendorf  has  been  considered  the  gold  standard  for  PO2 
measurement,  but  it  is  not  suitable  for  a  longitudinal  study 
of  the  type  reported  here.  In  vivo,  proton  M  Rl  (e.g.,  DCE- 
MRI  and  blood  oxygen  level-dependent  MRI)  provides  a 
noninvasive  approach  to  assess  tumor  vasculature,  particu¬ 
larly  useful  in  response  to  interventions,  but  neither  DCE 
nor  blood  oxygen  level -dependent  M  Rl  provides  a  straight¬ 
forward  correlation  with  PO2  (37,  38).  FREDOM  not  only 
provides  PO2  values  simultaneously  at  multiple  specific 
locations  within  a  tumor,  but  also  reveals  dynamic  changes 
at  individual  locations  with  respect  to  interventions. 

Central  necrosis  develops  in  13762NF  tumors  even  at  a 
small  size  and  becomes  substantial  when  tumors  become 
bigger,  which  has  been  confirmed  previously  by  histologic 
examination  (32).  In  some  of  the  tumors  used  in  this  study, 
a  low-intensity  central  region  was  observed  in  Ti-weighted 
images,  in  which  little  enhancement  was  noted  after  infu¬ 
sion  of  the  contrast  agent.  After  treatment  with  CA4P,  the 
area  of  these  regions  of  low-signal  intensity  increased 
slightly  (data  not  shown).  A  similar  observation  was  re¬ 
ported  by  Beauregard  etal.  (39)  in  human  colon  carcinoma 
xenograft.  In  common  with  our  previous  study  of  13762NF 
tumors  (31,  32),  we  found  considerable  intratumoral  heter¬ 
ogeneity  in  the  distribution  of  PO2  values,  ranging  from 
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Fig.  2.  (a)  Distribution  of  hexafluorobenzene  (HFB)  in  representative  tumor  (No.  4  in  Table  3).  Overlay  of  ^^F  signal 
density  on  anatomic  image  indicates  HFB  in  both  peripheral  and  central  regions,  (b)  PO2  maps  obtained  from  same 
tumor  comprise  two  separated  groups  of  voxels  that  correspond  to  locations  of  H  FB  on  anatomic  image;  25  individual 
voxels  were  traceable  from  pretreated  baseline  to  2  h  after  CA4P  with  oxygen  breathing;  and  31  voxels  could  be 
followed  24  h  after  treatment.  Significant  decrease  in  PO2  was  evident  for  all  individual  voxels  after  CA4P.  PO2  did  not 
respond  to  oxygen  inhalation  after  2  h;  24-h  maps  showed  improved  PO2  and  significant  response  to  oxygen  breathing 
in  peripheral  region  (right),  but  not  in  central  region  (left),  (c)  Mean  PO2  curves  shown  for  peripheral  (black  squares) 
and  central  (white  circles)  voxels  of  tumor.  *p  <0.05  from  baseline  air,  <0.05  from  24-h  air,  "^p  <0.05  from 
periphery. 
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Fig.  3.  pOj  histograms  based  on  reiative  numbers  of  voxeis  in  seven  tumors  showing  ieft  shift  after  CA4P.  Ordinate 
ranges  (e.g.,  15  refers  to  10  mm  Hg  <  pOj  <  15  mm  Hg).  x  =  mean;  m  =  median. 


hypoxic  to  well  oxygenated.  Compared  with  our  previous 
data  on  13762NF  tumors,  we  found  the  average  baseline 
PO2  values  were  lower  in  the  present  study  (mean  25  mm 
Hg  vs.  63  mm  Hg).  This  difference  likely  resulted  from  an 
altered  anesthesia  protocol,  because  we  now  use  a  baseline 
of  air/isoflurane,  instead  of  33%  oxygen  in  nitrous  oxide 
with  metafane  (32).  However,  oxygen  challenge  produced  a 
remarkably  similar  extent  of  Increase  in  PO2  values  ob¬ 
served  In  both  studies  (72  mm  Hg  vs.  73  mm  Hg). 

Previous  studies  showed  little  macroscopic  redistribution 


of  HFB  during  a  period  of  hours  (40),  but  it  did  clear  from 
the  tumors  with  a  typical  half-life  about  600  min  (41).  The 
clearance  of  HFB  normally  precludes  chronic  long-term 
studies  of  oxygenation,  unless  additional  doses  of  HFB  are 
administered.  For  all  three  control  tumors  in  this  study, 
additional  HFB  (50  ^aL)  was  required  for  the  24-h  follow-up 
study.  However,  no  additional  HFB  was  required  for  studies 
24  h  after  CA4P,  presumably  owing  to  vascular  shutdown. 

These  results  provide  the  first  insight  into  regional  tumor 
oxygen  dynamics  in  response  to  CA4P  in  a  syngeneic  rat 


Table  3.  Tumor  oxygen  dynamics  assessed  by  ^®F  M  Rl  with  respect  to  CA4P  treatment 

PO2  (mm  Hg) 


Time  course  post  CA4P 


Group 

Case 

No. 

Size 

(cm^) 

Baseline 

air 

30  min 

60  min 

90  min 

2  h 

Oxygen 

24  h  base 

Oxygen 

CA4P  (30  mg/kg) 

1 

2.0 

36  ±  1 

NA 

NA 

NA 

2  ±  3* * 

5  ± 

2* 

NA 

NA 

(n  =  1) 

2 

0.9 

12  ±  2 

NA 

NA 

NA 

-3  ±  2* 

-3  ± 

2* 

NA 

NA 

3 

1.2 

13  ±  0 

35  ±  11 

11  ±  3 

9  ±  4 

1  ±  2* 

2  ± 

2* 

5  ±  0* 

30  ±  6+ 

4 

0.7 

43  ±  5 

29  ±  3* 

21  ±  3* 

20  ±  2* 

13  ±  1* 

16  ± 

5* 

20  ±  2* 

48  ±  7+ 

5 

0.4 

12  ±  2 

5  ±  3* 

9  ±  3 

4  +  2* 

2  ±  1* 

7  ± 

3* 

NA 

NA 

6 

0.7 

15  ±  1 

19  ±  3 

9  ±  2* 

4  ±  2* 

1  ±  2* 

0  ± 

3* 

12  ±  1 

73  ±  12^ 

7 

0.3 

28  ±  6 

29  ±  2 

16  ±  4 

6  ±  2* 

1  ±  1* 

7  ± 

2* 

23  ±  2* 

79  ±  6+ 

Mean 

0.9  ±  0.2 

23  ±  5 

23  ±  13 

13  ±  5 

9  ±  3* 

2  ±  2* 

5  ± 

3* 

15  ±  4 

58  ±  11+ 

Saline  (n  =  3) 

8 

0.6 

34  ±  2 

47  ±  5 

40  ±  5 

36  ±  10 

31  ±  8 

88  ± 

10* 

37  ±  0 

86  ±  5+ 

9 

0.2 

47  ±  2 

62  ±  8 

65  ±  6 

64  ±5 

61  ±  4 

114  ± 

8* 

67  ±  1 

159  ±  7+ 

10 

0.6 

13  ±  1 

16  ±  6 

10  ±  4 

9  ±  4 

12  ±  6 

35  ± 

10* 

14  ±  1 

78  ±  15+ 

M  ean 

0.5  ±  0.1 

31  ±  10 

42  ±  14 

38  ±  16 

36  ±  16 

35  ±  14 

79  ± 

23* 

39  ±  15 

108  ±  26+ 

Abbreviations:  NA  =  no  measurement:  other  abbreviations  as  in  Table  1. 
Data  presented  as  mean  ±  SE. 

*  p  <0.05  from  baseline. 

p  <0.05  from  24  h  air  (24  h  after  injection  of  CA4P). 
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Fig.  4.  Perfusion  marker  Hoechst  staining  before,  2  h,  and  24  h  after  CA4P  (left).  Vascular  endothelium  of  same  field 
immunostained  by  anti-CD31  (red,  middle).  Good  match  (right)  between  Hoechst  and  anti-CD31-stained  vascular 
endothelium  found  in  pretreated  tumor.  Significant  reduction  in  perfused  vessels  detected  2  h  after  treatment  followed 
by  a  recovery  at  24  h.  Bar  =  100  jam. 


tumor.  A  distinct  similarity  was  noted  between  the  results  of 
the  PO2  measurements  and  the  more  traditional  DCE,  but 
the  quantitative  PO2  values  provide  the  potential  for  exploit¬ 
ing  synergy  with  other  oxygen-dependent  therapies.  The 
observations  also  demonstrate  the  value  of  FRE DOM  in 
assessing  dynamic  changes  in  regional  tumor  PO2  in  vivo  in 


response  to  intervention.  We  believe  that  dynamic  measure¬ 
ments  are  particularly  valuable  for  understanding  the  mode 
of  action  of  therapeutic  response  to  VTAs.  Most  signifi¬ 
cantly,  these  measurements  lay  a  foundation  to  optimize  the 
timing  of  combination  therapy  involving  fractionated  radio¬ 
therapy  and  multiple  doses  of  VTAs. 
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ABSTRACT  Bioluminescence  imaging  (BLI)  has 
found  significant  use  in  evaluating  long-term  cancer 
therapy  in  small  animals.  We  have  now  tested  the 
feasibility  of  using  BLI  to  assess  acute  effects  of  the 
vascular  disrupting  agent  combretastatin  A4  phosphate 
(CA4P)  on  lucif erase-expressing  MDA-MB-231  human 
breast  tumor  cells  growing  as  xenografts  in  mice. 
Following  administration  of  luciferin  substrate,  there  is 
a  rapid  increase  in  light  emission  reaching  a  maximum 
after  about  6  min,  which  gradually  decreases  over  the 
following  20  min.  The  kinetics  of  light  emission  are 
highly  reproducible;  however,  following  i.p.  administra¬ 
tion  of  CA4P  (120  mg/kg),  the  detected  light  emission 
was  decreased  between  50  and  90%,  and  time  to 
maximum  was  significantly  delayed.  Twenty-four  hours 
later,  there  was  some  recovery  of  light  emission  follow¬ 
ing  further  administration  of  luciferin  substrate.  Com¬ 
parison  with  dynamic  contrast-enhanced  MRI  based  on 
the  paramagnetic  contrast  agent  Omniscan  showed 
comparable  changes  in  the  tumors  consistent  with  the 
previous  literature.  Histology  also  confirmed  shutdown 
of  tumor  vascular  perfusion.  We  believe  this  finding 
provides  an  important  novel  application  for  BLI  that 
could  have  widespread  application  in  screening  novel 
therapeutics  expected  to  cause  acute  vascular  changes 
in  tumors. — Zhao,  D.,  Richer,  E.,  Antich,  P.  P.,  Mason, 
R.  P.  Antivascular  effects  of  combretastatin  A4  phos¬ 
phate  in  breast  cancer  xenograft  assessed  using  dy¬ 
namic  bioluminescence  imaging  and  confirmed  by 
MRI.  FASEB J.  22,  2445-2451  (2008) 

Key  Words:  breast  tumor  *  vascular  targeting  agent  *  vascu¬ 
lar  disrupting  agent  *  pharmacodynamics 


Bioluminescence  imaging  (BLI)  has  found  a  major 
role  in  small  animal  research  (1-3).  For  tumor  cells 
transfected  to  constitutively  express  luciferase,  there 
are  numerous  reports  examining  tumor  growth  and 
metastatic  spread  (4,  5).  BLI  has  also  found  widespread 
use  in  examining  changes  in  tumor  growth  over  a 
period  of  many  weeks  with  diverse  therapies  (6-10). 
BLI  requires  administration  of  luciferin  substrate, 
which  may  be  achieved  by  multiple  routes,  including 
intravenous  (i.v.),  intraperitoneal  (i.p.),  or  subcutane¬ 


ous  (s.c.)  (11,  12).  Luciferin  is  then  carried  throughout 
the  vasculature  and  has  been  shown  to  readily  perme¬ 
ate  every  region  of  the  body  (9,  13,  14).  Given  the 
importance  of  vascular  transport,  it  occurred  to  us  that 
the  measurement  of  the  light-emitting  dynamics  would 
be  related  to  vascular  delivery  of  the  substrate.  Thus, 
any  agent  causing  major  acute  affects  on  tumor  vascu¬ 
lature  could  influence  the  light  emission  kinetics.  In 
particular,  it  appeared  that  vascular  targeting  agents 
(VTAs;  also  referred  to  as  vascular  disrupting  agents), 
such  as  combretastatin,  could  be  assessed  based  on  the 
dynamics  of  light  emission  detected  by  BLI. 

Tumor  growth,  survival,  and  metastasis  depend  criti¬ 
cally  on  the  development  of  new  blood  vessels  (15). 
Therefore,  extensive  research  has  focused  on  develop¬ 
ing  strategies  to  attack  tumor  vasculature  (16-19). 
Tubulin-binding  agents  {e.g.,  combretastatin  A-4-phos- 
phate  (CA4P)  and  ZD6126)  represent  one  kind  of  VTA 
(17).  Promising  preclinical  studies  have  shown  that 
such  agents  selectively  cause  tumor  vascular  shutdown 
and  subsequently  trigger  a  cascade  of  tumor  cell  death 
in  experimental  tumors  (20).  Although  massive  necro¬ 
sis  can  be  induced,  tumors  usually  regrow  from  a  thin 
viable  rim.  To  better  understand  the  mode  of  action 
and,  hence,  facilitate  optimized  therapeutic  combina¬ 
tions,  in  vivo  imaging  approaches  have  been  applied  to 
monitor  physiological  changes  resulting  from  VTA 
administration  (21-23).  Dynamic  contrast-enhanced 
(DCE)  MRI  based  on  the  transport  properties  of  the 
small  paramagnetic  contrast  agent  gadolinium-DTPA 
(Gd-DTPA)  is  the  most  commonly  used  imaging  ap¬ 
proach  to  study  tumor  vascular  perfusion  and  perme¬ 
ability.  DCE  MRI  was  included  as  part  of  the  Phase  I 
clinical  trials  of  CA4P  (24,  25).  Results  of  preclinical 
and  clinical  DCE  MRI  studies  have  shown  a  reversible 
change  in  vascular  perfusion  in  the  tumor  periphery 
following  a  single  dose  of  VTA  (26-29). 

The  VTA  CA4P  causes  tumor  vascular  shutdown, 
inducing  massive  cell  death.  We  recently  showed  acute 
hypoxiation  within  90  min  following  CA4P  administra- 
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tion  to  rats  bearing  syngeneic  breast  13762NF  tumors 
using  MRI  (30).  Rapid  vascular  shutdown  in  tumors 
after  administration  of  CA4P  to  animals  and  patients 
has  also  been  observed  by  many  other  imaging  modal¬ 
ities,  including  positron  emission  tomography  based  on 
the  distribution  of  water  (31),  DCE  MRI  (24,  32, 
33),  DCE  computed  tomography  (CT)  (34),  ^^E  MRI  of 
tumor  oxygenation  (30),  laser  Doppler  flowmetry  (35), 
radiolabeled  iodoantipyrine  uptake  (27),  near  infrared 
spectroscopy  (36),  interstitial  fluid  pressure  (35),  and 
intravital  microscopy  (37).  By  comparison,  BLI  is  par¬ 
ticularly  inexpensive  and  easy  to  apply  in  animal  mod¬ 
els,  and  we  believe  it  could  be  an  effective  screening 
tool  for  evaluation  and  comparison  of  vascular  target¬ 
ing  agents  as  well  as  long-term  tumor  growth.  BLI  is 
exceedingly  sensitive  with  the  capability  of  detecting 
subpalpable  tumor  volumes. 

We  have  now  explored  the  ability  of  planar  BLI,  a 
widely  available  modality,  to  investigate  the  acute  ef¬ 
fects  of  CA4P  on  human  breast  MDA-MB-23I  xenograft 
tumors  and  provide  correlates  with  the  more  tradi¬ 
tional  MRI  approach  and  histology. 


MATERIALS  AND  METHODS 

Tumor  model 

Human  mammary  MDA-MB-23I  cells  were  infected  with  a 
lentivirus  containing  a  hrefly  luciferase  reporter,  and  highly 
expressing  stable  clones  were  isolated.  Tumors  were  induced 
by  injecting  10®  cells  subcutaneously  in  the  flank  of  13 
athymic  nude  mice  (BALB/c  nu/nu;  Harlan,  Indianapolis, 
IN,  USA).  Tumors  were  allowed  to  grow  to  ~6  mm  diameter 
(~I00  mm^)  and  then  selected  for  BLI  or  MRI. 

In  vivo  BLI 

Tumor-bearing  mice  were  anesthetized  (isoflurane/ Og  in  an 
induction  chamber;  isoflurane  from  Baxter  International 
Inc.,  Deerfield,  IL,  USA)  and  a  solution  of  D-luciferin  (450 
mg/kg  in  PBS  in  a  total  volume  of  250  |jl1;  Biosynthesis, 
Naperville,  IL,  USA)  was  administered  s.c.  in  the  neck  region. 
Anesthesia  was  maintained  with  isoflurane  (1%)  in  oxygen 
(I  dm^/min)  and  a  series  of  light  images  (60  s  each)  was 
acquired  using  a  single  reference  camera  of  our  multicamera 
Light  Emission  Tomography  System  over  a  period  of  20-30 
min  and  the  light  intensity-time  curves  evaluated. 

The  imaging  system  uses  CCD  (charge-coupled  device) 
cameras  selected  for  their  high  and  nearly  constant  sensitivity 
over  the  full  range  of  wavelengths  commonly  used  in  optical 
imaging,  from  blue  to  near  infrared.  The  CCD  is  a  noncolor, 
back-illuminated,  full-frame  image  sensor  with  512  X  512 
pixels  [(SITe)  SI-032AB  CCD;  Scientihc  Imaging  Technolo¬ 
gies  Inc.,  Tigard,  OR,  USA].  The  quantum  efficiency  of  the 
CCD  is  greater  than  85%  from  400-750  nm,  and  remains 
above  50%  up  to  900  nm.  The  CCD  has  a  pixel  size  of  24  X 
24  fjim,  providing  a  large  well  capacity  of  350,000  e“,  with  a 
sensitivity  of  2.6  |jLV/e~,  low  dark  current  (20  pA/cm^  at 
20°C),  and  low  readout  noise  (5  e~  RMS),  providing  a 
dynamic  range  of  75,000.  The  CCD  is  cooled  to  -40°C, 
reducing  the  dark  current  signal  to  less  than  O.I  e~/pixel/s. 
The  large  dynamic  range  of  the  detector  is  coupled  with  a  16 
bit  analog  to  digital  converter,  allowing  quantitative  detection 


of  both  high  and  very  low  signals  simultaneously.  The  CCD  is 
incorporated  in  a  self-contained,  cooled  camera  equipped 
with  electronic  circuitry  and  fast  optics  (25  mm  focal  length, 
f/0.95).  Each  camera  is  calibrated  using  a  low-intensity, 
diffuse,  flat-held  source  that  gives  a  known  radiance  (adjusted 
typically  at  3.0X10“^  W/cmVsr).  The  light  source  is  period¬ 
ically  checked  for  uniformity  using  a  NIST-traceable  research 
radiometer  (model  IL  1700;  International  Light,  Inc.  New- 
buryport,  MA,  USA).  By  imaging  this  source,  the  digital  units 
provided  by  the  camera  digitizer  can  be  converted  directly 
into  absolute  physical  units  (W/ cvo? / sr  or  photons/ s/ cvo? / sr) . 

Saline  control  or  combretastatin  A4P  (CA4P;  120  mg/kg; 
OXiGENE,  Inc.  Waltham,  MA,  USA)  in  saline  (120  |jil)  was 
injected  i.p.  immediately  after  baseline  BLI  and  the  mouse 
allowed  to  wake  up.  Two  hours  later,  the  BLI  time  course  was 
repeated  based  on  new  injections  of  luciferin.  In  some  cases, 
mice  were  examined  again  24  h  after  CA4P  administration. 

In  vivo  MRI 

A  second  cohort  of  nude  mice  (7z=3)  with  size-matched 
tumors  was  studied  by  MRI  using  a  4.7  T  horizontal  bore 
magnet  with  a  Varian  INOVA  Unity  system  (Palo  Alto,  CA, 
USA).  Each  mouse  was  maintained  under  general  anesthesia 
(air  and  1%  isohurane).  A  27  G  butterhy  (Abbott  Laborato¬ 
ries,  Abbott  Park,  IL,  USA)  was  placed  in  a  tail  vein  for 
contrast  agent  administration.  Pertinent  image  slice  positions 
were  based  on  fast  scout  images.  Tl-weighted  (TR=200  ms; 
TE=I5  ms;  slice  thickness  =1.5  mm;  FOV=  50X50;  in  plane 
resolution  390  |jLm)  and  corresponding  T2-weighted 
(TR=I500  ms;  TE  =  80  ms)  spin-echo  multislice  axial  images 
were  acquired.  For  DCE  MRI,  a  series  of  3  contiguous 
Tl-weighted  images  (TR=70  ms;  TE  =  I2  ms;  total  acquisition 
time  10  s  with  same  voxel  dimensions  as  above)  was  acquired 
before  and  after  i.v.  bolus  injection  of  the  contrast  agent 
Gd-DTPA-BMA  (O.I  mmol/kg  body  weight;  Omniscan^^, 
Amersham  Health  Inc.,  Princeton,  NJ,  USA)  through  a  tail 
vein  catheter.  DCE  MRI  was  performed  before  and  2  h  after 
CA4P  infusion  (120  mg/kg,  i.p.). 

Data  were  processed  on  a  voxel-by-voxel  basis  using  soft¬ 
ware  written  by  us  using  IDL  5. 3/5. 4  (Research  Systems, 
Boulder,  CO,  USA).  For  each  slice,  the  tumor  was  separated 
into  regions  of  center  and  periphery,  respectively.  The  tumor 
periphery  was  taken  to  be  a  1-2  mm  thick  rim  aligned  around 
the  whole  tumor.  Signal  intensity  vs.  time  curves  were  plotted 
and  relative  signal  intensity  changes  (A 57)  of  each  tumor 
voxel  were  analyzed  using  the  equation  A5/  =  {SI^  —54)/ 54, 
where  54  is  enhanced  signal  intensity  in  the  voxel,  and  54  is 
the  average  of  the  baseline  images.  The  area  under  the 
normalized  signal  intensity-time  curve  (lAUC)  was  integrated 
for  the  hrst  60  s  after  Gd-DTPA-BMA  injection. 

Immunohistochemistry 

Two  hours  after  saline  or  CA4P  injection,  the  blue  fluorescent 
dye  Hoechst  33342  (10  mg/kg,  Molecular  Probes,  Eugene, 
OR,  USA)  was  injected  into  the  tail  vein  of  anesthetized  mice, 
and  the  tumors  were  excised  I  min  later.  Tumor  specimens 
were  immediately  immersed  in  liquid  nitrogen  and  then 
stored  at  —  80°C.  A  series  of  6  |jLm  frozen  sections  from  several 
regions  of  each  tumor  was  immunostained  for  luciferase 
(Serotec  Inc.,  Raleigh,  NC,  USA).  Slices  were  imaged  for 
Hoechst  33342  under  ultraviolet  (UV)  wavelengths  (330-380 
nm) .  Perfused  vessels  were  determined  by  counting  the  total 
number  of  structures  stained  by  Hoechst  33342  in  4  helds  per 
section  selected  to  show  high  perfusion  and  calculating  the 
mean  number  of  vessels  per  square  millimeter.  On  the 
following  day,  the  same  slices,  as  well  as  the  adjacent  ones, 
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Figure  1.  Control  BLI  dynamic  measurement.  A)  Staining  for  luciferase  in  an 
MDA-MB-231-luc  tumor  tissue  section  showed  extensive  fluorescent  mAh  signal 
predominantly  located  in  cytoplasm  of  many  tumor  cells.  B)  Bioluminescent  image 
from  a  representative  MDA-MB-231-luc  tumor  (140  mm^)  acquired  in  1  min,  4  min 
after  luciferin  injection  (450  mg/kg  s.c.).  Q  Bioluminescent  image  from  the  same 
tumor  following  repeat  s.c.  luciferin  injection  2  h  after  saline  injection  (0.15  ml 
i.p.).  B)  The  BLI  intensity  was  integrated  over  the  whole  tumor,  and  the  time 
course  was  essentially  identical  for  the  2  investigations  conducted  2  h  apart  with 
respect  to  saline  administration,  as  a  control.  Maximum  signal  intensity  was 
observed  at  6  min  after  s.c.  luciferin  administration. 


Time  (min) 


were  immunostained  for  the  endothelial  marker,  CD31.  A 
primary  rat  anti-mouse  CD31  monoclonal  antibody  (1:20 
dilution;  Serotec)  was  added  and  incubated  overnight  at  4°C 
in  a  humid  box.  Slides  were  incubated  with  Cy3-conjugated 
goat  anti-rat  secondary  antibody  (1:100  dilution;  Jackson 
Immunoresearch  Laboratories,  West  Grove,  PA,  USA)  for  1  h 
at  37°C.  After  mounting  with  Vectorshield®  medium  (Vector 
Laboratories,  Burlingame,  CA,  USA),  the  slides  were  observed 
under  red  fluorescence  (530-550  nm  excitation)  to  detect 
anti-CD31  and  then  the  corresponding  Hoechst  33342  again 
under  UV  light.  Image  analysis  was  performed  using  Metaview 
software  (Universal  Imaging  Corp.,  West  Chester,  PA,  USA).  For 
luciferase  staining,  monoclonal  mouse  antiluciferase  mAh 
(1:150;  Serotec)  and  fluorescein  iso  thiocyanate-conjugated  goat 
anti-mouse  secondary  antibody  Qackson)  were  used. 

Statistical  analysis 

Statistical  signihcance  was  assessed  using  an  ANOVA  on  the 
basis  of  Fisher’s  protected  least  signihcant  difference  (PLSD; 
Statview;  SAS  Institute  Inc.,  Cary,  NC,  USA)  or  Student’s  t 
tests. 


RESULTS 

Histology  of  MDA-MB-231-luc  cells  showed  extensive 
expression  of  luciferase  detected  by  antiluc  mAb  in  all 
tumor  cells  (Fig.  1).  BLI  showed  intense  light  emission 
from  tumors  following  s.c.  administration  of  luciferin 


reaching  a  maximum  intensity  ~6  min  postadministra¬ 
tion  followed  by  a  gradual  decline  over  the  next  20  min. 
Repeat  measurements  following  administration  of  a 
second  dose  of  luciferin  2  h  after  injecting  saline  i.p. 
showed  highly  reproducible  results  (Fig.  1,  Table  1)  in 
terms  of  light  distribution,  maximum  light  intensity, 
and  time  to  maximum  light  emission.  Five  tumors  were 
examined  again  after  24  h,  and  light  emission  re¬ 
mained  highly  consistent,  with  mean  maximum  inten¬ 
sity  =  97  ±  6%  (range  49-130%).  By  contrast,  repeat 
BLI  2  h  after  treatment  with  CA4P  showed  a  signifi¬ 
cantly  lower  light  emission  (peak  intensity  2-1  OX 
lower)  and  delayed  peak  emission  (Fig.  2;  Table  1). 
Three  mice  were  examined  again  24  h  post-CA4P,  and 
signal  was  again  significantly  lower  than  baseline  (mean 
41  ±5,  range  21-66%),  though  with  some  recovery 
compared  with  2  h.  Histology  showed  a  well-developed, 
well-perfused  vasculature  in  these  tumors,  but  perfu¬ 
sion  was  essentially  eliminated  2  h  after  administration 
of  CA4P  i.p.  (Fig.  3). 

Signal  enhancement  observed  by  DCF  MRI  was  signif¬ 
icantly  less  in  all  3  tumors  with  a  mean  decrease  of  —70% 
in  perfusion/permeability  2  h  after  CA4P  (P<0.001,  Ta¬ 
ble  2,  Fig.  4) .  Little  change  was  observed  in  perfusion  of 
the  femoral  artery  (data  not  shown).  Muscle  in  1  mouse 
indicated  a  significant  reduction  in  lAUC,  but  others 
showed  no  significant  change  (Table  2). 


TABLE  1.  MDA-MB-23 1-luc  mammary  carcinoma  response  to  CA4P  observed  by  BLI 


Group 

Max  RLU  change  (%) 

Time  delay  (min) 

AUC  change  (%) 

Signal  change  at  7  min 

(%) 

Saline  {n=6) 

-16  ±19  (-48  to  3) 

0.1  ±  1.2  (0to2) 

11  ±  20  (-20  to 

31) 

-18  ±  25  (-55  to 

8) 

CA4P  (7^=10) 

-70  ±  15*  (-53  to  -92) 

+7.1  ±5.9*  (Otoll) 

-72  ±  16**  (-42  to 

-93) 

-76  ±  17*  (-97  to  - 

■53) 

Comparison  of  light  emission  at  baseline  and  2  h  after  administering  saline  (control)  or  120  mg/kg  CA4P.  Values  are  presented  as  mean 
±  SD  (range).  Changes  are  relative  to  baseline.  AUC,  area  under  curve  for  20  min;  RLU,  relative  luminescence  unit.  *P  <  0.05;  <  0.01. 
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Figure  2.  Tumor  response  to  CA4P  monitored  by  planar  BLI.  A  second  MDA-MB-231-luc  tumor  (130  mm^)  was  monitored 
sequentially  following  saline  or  CA4P  (120  mg/kg)  i.p.  infusion.  Each  image  was  acquired  in  1  min,  4  min  after  s.c.  luciferin 
injection.  In  contrast  to  saline,  CA4P  treatment  caused  a  significant  decrease  in  BLI  signal  intensity  2  h  after  treatment,  which 
remained  lower  24  h  later.  A)  Baseline  control,  5)  2  h  after  G\4P,  Q  BLI  intensity  curves  showing  intense  signal  pretreatment 
(solid  squares)  with  99%  less  signal  2  h  after  CA4P  (open  circles). 


DISCUSSION 

Following  administration  of  the  vascular  targeting 
agent  CA4P,  BLI  showed  significantly  decreased  and 
delayed  signal  in  the  human  mammary  MDA-MB-231 
carcinoma.  We  interpret  this  finding  as  due  to  induced 
vascular  shutdown.  Dynamic  contrast-enhanced  MRI 
confirmed  decreased  tumor  perfusion,  which  was  fur¬ 
ther  confirmed  by  histology.  Both  BLI  and  MRI  are 
noninvasive  and  depict  the  same  physiological  effects, 
but  BLI  is  much  cheaper  and  offers  a  high  throughput 
method  for  evaluating  novel  drugs  and  drug  combina¬ 
tions  and  scheduling. 

BLI  is  a  very  sensitive  technique  revealing  transfected 
cells  at  even  subpalpable  volumes.  Signal  is  attenuated 
at  depth,  but  the  technique  can  be  applied  routinely  in 
both  s.c.  models  and  for  monitoring  orthotopic  tumor 
development  (3-8,  12).  Several  groups  have  shown 


robust  correlations  between  detected  light  and  tumor 
volume  assessed  by  calipers  or  MRI  (5,  12,  38,  39).  In 
our  early  work  (12),  we  found  strong  correlations 
(r>0.8)  between  peak  light  signal,  initial  area  under 
the  curve,  or  signal  intensity  at  10  min  after  luciferin 
infusion  and  tumor  volume  following  i.p.  administra¬ 
tion  of  150  mg/ kg  D-luciferin.  However,  individual 
measurements  showed  discrepancies,  and  consecutive 
observations  could  vary  by  as  much  as  60%.  Overall 
estimates  at  95%  confidence  intervals  suggested  er¬ 
ror  <  20%.  We  later  established  better  reproducibility 
based  on  450  mg/kg  luciferin  administered  s.c.  It  is 
important  to  note  that  light-emission  kinetics  depend 
on  tumor  location.  Thus,  for  HeLa-luc  cells  growing 
subcutaneously,  we  found  maximum  light-emission 
around  10  (median)  and  12.7  (mean)  min  following 
i.p.  administration,  whereas  in  the  current  study,  using 
s.c.  mammary  tumors  and  s.c.  administration  of  lucife- 
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Figure  3.  Immunohistochemical  study  of  tumor  vascular  response  to  CA4P.  Perfusion  marker  Hoechst  33342  staining  pre-  and 
2  h  post-CA4P  (120  mg/kg).  Vascular  endothelium  of  the  same  field  was  immunostained  by  anti-CD31  (red).  A  good  match 
between  Hoechst  and  anti-CD31  stained  vascular  endothelium  was  found  in  the  pretreated  tumor  (overlay).  Two  hours  after 
treatment,  significant  reduction  in  perfused  vessels  was  detected.  Scale  bar  =100  |jLm. 
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TABLE  2.  DCE  MRI  data 


Tumor  ID 

Tumor  tissue 

Normal  muscle 

Pretreatment 

CA4P  2  h 

Reduction 

(%) 

Pretreatment 

CA4P  2  h 

Reduction 

(%) 

1 

0.74  ±  0.06 

0.49  ±  0.08* 

34 

0.76  ±  0.04 

0.67  ±  0.03 

12 

2 

1.65  ±  0.09 

0.48  ±  0.09* 

71 

1.17  ±  0.04 

0.73  ±  0.04* 

41 

3 

1.42  ±  0.08 

0.30  ±  0.07* 

79 

1.07  ±  0.03 

0.79  ±  0.02 

28 

Mean 

1.27  ±  0.27 

0.42  ±  0.06** 

67  ±  14 

1.00  ±  0.12 

0.73  ±  0.03 

27  ±  8 

lAUCgo  (initial  area  under  signal-time  curve  of  first  60  s)  assessed  across  regions  of  interest  corresponding  to  tumor  and  contralateral  thigh 
muscle  over  3  contiguous  MR  image  slices.  Values  are  means  ±  se.  <  0.05,  <  0.001  vs.  pretreatment. 


rin,  maximum  light-emission  occurred  after  —7.5  min. 
In  the  current  study,  BLI  repeated  2  h  after  saline 
injection  showed  light-emission  kinetics  highly  consis¬ 
tent  with  baseline  in  terms  of  maximum  light  emission, 
time  to  maximum  light  emission,  integrated  light,  and 
signal  observed  at  7  min  (Table  1;  Fig.  1).  In  5  mice, 
BLI  measurements  were  again  repeated  24  h  later,  and 
light  emission  remained  highly  consistent  (not  signifi¬ 
cantly  different  from  baseline  or  2-h  time  points) .  Two 
hours  after  CA4P  administration,  light-emission  kinet¬ 
ics  were  significantly  altered  with  both  delayed  emission 
and  reduced  light  intensity  (Table  1).  Light-emission 
kinetics  remained  low  24  h  later,  when  luciferin  was 
readministered  to  a  subgroup  of  mice. 

Changes  in  vascular  perfusion  were  confirmed  by 
histology.  These  tumors  show  extensive  vasculature  as 
revealed  by  anti-CD31  staining,  which  is  also  extensively 
perfused  (Fig.  3).  Two  hours  after  CA4P  administra¬ 
tion,  vasculature  was  detectable,  but  now  there  was  no 
evidence  of  perfusion. 

Dynamic  contrast-enhanced  MRI  has  been  applied  pre¬ 
viously  both  in  animal  models  and  the  clinic  to  many 
tumor  types  with  respect  to  combretastatin  administration 
(24,  30,  32).  The  changes  observed  here  (Fig.  4;  Table  2) 
are  consistent  with  previous  observations  (23,  30) .  Initial 
area  under  the  curve  reflects  tumor  blood  flow,  vascular 
permeability,  and  the  involved  fraction  of  interstitial 
space.  For  the  initial  control  study,  the  signal  enhance¬ 
ment  was  generally  considerably  greater  in  tumor  than 
contralateral  thigh  muscle  (Table  2;  Fig.  4F).  Following 
combretastatin  administration,  the  lAUC  decreased  sig¬ 
nificantly  in  each  tumor  with  a  mean  of  —70%  in  tumors 
(Table  2)  but  was  not  significantly  changed  in  2  of  3 
muscles.  One  mouse  did  show  decreased  perfusion  of 
normal  muscle  after  combretastatin,  as  also  reported 
previously  by  Tozer  et  al  (23)  in  rats. 

BLI  is  more  traditionally  used  to  assess  tumor  growth 
and  metastatic  spread  rather  than  physiological  effects. 
BLI  enables  detection  of  subpalpable  volumes  and 
deep-seated  tumors  in  mice  although  this  capability  has 
not  been  exploited  in  the  current  study.  BLI  is  nonin- 
vasive  but  does  require  administration  of  luciferin 
substrate.  Thus,  sampling  tumor  vasculature  is  only 
achieved  at  discrete  time  points  following  administra¬ 
tion  of  the  substrate.  The  requirement  for  a  reporter 
substrate  is  common  to  other  methods  of  interrogation 
also,  such  as  DCE  MRI  or  DCE  CT,  which  require 


administration  of  paramagnetic  or  radio  opaque  con¬ 
trast  agents  (24,  30,  32,  34).  Radionuclide  approaches 
can  be  applied  similarly  (31),  although  they  are  more 
often  used  with  autoradiography  postmortem.  Mea¬ 
surements  based  on  Laser  Doppler  flowmetry  (35), 
pressure  probes  (35),  and  NIR  (36)  may  allow  contin¬ 
uous  sampling.  Oximetry  based  on  ^^F  MRI  of  hexaflu- 
orobenzene  allowed  essentially  continuous  mapping  of 
pOg  changes  with  a  6.5  minute  time  resolution  (30). 

One  disadvantage  of  the  BLI  approach  to  assessing 
tumor  vasculature  is  the  need  for  luciferase-expressing 
cells.  However,  numerous  cell  lines  are  now  widely  avail¬ 
able  and  effective  stable  transfection  is  readily  achievable 
though  clonal  selection  could  lead  to  differences  between 
luciferase-expressing  and  parental  tumor  lines. 

While  combretastatin  (or  Zybrestat)  is  now  in  Phase 
II/III  clinical  trials  for  anaplastic  thyroid  cancer,  stud¬ 
ies  will  continue  to  be  required  in  preclinical  settings  to 
optimize  combinations  with  other  chemotherapeutic 
drugs  or  radiation  in  diverse  tumor  types.  We  believe 
BLI  provides  an  optimal  approach  for  examining  acute 
pharmacodynamics  together  with  the  potential  for 
long-term  chronic  assessment  of  tumor  control  or 
growth.  This  approach  could  be  equally  applicable  to 
other  VTAs,  which  cause  significant  acute  changes  in 
blood  flow,  such  as  ZD6126  or  5, 6-dime thylxanthe- 
none-4-acetic  acid  (19). 

There  is  increasing  realization  that  drugs  may  behave 
differently  with  various  tumor  types  and  sites  of  implan¬ 
tation.  The  observations  here  are  consistent  with  a 
previous  report  on  the  same  tumor  type,  which  showed 
extensive  central  necrosis  24  h  after  combretastatin 
administration  while  leaving  a  viable  peripheral  rim 
(20).  BLI  will  facilitate  effective  screening  of  many 
tumor  types  implanted  subcutaneously  or  orthotopi- 
cally  and  with  respect  to  monotherapy  or  combinations. 
We  believe  this  is  the  first  report  of  the  use  of  BLI  to 
monitor  acute  effects  of  a  drug  on  tumor  vasculature 
and  represents  an  additional  application  of  this  new 
technology. 

This  work  was  supported  in  part  by  the  U.S.  Department  of 
Defense  Breast  Cancer  Initiative  Idea  Awards  to  D.Z.  (DAMD- 
170310363)  and  to  P.P.A.  (W81XWH-04-1-0551),  the  Na¬ 
tional  Cancer  Institute  Cancer  Imaging  Program  (pre-ICMIC 
CA86354  and  U24  CAl  26608),  and  Simmons  Cancer  Center. 
NMR  experiments  were  performed  at  the  Advanced  Imaging 
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Figure  4.  DCE  MRI  monitoring  of  tumor  response  to  CA4P.  A-C)  Conventional  MR  images  of  a  nude  mouse  with  an 
MDA-MB-231-luc  mammary  tumor  (Tumor  2  in  Table  2):  Tl-weighted  (A),  T2-weighted  {B),  and  Tl-weighted  (C)  contrast- 
enhanced  (Gd-DTPA-BMA).  D,  E)  Dynamic  contrast-enhanced  MRI  was  performed  in  the  mouse  before  {D)  and  2  h  after  (£) 
CA4P  i.p.  injection.  Color  scale  map  =  normalized  contrast  enhancement  30  s  after  a  bolus  injection  of  Gd-DTPA-BMA  is 
overlaid  on  the  Tl-weighted  images.  Significantly  decreased  signal  enhancement,  compared  with  pretreatment,  was  observed  2  h 
after  i.p.  injection  of  GA4P.  F)  Ratio  of  tumorimuscle  signal  enhancement  following  infusion  of  Gd-DTPA-BMA  in  the  single 
image  slice  shown  in  D  and  E.  Solid  symbols  =  baseline;  open  symbols  =  2  h  after  GA4P. 


Research  Center,  a  National  Institutes  of  Health  Biological 
Threat  Reduction  Program  facility  (P41-RR02584) .  Com- 
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